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An Investigation of Separated Flows— 


Part II: 


Flow in the Cavity and Heat 


Transter 


A. F. CHARWAT,* C. F. DEWEY, JR.,** J. N. ROOS,? ano J. A. HITZtt 
University of California at Los Angeles 


Summary 


The first portion of this paper describes studies of the internal 
structure of the separated flow in a notch at a free-stream Mach 
number of 3. Observations include: flow visualization, spark- 
schlieren pictures of the fluctuations of the free shear layer, and 
studies of the diffusion of heat from sources placed in the sepa- 
rated region. The second part describes measurements of local 
heat transfer to the wall 

The external Mach number, the length-to-depth ratio of the 
cavity, the ratio of the oncoming boundary layer thickness to the 
notch depth (in the turbulent flow region), the thermal-to- 
momentum thickness ratio of the boundary layer and, finally, 
the geometry of the internal boundary of the separated region 
are varied as systematically as possible. On the basis of these 
observations, a simple model of the flow in and the heat transfer 
across the separated region is formulated. 
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Symbols 


local heat-transfer coefficient 

depth of notch 

length of notch 

critical closure length 

Mach number in the free stream 

pressure 

heat flux per unit area to wall 

Reynolds number 

Stanton number 

distance to virtual origin of the boundary layer 

temperature 

adiabatic wall recovery temperature 

streamwise coordinate (along notch floor from separa- 
tion corner 

coordinate normal to free stream (origin at notch wall) 

boundary-layer thickness just ahead of notch 

thickness of the free shear layer at reattachment 


Subscripts 
reference flat plate 
free-stream conditions 
wall condition 
conditions at separation point 


conditions at recompression point 


Introduction 


HE {INFORMATION on heat transfer to and across 
separated regions is difficult to interpret because 
of differences in test conditions and in the geometry of 


the models—i.e., notches,*'''* wakes, and separation 
steps in supersonic!*'* and subsonic‘ flow, axisymmetric 
spike-induced compression wakes, and various other 
complex separations”: '*. General results indicate a 
drop in the heat transfer coefficient h after separation 
as compared to the value ahead of separation, and an 
increase of h near and behind recompression. The 
average over the whole cavity floor varies from 5() to 
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120 percent of that for an equivalent length of attached 
flow. 

A theoretical analysis was proposed by Chapman? 
and applied by Larson.'* This model treats the cavity 
as an isothermal “‘dead air’’ sink to which heat is trans- 
ferred through the shear layer. The entire problem 
is assumed to lie in the transfer properties of this layer, 
which are estimated from an extension of an asymptotic 
solution of free jet mixing. 

The outstanding characteristic of this model is that 
it indicates a pronounced change of heat transfer to the 
cavity upon transition from laminar to turbulent flow 
in the shear layer. This is a consequence of the change 
of the shear layer profile, the sole transfer-rate-con- 
trolling element in the analysis. The absence of such 
an effect would constitute, in our opinion, a strong 
argument against the validity of the model. 

Chapman’s model disregards the resistance between 
the inner wall and the fluid in the cavity. There have 
been attempts to formulate the problem with more 
emphasis on the nature of the internal flow; Cole” 
considered the convection of heat around a vortex, 
but it is difficult to apply his analysis to a real case. 
Carlson?! modified Chapman’s work by assuming a 
velocity profile (in the form of a polynomial) both over 
a free shear layer and in the recirculated flow along the 
floor of the cavity. This assumption is entirely arbi- 
trary and not substantiated by data. 

None of the analyses mentioned above take into 
account the possibility that the flow may be unsteady 
and may involve a periodic ‘‘filling’’ and ‘“‘emptying”’ 
of the cavity with fresh fluid. If this were the case, 
the transfer across the shear layer would no longer 
control the transfer of heat to or across the cavity. 
Furthermore, it would not be necessary that any of 
the aforementioned three mechanisms dominate over 
the others. 

The foregoing arguments guided the planning of the 
present work. We attempted to investigate system- 
atically many aspects of the flow and obtain insight 
into its true nature by a critical comparison of the re- 
sults. 

It was shown in Part I of this series,! in which the 
pressure field was discussed, that the flow will jump 
across a cutout to form an “‘open’’ cavity when the 
length-to-step-height ratio of the notch is less than a 
‘critical’? value, which is about 11 for turbulent flow in 
the Mach number range between 1.5 and 4. The 
similarity parameter which governs the pressure dis- 
tribution in the “open” cavity regime is the ratio of 
cavity length to critical closure length, L/L,,.. There 
is a large decrease in L,, upon transition from laminar 
to turbulent flow in the shear layer. Consequently, a 
cavity which is “open” in laminar flow may ‘‘close”’ 
upon transition to turbulence. 

The present paper describes investigations of the 
internal flow in an “open” cavity. The results com- 
plement and explain the dynamic pressure traverses 
discussed in Part I and prove that there is a periodic 
inflow into the cavity. Also, measurements of local 
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heat transfer to the wall are presented and discussed in 
terms of a model based on the observed unsteady struc- 
ture of the flow. 


Test Equipment and Measurements 


Studies were made of the internal flow in a notch 
in the test section floor at a Mach number of 2.9. The 
boundary layer at separation was fully turbulent. 
The equivalent flat-plate Reynolds number was 1.5 X 
10° and the ratio of the boundary layer thickness to the 
step height was approximately 0.4.* 

The heat-transfer measurements were conducted in 
a specially designed wind tunnel. The test section 
was located in a flat plate replacing the plane of sym- 
metry of a two-dimensional wedge-flow expansion 
nozzle. The nozzle was not contoured, but the ex- 
pansion angle was small enough to ensure satisfactorily 
small gradients. The streamwise pressure variation 
per inch was always below 4 percent of the static 
pressure at the separation point. 

Fig. 1 gives the configuration of the test block, 
which was made up of mutually insulated stainless 
The insulation consisted of a 
Matching 


steel modular elements. 
0.020-in. air gap and a 0.005-in. mica sheet. 
cutouts in the modules formed an internal coolant 
passage. Each module was provided with an iron 
constantan thermocouple junction located 0.003 in. 
below the test surface at the center of the element. 

The wind tunnel operated at atmospheric stagnation 
pressure and stagnation temperatures between 50 and 
600°F. Water was used in the internal coolant passage 
(heat flow to the wall); a few special tests (heat flux to 
the air) were run with steam circulated through the 
internal passage and atmospheric stagnation air tem- 
perature. The thermal resistance of the modules was 
such that the temperature drop across each could be 
maintained, typically, at around 100°. Under this 
condition instrument calibration errors were below 5 
percent at JJ = 2 and below a maximum 10 percent 
at 7 = 3.5. Note that the use of the ratio h/h,, where 
h; is measured independently for each test, tends to 
eliminate such errors. Heat conduction between the 
modules and losses to the sidewall were neglected. An 
estimate of the possible error arising from this source 
indicates that the results are correct within 3 percent 
on the average, and within 10 percent under the worst 
conditions, that is, in regions of maximum streamwise 
temperature gradient. The modular construction was 
continued five modules beyond the test section proper. 

All reference quantities (such as Mach number, 
stagnation temperature, Reynolds number, etc.) were 
determined in the free flow at the station corresponding 
to the separation point (beginning of the notch). The 
Reynolds numbers were calculated from measured 
boundary-layer thicknesses (equivalent flat-plate Reyn- 
olds numbers) by use of the relation: 


6/S = 0.37Re-"5 = 0.37(UyS/m)-™ (1) 


* This arrangement is the same as that used in reference 1 
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Fic. 1. Photograph of the test section (exploded view) and schematic diagrams of its construction 


The boundary layer was turbulent, and could be cor- 
related to a good accuracy with a !/;7-power velocity 


profile. 


Flow Visualization 
A) The Flow Near the Wall 


An oil film smeared on the floor of the cavity or a 
suspension of lampblack in kerosene injected into it 
(Fig. 2) were used to show the gross nature of the in- 
ternal flow. It was found that all but the shortest 
notches (L/H < 1.5) contain a vortex at the separation 
corner and a distinct stagnation line separating it from 
This vortex opposes the 
‘diameter’ re- 
Another vortex 


the remainder of the flow. 
mean circulation in the cavity. Its 
mains equal to 7 regardless of L. 

exists at the recompression corner. 

equal to H. There is no stagnation line separating it 
from the central section of the cavity because the flow 
in that section is in the same (upstream) direction near 
the floor. Nevertheless, the complete lack of lamp- 
black in the corner, as compared to the central section, 
bears evidence to the much larger velocity in this com- 


Its diameter is 


pression ‘‘vortex.”’ 

The fact that the velocities in the central zone are 
low and unsteady was demonstrated by introducing 
small bearing balls or drops of mercury into the open 


cavity near recompression. They rolled upstream 
slowly (at rates of the order of one inch per second), up 
to the separation vortex. (Sometimes they would 
then creep over to the corners driven by the three- 
dimensional flow pattern to be discussed.) The diam- 
eter of the balls had no noticeable effect on their ac- 
celeration as long as they were less than about '/; of 
The motion of the balls was uneven 


The structure of the flow in an open 


the notch depth. 
and unsteady. 
cavity suggested by these observations, and confirmed 
by more elaborate measurements, is shown  sche- 
matically in Fig. 18. 

Fig. 2 shows an accumulation of the oil near the 
corners in longer notches, which indicates a _ three- 
dimensional shedding of the separation vortex. The 
location of the disturbance suggests a sidewall effect, 
which can be explained by considering the centrifugal 
force balance on a streamtube in the sidewall boundary 
layer, under the influence of the pressure distribution in 
A special model with fences separating out 
Fig. 2c shows 


the notch. 
the sidewall boundary layer was tested. 
a run in which oil is collected along a well-defined and 
uniformly two-dimensional line (vortex attached to 
the fences). However, this two-dimensional flow was 
found to be extremely sensitive to the alignment of 
the separation or recompression edges, to upstream 
nonuniformities, and even to transient disturbances. 
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Similar three-dimensional patterns were observed 
by all investigators who looked for them in notches,’ 
free wakes,® and steps.’ The question was raised 
whether two-dimensional cavity flow is fundamentally 
possible; our results indicate that it is, but only under 


extremely carefully controlled conditions. 


(b) Mass Exchange and Unsteadiness of the Flow 

Any detectable gas (ammonia, smoke) injected into 
the free flow upstream of the notch was found to pene- 
trate into the cavity almost instantaneously. When 
injected into the cavity, the gas was immediately de- 
tected in the boundary layer downstream of recom- 
pression. When injection was interrupted, the con- 
centration of the indicator in the cavity became zero 
very rapidly. There seems to be no doubt that there 
exists a strong exchange of mass between the cavity 
and the external flow. 

Unsteady flow is demonstrated by the series of spark- 
schlieren pictures shown in Fig. 3. These are taken 
from a study to be published elsewhere,’ which indicates 
that the shear layer forms a stationary space wave 
with wavelength approximately equal to H for all notch 
lengths and frequency proportional to Mo. This photo- 
graphic technique did not yield clear results for the 
(turbulent) thick boundary-layer case, Fig. 3c; the 
pictures show that there is some wave radiation present, 
but it is too diffuse and random to have a definable 
space structure. On the basis of preliminary results, 
there is a limiting ratio of the thickness of the oncoming 
boundary layer to notch depth beyond which no 
radiation is observed. For laminar flow it was found 
to be 6/H = 0.3. A study of turbulent flow is under 
way. 

The results of which Fig. 3 is a sample show that the 


OlL_ FILM ACCUMULATION 
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Stroudhal number based on notch depth (fly //), 
rather than (fU)/L),”° is constant. As Z increases the 
space-wave structure changes abruptly from one regime 
to the next which contains one more wave (of constant 
length). 

Note the way in which the root of the shock segment 
at recompression in Fig. 3b measurably penetrates into 
the cavity. Such pictures indicate again that this is 
a periodic flow in and out of the cavity past the re- 
compression step, rather than a potential wave motion 
like that on a stretched membrane. 


Heat Diffusion Studies 


For further study a method which is directly related 
to the fluid motion was chosen. The arrangement is 
A thin Nichrome wire was stretched 
A thermo- 


shown in Fig. 4. 
across the cavity and heated electrically. 
couple junction was rotated at a fixed radius around the 
heat source, and the temperature difference (with and 
without current supply to the wire) was recorded as 
function of the polar angle around the source. Alter- 
nately, a separately supported thermocouple probe 
was introduced into the cavity and temperature dis- 
tributions in Cartesian axes (¥ and X) were obtained 
for selected positions of the heater wire. 


(a) Characteristics of the Method 
The diffusion of heat from a two-dimensional source 
in uniform turbulent flow results in a thermal wake 

described by”” 
(T — T) = x—"? exp (— Y?/2Y?) - 


y? = 2xe Up 


LAMPBLACK STREAKLINES 
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separation 
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Fic. 2. Lampblack streaklines on the notch floor and oil-film visualization pictures for various configurations. (The light diagonal 


streaks in the separate wake picture are accidents of handling. ) 
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Spark-Schlieren (5 ywsec) photographs of wave radiation 
from cavities. 


Fic. 3. 


The coefficient of proportionality is evaluated from a 
heat balance, and ¢€ is the turbulent ‘‘eddy diffusivity. ** 

Fig. 5 shows (polar) traverses in a uniform flow 
at two Mach numbers at a constant radius of 0.060 
in. from the wire (0.010 in. diameter) for a fixed power 
dissipation. This constitutes a “‘calibration’”’ of the 
system. 

When the wire is placed in a steady but nonuniform 
flow the profile of the wake is different. With identical 
assumptions concerning the phenomenon, but with the 
mean velocity represented by 


Up — U,’ + y(Y) (3) 
one obtains for the profile of the thermal wake** 
(T — To) ~ x78 exp [—(Uo’ + vY)*/qervx] (4) 


The wake is skewed and the temperature maximum is 
found to occur at a polar angle from the source 


a = tan—! U)'/yx (5) 


measured from the vector flow direction. Note that 
it is conceptually possible to reconstruct the velocity 
profile from sequences of such wake surveys by integra- 
tion. 

In vortical flow one must also consider the possible 
effect of recirculation. The argument can be set forth 
most clearly by discussing a specific example. At 
the top of Fig. 6 are shown polar traverses around 
both a heated and an unheated wire placed at the 
centerline of a long (L/H = 10) notch close to the 
floor (0.015 in. above it). The unheated traverse is 
symmetrical; in this case, the thermocouple records 
the adiabatic recovery temperature on (approximately) 
a cylinder (7). The temperature increased by about 
10 degrees between the floor and the top of the traverse, 
indicating a higher velocity there, as expected. Also 
shown is (7. — 7), the local difference between the 


CAV try 
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flow. 
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readings with the source heated and unheated. A 
strong temperature lobe in the direction opposite to 
the external flow is found, indicating an upstream flow 
near the wall in the cavity; however, a relatively large 
lobe in the direction of the external flow is also found 
The temperature must decay in the direction of local 
If the flow is steady and the lobe on the down- 
stream side of the wire is due to heated air being re- 
circulated around the cavity, then the temperature 
must decrease upstream of the wire, but increase down 
This is not the case, as shown by the 
lower portion of Fig. 6 that the 
back lobe is not due to recirculation but to a reversal 


stream of it. 
One concludes 
of the flow during part of the time 

This striking confirmed by 
proves that the entire cavity is profoundly unsteady, 
not only in the sense of random “‘turbulent’’ fluctua- 


result, further tests, 


tions, but in the sense of definite reversals of the direc- 
tion of the flow. 

Consider the effect of the 
When it occurs, the magnitude of the temperature dif- 


flow reversal further 
ference measured at a given distance from a (constant 
power input) source can be, under proper conditions, 
larger in the direction of the smaller velocity. Tf one 
imagines a pulse with velocity |’, followed by one in 
the opposite direction with velocity V, (for equal 
time), one finds from Eq. 2 that the magnitude of the 
average temperatures in the plane of the flow at equal 
distances from the source in one direction (subscript 


F) and in the other (subscript 8) is: 
(Tr — 1o)/(Ta — 10) ~ V Vp Vy (6) 


Consequently, the direction of the larger diffusion 
lobe does not necessarily indicate the direction of mean 
mass flux. To ascertain this one must compare the 
results with measurements of dynamic pressure; the 
larger pressure reading will indicate the direction of the 


time-averaged flow. 
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Fic. 5. Polar temperature traverses around a heated wire in free 
tunnel (subsonic flow) at two Mach numbers. 


(b) Heat Diffusion Tests at the Centerline of the Notch 

Fig. 7 shows a series of heat-diffusion polar surveys 
at different heights above the notch floor. They are 
presented in two groups for clarity. These results 
must be interpreted in conjunction with stagnation 
pressure surveys (see Fig. 18 of reference 1) keeping in 
mind the characteristics of the heated wake in non- 
steady flow. 

There is a layer of reverse flow near the floor (Fig. 
7a) evidenced by the pronounced temperature lobe 
pointed upstream and the stagnation pressure readings 
which show flow in that direction. This flow is un- 
steady, and reverses part of the time, as evidenced by 
the existence of a smaller downstream-pointing lobe. 
The magnitude of the main lobe increases to a distance 
of 0.066 in. above the floor, indicating a decreasing 
velocity. At about 0.100 in., the forward and back 
lobes are almost equal in size*; this indicates that there 
is essentially no preferred direction of the flow. The 
height at which this happens coincides approximately 
with the height at which no definite dynamic pressure 
was recorded. 

Above this layer one finds another layer (0.100 to 
0.150 in. above the floor) in which stagnation pressure 
measurements show a downstream mean flow but the 
temperature surveys show an upstream-pointed lobe 
(Fig. 7b). This is not contradictory if one considers 
Eq. (6). The flow in this ‘‘buffer’”’ layer reverses direc- 
tion periodically. Physically, this suggests pulsations 
of the shear layer in and out of the cavity; when the 
shear layer dips in, the flow in the upper half of the 
cavity is entrained downstream; when it bends out, 
the flow in the cavity reverses. It confirms the spark 
photographs discussed previously. 

Above the buffer layer, both pressure and _heat- 
diffusion traverses indicate a strong downstream flow; 
this is the free shear layer proper. 

Similar studies of heat diffusion were also made 
near the recompression and the separation step. How- 
ever, the results were not interpretable. In both 
corners of the notch an almost uniform temperature 


* The asymmetry in the temperature traverses is due to the 
gradients over the finite distance between the source and the 
thermocouple. 
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was measured around the wire. These corners con- 
tain strong vortexes: recirculation and instability both 
contributed to the diffusion of heat. Consistent 
pressure surveys with which to monitor the interpreta- 
tion could not be obtained. 

Traverses were also made at several streamwise 
stations in the notch for fixed locations of the heated 
wire (at a fixed distance from the floor). From Fig. 8 
we see that: 

(a) The heat diffused from a wire in the recom- 
pression zone remains contained in a well-defined re- 
gion occupied by the recompression vortex. It escapes 
into the free flow past the edge of the notch. 

(b) The diffusion of heat from the wire at the center- 
line proceeds in both directions about equally, indi- 
cating a highly disturbed, direction-reversing flow. 

(c) The heat from the wire in the separation corner 
remains in the low velocity separation vortex and 
diffuses from there downstream through the cavity. 


(c) Summary of Observations 

The diffusion studies did not yield quantitative 
information. They did, however, prove that the flow 
in the cavity is definitely unsteady. The inability to 
obtain a mass balance in the cavity from pressure 
measurements alone (see Fig. 18 of reference 1) is 
clarified. The cavity appears to contain a strong 
vortex at the recompression corner, and also a weaker 
one of opposite sign at the separation corner. The 
flow in the central section exhibits three vertical 
layers: a reverse flow layer near the floor, unsteady, 
relatively weak, but definitely in the upstream direc- 
tion on the average; a buffer layer characterized by 
strong flow reversals with essentially no definite mass 
flux; and finally the shear layer. 
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Fic. 7. Polar heat diffusion traverses at notch centerline at various distances from the floor (V), Jo = 2.8,and L/H = 9: (a) return flow 
layer, (b) buffer layer. 


Heat Transfer Measurements 


a) Experimental Method 

Practical considerations and analysis suggest a com- 
parison of the transfer characteristics of separated re- 
gions with those of an “‘equivalent’’ attached boundary 
layer. In the present investigation the heat-transfer 
properties of the separated region are always referred 
to those of the attached boundary layer two to five 
notch depths ahead of separation, which were measured 
in each test. This station was chosen because the 
heat transfer is affected slightly by the expansion 
around the corner of the notch. The virtual distance 
from the origin of the oncoming boundary layer, S, 
as determined from traverses and Eq. 1, was long 
(about 30 in.) compared to the length of the notch 
(maximum 4 in.). Consequently, the Reynolds num- 
ber and the heat transfer characteristic of the attached 
boundary layer (which depend on the one-fifth power 
of S) were almost constant over the length of the 
notch and could be unambiguously defined by their 
value at the reference station ahead of separation. 
There was a slight change in the boundary-layer thick- 
ness associated with each Mach number. Table | gives 
the pertinent conditions. In the subsequent discus- 
sion, the tests are referred to by the Mach number 
and the boundary-layer classification “thin” or “‘thick.”’ 

The use of the heat-transfer coefficient for the flat 
plate (i,) just ahead of the notch for reference is ob- 


viously not universal because the history of the on- 
coming boundary layer is variable. In order to ascer- 
tain its possible influence, a comparison of the results 
for different upstream wall-temperature histories is 
shown in Fig. 9. The upper portion of the figure 
defines the conditions (obtained by changing the 
coolant distribution in the wall) while the lower por- 
tion shows corresponding ratios of local notch-floor 
heat-transfer coefficient to that of the standard station 
ahead of separation for a particular notch geometry and 
Mach number. This example shows a difference of 
approximately 10 percent in the ratio. It indicates 
that in the case of an adiabatic approaching boundary 
layer, such as in Larson's investigations,'® a still lower 
ratio is to be expected and that, in general, the previous 
thermal history of the approaching flow does bear some 
influence on the problem. 

It was questioned whether there is a difference be- 
tween flow with heat transfer to the wall and flow with 
heat transfer fo the fluid (such as was investigated in all 


TABLE 1 








Conditions for Heat Transfer Tests 
ons 6/H Reynolds no. 
Mach no. at Thin BL Thick BL at separation 
separation (Mo) (H = 0.50in.) (H = 0.25 in.) x 10-* 
2.1 0.85 Pres 4.8 
2.9 1.08 2.16 3.6 
3.5 1.14 2.27 2.3 
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other studies!'.!2). Two tests were conducted, in one 
of which the air was heated and the wall cooled (stand- 
ard configuration) while in the other the air was at room 
stagnation temperature and the wall was heated (with 
saturated steam, 300°F). The heat transfer coefficient 
ratios of both tests were found to lie within the data 
scatter. However, the accuracy of the measurements 
was relatively low with so low an overall temperature 
potential (lower than in any standard tests by a factor 
of approximately 3 to 4). 

All tests discussed subsequently were conducted 
under upstream conditions marked ‘“‘standard”’ in 
Fig. 9a. Fig. 10a shows a summary of all the measure- 
ments of the flat-plate heat-transfer coefficient, defined 
as follows: 

dr . Lo F,) . 
hy = — ; = KkK-— = (4) 
(Few — Fy) Fa Fe) 
K is a proportionality constant involving the module 
geometry and the conductivity of the wall material. 
Tw is the adiabatic wall recovery temperature com- 
puted from the relation: 


fem ; a I 9 ; 

Tr = pr +t (1 — ge) | 1 + m M,? (S) 
in which the recovery ractor nr was assumed equal to 
the one-third power of the Prandtl number: 


ne = (Pr)"3 (9) 


Compared, for instance, with the Reshotko—Tucker™ 
correlation, these data exhibit a somewhat larger de- 
crease with M/y; at My = 2 h,;is about 15 percent above 
and at MJ, = 3.5 it is about 30 percent below the 
values given by Eq. (34) of reference 25. The spread 
of the data includes all notch geometries, but there is 
no consistent effect of the notch L/H. 

Rewriting Eq. 7, it is possible to check the unknown 
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Fic. 8. Heat diffusion traverses at constant X/L, My = 
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Fic. 9. The effect of the temperature history of the on- 
coming boundary layer on the heat transfer to the cavity at M 


= 29and L/H = 6. 


adiabatic wall recovery temperature* from tests at 


different stagnation temperatures: 


q h te ns 5 h (: iar) 
T, 1+ (h/k) 7 1+ (h/k) To 


(10) 


All the data for a particular Mach number are shown 
in Fig. 10b plotted in this form. The correlation line 
shown is constructed using Eq. (9) to determine the 
intercept, and the mean value of the measured transfer 
coefficients (Fig. 10a) is used to determine the slope; 
it is seen that Eq. (9) is valid within the experimental 
accuracy. 

However, Eq. (9) applies only to an attached bound- 
ary layer. It has been found in other investigations*.‘ 
that in the separated regions ne differs somewhat (by 
2 to 5 percent on the average), although the available 
data are not conclusive. In view of this uncertainty, 
we used a simpler heat-transfer coefficient h defined 
on the basis of the adiabatic recovery temperature 
before separation, Ty ;. One has: 

h = hye. Taw _ le (11) 

awf ~~ lw 
It is evident that for small heat-transfer rates (7,, > 
Taw) the ratio of this to the true local h in the cavity is 


* The test set-up could not be insulated adequately so as to 
permit 7, to be determined experimentally. 
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very sensitive to small variations in 74y/ Tous. How- 
ever, under the present test conditions, one would not 
expect the heat-transfer coefficient / to be in error by 
more than 5 to 10 percent of Ajoc. in the cavity. 


b) Results 


Fig. 11 shows the distribution of the heat-transfer 
coefficient ratios over the floor of a (rectangular) notch 
at a fixed Mach number as a function of L/H for the 
thin and thick boundary-layer tests. 

Individually, the curves show a drop below unity at 
the separation point, a minimum in the first half of 
the notch, and then a rise well above unity at the re- 
compression step. The heat-transfer coefficient ratio 
increases with L/H near recompression, but decreases 
near separation. The ‘‘minimum”’ decreases with L/H 
and also shifts to a lower value of x L. At L/H = 9 
a “‘plateau’’ appears in the center of the notch, which 
is associated with the beginning of reattachment of the 
flow to the floor, and the ‘‘closure’”’ of the cavity (the 
critical closure length ratio is about 11 for the small 
6/H and about 12 for the large 6/H, see reference 1). 

The effect of increasing 6/H is to shift the curve to 
higher h/h,, while the curves remain qualitatively 
similar. 

Fig. 12 shows a parametric representation of the 
behavior of the heat-transfer coefficient ratio in the 
domain of My and L/H. The curves are for constant 
x L except for the curves marked “minimum,” which 
were drawn from values of the lowest the ratio h/h, in 
the notch. This value is about 5 percent lower than 
that at the separation corner.* It is seen that the 
general level of the heat-transfer ratio is affected 
slightly by Mach number, less in the thin boundary- 
layer than in the thick boundary-layer tests. This is 
particularly noticeable at the station at 70 percent of 
the notch length. When the boundary layer is thin, 
this station is still essentially in the low-heat-transfer 
region of the notch, but when the boundary layer is 
thick, it is already affected by the rise towards the re- 
compression point. This reflects strongly on the 
average heat transfer over the entire floor. 

Fig. 13 shows data for a downstream-facing separa- 
tion step (infinite 1/H). The minimum value of the 
heat-transfer ratio occurs immediately behind separa- 
tion. It is followed by a rise to a maximum of about 
1.1 at a value of x/H at which one expects reattach- 
ment of the separated shear layer (see reference 1), and 
then decays to a value somewhat lower than unity as 
the equilibrium boundary layer is_ re-established. 
The final value below unity is to be expected, because 
the reattached boundary layer is thickened. There is 
a slight Mach number effect which may, however, be 
due to the 6/H effect associated with changes in J) 
(Table 1). 

* The slight rise of h/h; at separation is thought to be caused 
by the heat transfer to the cavity from the (noninsulated) inner 
vertical face of the separation step. The separation module is 
exposed to the attached boundary layer on its top surface and 
its temperature is higher than that of the first floor module. 
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The minimum h/h, of Fig. 13 is consistent with the 
trend exhibited by the minimum lines of Fig. 12b with 
increasing L/H. The downstream-facing step under 
the thin-boundary-layer conditions could not be tested, 
but it is to be expected that the same general trend 
would be found with a minimum value at separation as 
low as h/h; ~ 0.2 (see Fig. 16a). 

Fig. 14 shows an integrated average heat-transfer 
coefficient ratio for the floor of the notch. It is seen 
that when the boundary layer is thin, the separated 
flow definitely “‘shields” the wall by a factor of about 
2. This is drastically different in the case of a thick 
boundary layer. The changeover seems to occur for 
a value of the ratio 6/H of about one, the same limit 
that was found to differentiate between “thin” and 
“thick” boundary layer behavior in the pressure dis- 
tribution. A slight Mach number effect is evident; 
again, this may reflect the influence of the variation in 
6/H rather than in the free-stream Mach number 
proper. 

In the foregoing, only the floor of the notch was con- 
sidered. The overall heat transfer per unit projected 
area to a square (H X H) section of the recompression 
step is presented in Fig. 15 for varying L/H and M,. 
The step was made of high-conductivity material 
(brass) so that the surface temperature was nearly 
uniform,* and it was insulated from the backup plate 


* The measurement does not, however, depend on this as- 


sumption 
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Fic. 10. Calibration tests: (a) heat-transfer coefficient 
ahead of separation; (b) correlation with assumed adiabatic 
recovery factor. 
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as shown in the sketch. The local temperature dis- 
tribution and heat transfer over its surface could not 
be measured. We note that, while gr/gr based on 
projected area is higher than unity, when based 
on exposed area, it is less than unity for most of the 
tests. 

A further study is shown in Fig. 16 for a specific 
L/H and My, including two different geometries of 
the step. The upper curve shows the heat transfer 
along the notch per unit projected area for both ge- 
ometries, and the wedge-shaped step is seen to reduce 
this ratio. The lower curve shows the same data 
based on exposed surface area; on this basis the ratio is 
less for the square step. It seems that there exists a 
separation ‘bubble’? immediately downstream of re- 
compression, which is more pronounced for the rec- 
tangular corner and partly ‘‘shields’” the wall. The 
use of projected surface is practically more significant: 
in terms of it the heat flux to the recompression zone 
It is decreased for slanted recom- 


’ 


is relatively high. 
pression steps. However, the shape of this step does 
not seriously affect the heat transfer either upstream 
in the notch or downstream. 

Some measurements were also obtained in a double 
notch by removing the backup plate. The second 
notch is evidently somewhat affected by the first notch 
when the reattachment surface is short. The general 
level of heat transfer in the second notch is higher and 
the distribution exhibits a ‘‘bulge’’ near the center of 
the notch which does not exist when the oncoming flow 
is uniform. The difference between successive notches 
could be eliminated by lengthening the reattached zone 
between them. 


Subsonic Flow 


A limited number of tests were mace in uniform sub- 
sonic flow with the same model. In general, there are 
no basic differences in the heat transfer to the notch in 
subsonic and supersonic flow, except that the transition 
from the “open” to the “closed” cavity regime is 
smoother, and the magnitudes of the coefficient ratio 
L/H are lower on the average, in subsonic Mow. The 
same was found in the study of the pressure field. Fig. 
17 shows a parametric map similar to that of Fig. 12, 
for My) = 0.42 and a thin boundary layer. 


Discussion 


The results of this study of heat transfer to cutouts 
in the boundary are in general agreement with trends 
already established in previous investigations (in par- 
ticular, reference 12). Direct comparison of the data 
is difficult because of adjustments needed to compen- 
sate for the different experimental models and methods. 
However, when these differences are considered* one 

* For instance, reference 12 does not distinguish between the 
“floor” and the recompression step; the oncoming boundary layer 
is very thin; and caution must be exercised in interpreting the 
definition of the reference attached boundary layer. 
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finds good agreement in the magnitudes of the numerical 
results. 

The present work crystallized some of the secondary 
trends also, and threw light on the detailed structure 
of the flow. Between this work and the NACA pro- 
gram, the behavior of separated cavities in a restricted 
domain of Mach number, Reynolds number, cavity 
geometry, and initial conditions has become em- 
pirically fairly well known. 

However, experiment failed to substantiate the valid- 
ity of the shear layer model. Good agreement between 
the theory and experiment is found for the /aminar two- 
dimensional flow, but this does not constitute a satis- 
factory argument. A 
associates the entire heat transfer problem with the 
nature of the shear layer must be judged by how it 
responds to changes in the properties of the shear layer 
Chapman’s analysis can be said to provide a good corre- 
lation equation for /aminar separations, but not for 
transitional and turbulent heat transfer through the 
observed heat transfer 


theory which emphatically 


cavity. Furthermore, the 
gradients along the floor of a cavity are quite foreign 
to the model, and noting their magnitude one finds it 
difficult to neglect them. 

As already mentioned in the Introduction, there are 
three transfer mechanisms at play: (a) conduction 
between the wall and the ‘dead air,” (b) conduction 
through the shear layer, and (c) mass convection by 
exchange of fluid between the cavity and the external 
flow. 
dominate the others either throughout the range of 
variables or at all, but such a hypothesis is the basis 
for postulating a simple “‘model.’”’ The Chapman 
Larson assumption that (b) dominates does not check 
with experimental measurements, at least not in gen- 
eral. Carlson’s analysis, which stresses (a), is simply 
too arbitrary; besides, the present measurements of 
the internal flow definitely disprove his particular 
assumption of the return-flow velocity profile. They 
indicate, to the contrary, that the conductivity of the 
wall film is very high due to an extremely turbulent 
unsteady and vortical flow. The mechanism (c) is 
introduced here for the first time. The observations 
reported earlier show that it is not a negligible effect, 
and it is reasonable to study the possibility that it 
dominates the heat transfer problem. 


It is not necessary that one of these mechanisms 


The Mass-Exchange Model 


The cavity “breathes’’ at ultrasonic frequencies. 
During one half of the pulse mass is injected into it. 
Since the heat transfer to the wall inside the cavity is 
high, the excess fluid heats up to wall temperature and 
carries this heat away when ejected from the cavity 
during the ejection half of the pulse. The rate of heat 
exchange depends thus on the unsteady mass exchange. 
Schematically, the internal flow is as shown in Fig. 
18. The mechanism is this: the separating stream- 
line fluctuates to an extent + A about its mean position 
so that a mass dw is periodically fed into the cavity. 
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Fic. 11. Variation of the heat-transfer distribution on the floor of the cutout as function of L/H at Mo = 2.9: (a) thin oncoming bound- 


ary layer; (b) thick oncoming boundary layer. 


This fluid is recirculated around the recompression vor- 
tex, picked up in the shear layer, and discharged past 
the recompression point. The momentum impulse 
associated with the mass dw must be proportional to 
the centrifugal impulse of the mass dw whipped around 
the recompression vortex. This relation fixes the 
magnitude of the pulsation, assuming that the profile 
in the neighborhood of the mean streamline is known. 
A simple formulation of this mechanism is outlined 
in the Appendix. It leads to the following result for 
the mean heat transfer to the separated cavity: 
g 1 6s 
aa 


St = ee = 
Gu of em = (09) Lb 


(12) 


where 6, is the thickness of the shear layer just ahead 
of recompression and ¢ is a function of the shear layer 
form describing its slope at the mean separating 


d(U Uo) P 
— —_ eae 13 
kz 6.) i. ne 


Comparison With Experiment 


streamline: 


Provided that the oncoming boundary layer is thin 
relative to the notch, the effect of it on the development 
of the free shear layer is small. The function ¢ is a 
universal function for either laminar or turbulent flow. 


From Chapman’s laminar solution one finds: 
gd, = 1.06 (14) 


For turbulent flow one finds, using Goertler’s solution 
of free-jet mixing: 


U/Us = '/2{1 + erf (c Y/x)} (15) 
¢=0.8 (16) 


The numerical values depend on how one defines 6,, 
but it is obvious that the turbulent and the laminar 
values of ¢ are essentially equal. The entire effect of 
transition and Reynolds number is contained in 6,/L, 
the growth rate of the free shear layer: 

(a) Since upon transition to turbulence (at a fixed 
Re) 6,/L certainly increases, one expects an increase in 
heat transfer. However, it seems reasonable to pre- 
sume that the transitional thickening of a free shear 
layer would be of the same order of magnitude as the 
thickening of an attached boundary layer; therefore, 
the ratio h/h, would not be expected to change very much 
upon transition. These results are borne out by the 
data.'” 

(b) In laminar flow 6,/L ~ Re'/?, therefore 


Stiam. ~ Re~? and (h/hy)iam. = const. (17) 


as verified by Larson’s data. 
(c) The rate of growth of a free shear layer in turbu- 
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lent flow is not well known. However, referring to 
the measurements of the thickness of the recompression 
profile presented in Part I of this series,! one finds 
that the ‘‘thin’’ boundary-layer data follow the trend 
6,/L ~ Re-?!> very closely, so that Larson’s result, 


9/ 


Stu. ~ Re~*> and (h/hy)turn. ~ Re~™® (18) 


is consistent with this model as well as with the inde- 
pendent pressure measurements. 

(d) Eq. 12 predicts no Mach number dependance 
in either laminar or turbulent flow, as long as the 
Mach number is low enough so that 6, is not altered by 
compressibility effects. This is essentially in accord 
with the observations of Larson and with Fig. 14. 

(e) Few theoretical arguments can be used in the 
case where the oncoming boundary layer is_ thick 
enough to influence the development of the shear layer. 
However, it seems certain that a thick boundary layer 
at separation will result in a thicker shear layer at re- 
compression, that is, 6, will increase with 6. Conse- 
quently, the heat transfer will increase. Fig. 14 bears 
out this conclusion. 

(f) The trends shown by Fig. 9 relating to changes 
in thermal history of the oncoming flow are in qualita- 
tive agreement with the model. When the oncoming 
boundary layer is adiabatic, the fluid which penetrates 
the cavity is at a higher initial temperature, absorbs 
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less heat from the cavity wall, and the heat transfer 
rate is decreased. 

(g) Finally, the proposed structure of the flow inside 
the cavity is compatible with the observed distribution 
of heat transfer over the floor of the notch. One would 
expect the highest rate of heat transfer from the 
boundary to occur near the recompression zone, the 
remainder of the cavity contributing little to the 


exchange. 


Limitation of the Model 


The mass-exchange heat transfer model exhibits all 
the trends indicated by experimental data. Should 
one wish to extend the argument, the numerical values 
obtained from Eq. 12 (constant of proportionality in 
Eqs. A5 and A6 equal one, Chapman’s laminar solution 
for 6,/L, and Eq. 14) also check the measured Stanton 
However, these confirmations must be 
The present model depends on 


numbers. 
taken with reserve. 
the unsteadiness of the flow at the exclusion of the other 
mechanisms. Evidence of the existence of pulsation 
has been presented here (Fig. 6) and elsewhere,’ but 
there is also evidence that not all cavities pulsate 
For instance, preliminary work with cavities having 
slanted (wedge-shaped) recompression faces indicated® 


that no wave radiation occurred. In such cases, the 
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boundary layer; (b) thick boundary layer. 


(Note: points are averages over all tests.) 
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Fic. 13. Heat-transfer coefficient ratio behind a downstream- 
facing step. Thick boundary layer 


analytical model under discussion is clearly question- 
able. 

One suspects that neither the Chapman—Larson 
‘“conduction’’ model nor the present mass-exchange 
model are sufficiently complete. Perhaps both mech- 
anisms share in the process in varying proportions, 
depending on the geometry, and since their contribu- 
tion to the transfer of heat is of the same order, it is 
difficult to distinguish between them. However, the 
present analysis seems to yield a more accurate semi- 
empirical correlation of the measurements over the 


entire range investigated to date than other models. 


Concluding Remarks 


The ratio of the heat transfer to an “open” cavity to 
that of an equivalent attached boundary layer is pri- 
marily a function of the thickness of the oncoming 
boundary layer. If it is thin (that is, less than one 
notch depth) the heat transfer ratio is reduced by a 
factor of approximately two; if it is thick, the heat 
transfer can be even increased. 

However, the heat transfer per unit surface area is not 
necessarily the important criterion. The drag of a 
notched wall is an order of magnitude higher than that 
of an equivalent flat plate. Therefore, the heat transfer 
to a notched wall per unit drag is around ! » of that to 
a plate. It is possible, subject to an extension of these 
findings to hypersonic Mach numbers, that this fact 
can be useful in re-entry problems. 

It is proper to add that preliminary investigations 
of notches in the boundary of both supersonic and sub- 
sonic accelerating flows showed that the external 
pressure gradient has a strong influence on pressure 
and heat-transfer distributions. Some pressure meas- 
urements are reported in reference 26. Other meas- 
urements showed a twofold increase in heat transfer to 
the wall under these conditions. This is currently 


under investigation. 


Appendix 


If the shear layer dips in and out by a distance 


y + A (see Fig. 18) about its mean position (y = 0, 
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u = Ux) as suggested by the results reported in the 
text, then there is an excess of momentum per pulse 
fed into the cavity, which is 


TA 
dM = p| u*dy — 2pux*A = ?/3pf72A? (A1) 

—A 
where the velocity in the neighborhood of the origin 
is assumed to vary linearly with distance: 


u= us +Ty (A2) 


This is associated with a periodic mass pulse 
A 
dw+ = —dw- = of udy — pudA = pl'A? (A3) 
0 


The impulse associated with the momentum transfer 
to the cavity and the centrifugal force due to the un- 
balance mass whipped around the recompression vor- 
tex must be in balance. The measurements of acoustic 
radiation suggest that the rim speed of the cavity vortex 
is proportional to the external flow velocity while its 
radius is of the order of the notch depth. Then on 
dimensional grounds”*: 

Momentum impulse dM 

= const. - (A4) 


(dw Uo?/H)r 


centrifugal impulse 
where 7 is the duration of the pulse governed by the 

circular frequency of the vortical motion, 
t~ H/U) (A5) 


* This dimensional argument does not take into account the 


effects of compressibility. 
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Fic. 14. Integrated average heat-transfer coefficient ratio for 
the floor of the cavity in the Mach number and L/H domains for 
thin and thick boundary layers 
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Fic. 15. Heat flux per unit projected area to the recompression 
step. Thick boundary layer. 


so that one can relate the gradient of the shear-layer 
profile to the shear-layer thickness: 


i A _ Uo d(U Up) Uy 
Y=0 


= = (A6) 
dY ? 


5, d(Y/6,) 6, 


where 6, is the thickness of the shear layer while ¢ is a 
shape factor discussed in the text. 

Introducing these into Eq. (A4), one can solve for 
the‘unknown which is 


A ~ 6,/¢ (A7) 


The unsteady mass exchange per pulse (note that the 
time-averaged mass flux per unit span to the cavity is, 


of course, zero) is: 
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Fic. 16. Heat-transfer coefficient ratio behind recompression for 
a thick boundary layer and various geometries. 


The heat absorbed and carried away by the fluid 
flowing through the cavity per unit span per second is 


g = CAW( Lan — Tp) (A9 


where 7, ‘s the temperature of the incoming fluid 
For a cavity of the length Z one has, using Eq. (AS): 


q l 6, 


— ————_ = (A10) 
ee . T)CppUo r79) E 


St = 
or, relating the transfer properties of the cavity to those 
of an attached boundary layer just ahead of separation 


for which 


St, ~ Re 


Lo 
h/hy 


os 
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Fic. 17. Map of heat transfer to a cavity in subsonic flow: 
M = 0.42, 6/H = 0.28, turbulent flow. 


where 7 is 1/2 and 1/5 for laminar and turbulent flow, 


respectively, one has: 


St h l (¢ 
a, b+ 6 
It should be remembered that 6, is the thickness of the 
shear layer and not the thickness of the approaching 


boundary layer, while Ke relates to the characteristic 
of the boundary layer at the conditions selected as 


Re” (Al l ) 
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Fic. 18. Structure of the internal flow for the mass-exchange 
heat transfer model. 
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On Aircraft Maneuvers in the 


Vertical 


Plane Under Assumptions Leading to Easy 
Integrability of the Fliight-Path Equations 


F. EDWARD EHLERS* ann GEORGES BRIGHAM*#* 
Boeing Airplane Company 


Summary 


The differential equations for the flight path of an airplane in 
the vertical plane with constant acceleration of gravity are 
integrated in closed form under the following sets of assump- 
tions: 

(1) Net thrust-to-weight ratio and lift-to-weight ratio are 

constant. 

(2) Centrifugal acceleration and net thrust-to-weight ratio 

are constant. 

(3) Velocity and lift-to-weight ratio are constant. 

(4) Angular velocity and net thrust-to-weight ratio are con- 

stant. 
Graphs and tables of some of the solutions are provided. 


Symbols 
a = a/K 
D = drag 
g = acceleration of gravity 
h = B — cos# 
ho = B-1 
K = number of g’s of centrifugal acceleration 
L = lift 
t = time 
r = thrust 


velocity of airplane 
initial velocity for 6 = 0 
= initial velocity for upside-down position, 6 = 7 


* 
ll 


WwW = weight of airplane 

(x, y) = Cartesian coordinates of airplane position 

z = 1/K 

a = (T — D)/W, net thrust factor 

B = L/W, load factor 

o = g/vw 

o = angle in radians defined by Eq. (9) 

pi = angle in radians defined by Eq. (22) 

¥ = (a-0)V2/2 

Vo = @ 2/2 

& = gx/m?) 

n = gy/v?> dimensionless variables 

T = gt/u \ 

60 = angle of flight path measured with respect to x-axis 
as horizontal position 

w = angular velocity of airplane 


Introduction 


N THE PRELIMINARY DESIGN OF AIRCRAFT, it is con- 
venient to have simple means of estimating the 
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or 
1° 4) 


flight path, angle, and speed of the aircraft undergoing 
some maneuvers in the vertical plane. For example, 
such problems may arise in the study of bomb-dropping 
techniques and of evasive tactics of fighter craft. Also, 
one may be concerned with the loss of altitude which 
an airplane must undergo in order to right itself from 
an upside-down position. To facilitate such investi- 
gations, some closed-form integrations of the equations 
of the flight path of an airplane have been obtained 
for the following sets of assumptions: 


(1) Net thrust-to-weight ratio and load factor (lift- 
to-weight ratio) are constant. 

(2) Centrifugal acceleration and net thrust-to-weight 
ratio are constant. 

(3) Velocity and load factor are constant. 

(4) Angular velocity and net thrust-to-weight ratio 
are constant. 


These solutions may also be used in estimating more 
complicated flight paths by considering the load factor or 
net thrust factor, for example, as varying in a piecewise- 
constant manner, to approximate the desired con- 
tinuous relation. Graphs and tables of some of the 
solutions are provided in dimensionless variables to 
facilitate use by the engineer. 


Equations of Motion 


Let @ be the angle between the flight path and the 
horizontal direction taken as the x-axis. Defining the 
lift L as the force on the airplane normal to the flight 
path and the thrust 7 and drag D as always tangent to 
the flight path yields the following equations of motion 
for the airplane: 


dv/dt = g |((T — D)/W — sin 6] (1) 
v d0/dt = g(L/W — cos 8) (2) 
dx/dt = vcos6, dy/dt = v sin @ (3) 


where v is the speed of the airplane, W the weight, g 
the acceleration of gravity, and x and y are the 
Cartesian coordinates of the airplane position. 


(I) Maneuvers With Constant Load Factor 
and Constant Net Thrust Factor 
The Velocity as a Function of the Angle of Flight 


If we assume that the airplane is flown in such a way 
that the load factor L/W and the net thrust factor 
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(I. — D)/W are maintained constant, then Eqs. (1) 
to (3) can be integrated in closed form. Let 


(T —D)/W=a 
be the net thrust factor, and 
L/W = 8 


the load factor. The equations of motion then become 


dv/dt = g(a — sin 6) (4) 
and v d6/d! = g(8 — cos @) (5) 


Solving Eq. (5) for dt yields 
dt = v d0/g(8 — cos 6) (6) 


Eq. (4) now can be integrated if we replace dt in terms 
of @and v by Eq. (6). This gives 


dv/v = [(a — sin 0)/(8 — cos 6) |dé 


dv/v = [add/(8 — cos 6)| — 
[(sin 6 d@)/(8 — cos @)| (7) 


Integrating from 0 to 6 yields 


v/v = [(8 — 1)/(B — cos 6)|e*” (8) 
9 ‘ 
- B+1 g 

where ¢ = Vv tan—! | ; tan | (9) 
B? aad l B _ l = 


and v = vy até = 0. 


The Coordinates as Functions of the Flight-Path Angle 

The coordinates x and y are found as parametric 
equations in @ when dt in Eq. (8) is replaced with dé by 
Eq. (6). Thus we have 


6 
o * f [v? cos 6/g (8 — cos @)] dé 
0 


6 
y= f [v? sin 6/g (8 — css 6)] dé 
0 


Introducing the dimensionless variables 


Vo 


/ 9 f 
— = gx/u", n = gy/u*, T= gt, 


2 


LOSS IN ALTITUDE 92 


w 
Fic. 1. Loss in altitude for an airplane righting itself from an 
upside-down position. Constant load factor 8 and constant net 
thrust factor a. 


and eliminating v by Eq. (8), we obtain 


6 
t= it f [e?** cos 6/h*] do (10) 
0 
° € 
n = hi? f e~“* sin 6 d6/h? (11) 
0 
6 
r= inf e**d6/h? (12) 
0 
where h = 8 — cos @ and fy = B — 1. 


The first step in evaluating the integrals in (10), 
(11), and (12) is to express them in terms of the func- 
tions of the form 


6 
F(a, 0) = f e** dé/h" 
0 


xg 
“c 


Fic. 2. Comparison of loop flown at constant load factor with 


= = )- 


loop flown at constant centrifugal acceleration. 6 = K = 
a = 0.50. 


We then find a recurrence formula for the F,(a, @) by 
repeated integration by parts, using the relation d@ = 
d6/h. Since (cos 0)/h = (8/h) — 1, Eq. (10) can be 
written 


& = h?(BF;(2a, 0) — F.(2a, 6)] (13) 


6 
+ 2 . 
Now f e~** sin 6 d6/h? = 
0 


1] » . © sin 0 
—| e** sin 6/h? — f ed 
2a 0 h? 


Differentiation of (sin @)/h? and use of the relation 
(cos 6)/h = (8/h) — 1 yields, for Eq. (11), 


n = (ho?/2a) (e?a* sin 0/h? — F,(2a, 0) + 
38F2(2a, 0) — 2(8? — 1)F;(2a, @)] (14) 


and for Eq. (12) t = hoF:(a, 0) (15) 
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Fic. 3. Flight paths flown at constant load factor 8 and constant 
net thrust factor a. B = 2,3; a = 0, 0.25. 0.50, 0.75. 


Consider first the function 


, 0 
Fi(a, 6) = { e* d/h = { e~? dob 
0 0 


This integrates easily to 
Fi(a, 0) = (e** — 1)/a@ (16) 


The second integral becomes 


6 
fe d6/h* 
0 
1 6 
= |<**/h — 1/ho +f ee sin 9do/h | 
a 0 
. ’ i... 
= e*/h — 1/ho}) + —| e** sin 0/h — 
a an 
6 " 
f ot] (= ‘) 
0 h 


Differentiation of (sin @)/h and use of the relation 


II 


F(a, #) 


(cos 0)/h = (B/h) — 1 
yields 


F(a, 0) = (1/a) [(e**/h) —(1/hy)| + (1/a?) le? sin 0/h + 
BFy(a, 0) — (8? — 1) F(a, 6) | 


Solving for /2(a, 6) we finally obtain 


a=3/4 
1/2 
roe ~ 
a=l/4 
=0 
n = \ 
B=4 ) 
err 
Fg ai 
Se 
=a 


Fic. 4. Flight paths flown at constant load factor 8 and constant 
net thrust factora. B = 4; a = 0, 0.25, 0.50, 0.75. 


F(a, 0) = { (e*® ‘h)(a + sin #) + BFi(a, 6) — 
[a/(8 — 1)|) {/(e? +82 — 1) + (17) 


Similarly, 


F3(2a, 0) = }(e°**/h?)(a + sin 0) — F,(2a, 0) + 
3BF.(2a, 0) — [a/(B — 1)*| }/2(a2 + B? — 1) (18) 


The flight path for which the thrust just balances the 
drag has special interest because of the simplicity of the 
formulas. When a = 0, then 


F,(0, 6) = @ 
F,(0, @) [sin 0/(8? — 1)h] + [6¢/(6? — 1)] 
and F;(0, 0) = [(sin 0/h?) — @ + 36F,(0, @) ]/2(6? — 1) 


II 


Substituting in the relation for £ and 7 we obtain, after 
some simplification, 


£ = }[8(6 — 1)/A] + [(62 + 2)/(8 + 1)] {sin 8 
2(8 + 1)h + [38¢/2(6 + 1)7] (19) 


and 7 = [sin 0/h(8 + 1)] + [6¢/(8 + 1)] (20) 


The expression for 7 can be found more easily by 
setting a = Oin Eq. (11). Hence, 


6 
7 = nef sin 6d0/h® = (1 — ho?/h?)/2 (21) 
0 


Downward Loop Under Constant Load Factor 
and Net Thrust Factor 

A problem of practical interest is finding the loss in 
altitude and the increase in velocity for an upside- 
down airplane righting its position by means of a half 
loop. This information could be computed from the 
formulas of the previous sections, but for greater con- 
venience we shall derive new formulas. We integrate 
Eqs. (7), (10), (11), and (12), replacing the lower 
limit by 7, and fy by 8 + 1. The velocity is then 


given by 


v = v;[(8 + 1)/(8 — cos 6)Je™ 


where @ : tan~ ah. s= 22) 
where = é —~4/— x 22 
Vg? 1 +i 


and v, is the velocity of the airplane at the horizontal 
upside-down position. The angle ¢, will vary from 0 
to 7/V/g2 — 1 as 6 varies from 7 to 27. The relations 
for £, y, and 7 are those of Eqs. (13), (14), and (15) 
except that /» is replaced by 6 + 1 and the functions 
F(a, @) are given by 


6 
F,(a, 0) = f e™ do/h” 


with ¢, defined in Eq. (22). Performing the integra- 
tions yields the same forms for F;(a, 0), Fo(a, @), and 
F;(2a, 0) as Eqs. (16), (17), and (18) with ¢; replacing 
@ and a/(8 + 1) and a/(8 + 1)? replacing a/(8 — 1) 
and a/(8 — 1)’, respectively. 
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The altitude lost by an airplane in righting itself is 
found for @ = 27. Hence 


(8 + 1)? 





,= — [— F,(2a, 27) + 38F2(2a, 27) — 
2a 
2(8? — 1)F;(2a, 27)| 
where 
e2anr/V/ B= 1 — 1 
F\(2a, 27) = enietiek Aaa 
2a 





oe jae 2a 
F,(2a, 27) = — a a Se ae 
om ae @+pij* 
BF,(2a, 2m) | /4a? + p? — 1 


‘i= Vv | _ F 5 | = 
iL @- D? (8+ 1)? 


F,(2a, 2m) + 38F,(2a, 27)}/2(a?+ 8? — 1) 


F;(2a, 237) = 


For a = 0, the loss in altitude reduces to 


. ceil '/, (Gt) + '/» 


A graph of altitude loss vs. 8 is shown in Fig. 1 for a = 
0, 1/4, and !/>. The loss in altitude decreases with in- 
creasing 6 and increases with thrust, as we should ex- 


pect. 


(II) Maneuvers With Constant Centrifugal 
Acceleration and Constant Net Thrust 


The Velocity as a Function of Flight-Path Angle 
If the centrifugal acceleration is maintained constant 
during flight, then Eq. (2) becomes 


v d0/dt = Kg (23) 


where K is a dimensionless constant. Eliminating di 
from Eq. (4) by means of Eq. (23) yields 


v 6 
f dv/v = f (a — sin 6)d6/K 
M% 0 


Integrating gives 


v/vo = exp | ab/K + (cos @ — 1)/K} 
or v/vo = exp {a0 + 2(cos 6 — 1)} (24) 


where 1 = a/K and z = 1/K. 


The Time of Flight as a Function of Flight-Path Angle 
The time is found by integrating Eq. (23) with v 
given by Eq. (24). Thus, with r = gt/uvo, we get 


T= : | exp {a0 + 2(cos 6 — 1)} d0 
0 


Tl . . zcos@ : 
To evaluate this integral we first expand e°“*" in a 
~ ° ° — zcos 0: e ° 
Fourier series. Noting that e*“*" is an even function 


of 6, we have 


ao 
F4 6 
e® *° = Fa, cos nb 
n=0 
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Fic. 5. Flight-path angle @ vs. dimensionless time 7- for flight 
at constant load factor 8 and constant net thrust factora. 8B = 
2.3; a = 0, 0.25, 0.50, 0.75. 


where 


1 T 9 r 
>os 6 = zcos @ 
a = f "dd, a. = f e* ©°* “cos n6 dé, n> 0 
T/0 7T</0 


These integrals can be expressed in terms of Bessel 
functions of imaginary argument by use of Eq. (4), 
p. 1S1, reference 2. Thus 

a = I(z), ay, = 21,(2) 


The integral then becomes 


e6 
= 2 =| 102) | e”’ dé + 
0 


, , ee 
2 yi I,(2) f e” cos nO ao | 
n= 0 


Since the series is absolutely convergent for all finite 
A > 1, we can integrate term by term. Hence, 


1)/a + 


T= ze~*{Io(z)(e”” _ 


2 >> 1,(z) [e* (a cos n6 + nsin n0) — al(a? + n?)} 





T 

Fic. 6. Flight-path angle @ vs. dimensionless time r+ for flight 
at constant load factor 6 and constant net thrust factor a. 8B = 
4; a = 0, 0.25, 0.50, 0.75. 
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TABLE 1 











ra = (2/Vp  — 1) tan [V (8 + 1)/(8 — 1) tan (¢/2)] 




















g | B-2 B=3 B=4 

0° .0000000 |} . 00000000 . 00000000 
10° . 1736561 . 08704633 . 05807968 
20° . 3422650 | .1728030 1155818 
30° . 5017389 . 2561270 . 1719792 
40° . 6494999 . 3361326 . 2268358 
50° . 7845022 41224037 . 2798304 
60° . 9068986 . 4841696 . 3307644 
70° 1. 0176002 . 9518939 . 3795529 
80° 1.1178973 | .6155770 . 4262086 
90° 1.209196 . 6755109 . 4708197 
100° 1.292872 . 7320634 . 9139302 
110° 1.370185 . 7856374 . 9545199 
120° 1.442270 . 8366488 . 9939929 
130° 1.510062 . 8854611 . 6321304 
140° 1.574592 . 9326010 . 6692512 
150° 1. 636566 .9783155 . 7054765 
160° 1.696763 | 1.023010 . 7410589 
170° Ll, FOOT 1. 067032 . 7762137 
180° 1.813792 1. 19C721 .8111558 








The Altitude as a Function of Flight-Path Angle 


The altitude is found by evaluating 


3 
Il 


6 
[ (v/v) sin 6 (dr/d@) dé 
Jo 


0 5 
= 2 f exp {2a + 22(cos @ — 1); sin 0 dé 
0 


Integrating by parts, after noting that 


d(e7#°s °—)) /d@ i — 2ge72(cos 6—1) sin 6 


we obtiin 


ACE 





7 = (1/2) {1 — exp[2a@ + 2z2(cos @ — 1)}}+ 
t] 
af exp [2a0 + 2z2(cos @ — 1)]| dé 
0 
This integral has the same form as the integral for the 
time 7, with 2z and 2a replacing z and a respectively. 


Then, following the same procedure used in evaluating 
7 leads to 


n = (1/2)[1 — (v/m)?] + 7) (1/2)1 (22) ~~ e 





2 ¢ 2a0/¢ : So of 
>> T,, (22) e*"(2a cos nO + n sin @) _ (26) 
4a? + n? 


n=1 


The Range as a Function of Flight-Path Angle 


Now the horizontal distance or range is given by 


t= f (v/v) cos 0( dr/dé) dé 
0 


6 
or gf=2 fexp [2a6 + 22(cos @ — 1)] cos 6 dé 
0 


22 cos 6 


Using the Bessel expansion of e , we obtain 
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t= ze f et) Ia(22)e0s 6+ 
0 


> I,,(2z) [cos(n — 1)@ + cos(n + alt dé 


n=1 


or combining like cosine terms, we get 


6 
é= zen f eh (2s) + 
0 


de [n-1(22) + Inss(22)] cos no d§ 


n=1 


Finally, integrating term by term yields 


Ss 


t= ze~**4(1/20)Iy(25) —1)+ 


dD (n-1(22) + Tngs(22)] [e?*” (2a cos n6 + 


n=1 
n sin n0) — 2a] /(4a? + nh (27) 


When K > 1 (zg < 1), then the series for £, 7, and 7 
converge very rapidly and only a few terms are needed, 
since 


I,(%) ~ 2"/n 


The calculations of £, 7, and 7 are not as formidable as 
it might seem, since tables of /,(x) are given in refer- 
ence 3. The functions Jp(x) and J;(x) ate also found 
in reference 2, and higher-order Bessel functions may 
be computed from the recurrence relation 


In4i(x) = T,-1(x) = (2n/x)I,(x) 


The Load Factor for a Turn at Constant Centrifugal 
Acceleration 

The equations for the trajectory do not depend ex- 
plicitly on the load factor 8. Under constant centrif- 
ugal acceleration the load will vary as a function of 0. 
This is given by combining Eq. (23) with Eq. (5), which 
yields 


8—-cosA@=K 








Flight paths at constant centrifugal acceleration K and 
K = 1,2; a = 0,025. 0 50, 075. 


Fic. 7. 
constant net thrust factor a 
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or L/W =K + cos @ (28) 


The cos @ term gives the influence of gravity. The 
load factor is seen to vary between the limits of A — 1 
in the upside-down horizontal position and K + 1 in 
the upright horizontal position. 


The Long Climb Under Constant Centrifugal Acceleration 


If K < 1, then the load factor for small angles is 
slightly larger than unity and the parameters a and 2 
are large. The velocity is then approximated by 


, e (a9 — 20 )/2 


sin e we are concerned only with small angles. The 


time of flight is given by the integration of 


8 
sa?/2 —s(a—6)2/2 
r = ze™ f co dé 
0 


This integral can be expressed in terms of the error 
function 


and becomes 

1 V 92/2 &?[Ex() * E2(Y)] 
Y = (a — 0)V/2/2 and = avV/z/2 
and the upper sign is taken when y > 0 and the lower 


sign when y < 0. 
Similarly, 


(29) 


where 


Str 


8 
= f (v/v) cos 6 (dr/dé) dé 
0 


6 
0 
I F avz a 
/2 e* ec dt— 
(a—0) Vz 


8 9 
0 


bho] = 


II2 


l avz -72 oe 
5 _e (a — t/Vz)? dt 
“ J(a-avVz 
ga (a 
ae 
ag 
en 
—— =. ™~ 
st 
= — ? 
ies 
—__ ‘’ 
n =O \ 
ay \ }—/ 
a \ / tA 
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Fic. 8. Flight paths at constant centrifugal acceleration K and 
constant net thrust factor a. K = 3; a = 0, 0.25, 0.50, 0.75. 
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a , Y 
K=2 y, “Kei 
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Fic. 9. Flight-path angle @ vs. dimensionless time 7 for flight 
paths at constant centrifugal acceleration A and constant net 
thrust factor a. K = 1,2; a = 0, 0.25, 0.50, 0.75 
Integrating yields 
f> Vz [1 — (a?/2) — (1/42)] X 
[E(V2yo) F E(o/2y)] + 
(1/4)(a + de?" —™ — (a/4) (30) 
Also, 
6 
n= f (v/vo) 6 (dr/d0) dé 
0 
Za? = —t? 
~ze e (a — t/~/z) dt 
(a-OA)Vz 
> a2 e™ [E(V/2 he) F E(V/2y) + 
(1 — e7*°~")/2 (31) 


These formulas are more amenable to calculation than 
the series of Bessel functions since they contain fewer 
terms and good tables of the error function are readily 
available (see reference 4, for example). 


(III) Flight Path With Constant Velocity 
and Constant Load Factor 


We shall now consider the maneuver with the ve- 
locity maintained constant. This is a realistic flight 
for a high-powered airplane over a limited portion of an 
upward loop, and for an airplane with speed brakes over 
a limited portion of downward loop, or for a transport 
plane which can use spoilers and extend its landing gear. 

The assumption of constant velocity leads to a rela- 
tion for the difference between thrust and drag in terms 
of the flight-path angle. From Eqs. (1) or (4) we have 


(T — D)/W = sin @ 
Eq. (5) becomes 
g dt/vy = d6/(8 — cos 8) 


If 8 is assumed to be constant also, then for initial level 
flight this integrates to 


(32) 
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Flight paths for constant load factor 8 and constant net thrust factor a 





















































(8 = 2; a = O, 0.25, 0.50, 0.75) 
B=2, a=0 B=2, a=025 

8 € 7 r WG € 7 r WM 

0° 0. 000 0. 000 0. 000 1.000 .00000 | .00000 | .00000 | 1.0000 
10° 0.189 0.015 0.172 .992 . 17865 . 01574 .17658 | 1.0288 
20° 0. 322 0. 055 0.335 .972 .35141 .06190 | .35033 1.0274 
30° 0.444 0.111 0.462 . 943 . 50044 . 13107 .51230 .9996 
40° 0.532 0.171 0. 609 .910 .61630 | .21175 . 65681 . 9532 
50° 0.589 0. 228 0.712 . 879 . 69831 . 29325 . 78180 . 8964 
60° 0.623 0.278 0.801 . 850 . 75129 . 36835 . 88791 . 8363 
70° 0. 643 0. 322 0. 869 . 823 . 78199 . 43357 . 97728 .7779 
80° 0.649 0.351 0.926 . 806 . 79680 .48811 | 1.05260 .7241 
90° 0. 653 0.375 0.973 .791 . 80085 .53264 | 1.11652 . 6765 
100° 0.652 0. 394 1.015 .778 . 79782 .56846 | 1.17133 . 6356 
110° 0. 646 0.409 1. 047 . 769 . 79026 .59699 | 1.21918 6014 
120° 0.641 0.420 1.079 . 762 . 77982 .61951 | 1.26151 . 5737 
130° 0.635 0.429 1.103 . 756 . 76804 . 63687 | 1.29980 .5519 
140° 0. 629 0.435 1.128 .752 . 75409 .65057 | 1.33469 .5359 
150° 0.623 0.439 1.148 .749 . 73988 .66056 | 1.36757 .5253 
160° 0.617 0.442 1.169 .747 .72515 .66751 | 1.39906 .5199 
170° 0.610 0.444 1.193 . 746 . 70974 .67160 | 1.42951 .5196 
180° 0.601 0.445 1.210 .745 . 69396 .67297 | 1.45991 . 5246 

B=2, @=0.50 B=2, a=075 

7] "'§ 2 as WN € 2 T WVo 

0° . 00000 00000 | .00000 | 1.0000 . 00000 . 00000 00000 | 1.0000 
10° . 18669 01668 . 18049 1. 0744 . 19522 . 01769 . 18451 1.1220 
20° . 38323 06946 | .36583 1.1191 . 41893 . 07806 . 38224 1.2192 
30° . 56714 15508 . 54581 1.1333 . 64602 . 18411 . 58227 1.2847 
40° . 72148 26280 | .71267 1.1213 . 85169 . 32798 . 77497 1.3190 
50° . 83868 37948 . 86219 1.0908 | 1.01924 .49512 . 95383 1.3271 
60° . 91939 49412 . 99326 1.0492 | 1.14225 .67017 | 1.11575 1.3162 
70° . 96894 59963 | 1.10692 1.0032 | 1.22227 .84088 | 1. 26031 1.2938 
80° . 99413 69261 | 1.20527 .9575 | 1.26510 .99945 | 1.38874 1. 2663 
90° 1.00131 77225 | 1.29075 .9153 | 1.27784 | 1.14191 | 1.50308 1. 2383 
100° . 99562 83918 | 1.36577 .8781 | 1.26713 | 1.26698 | 1.60563 1.2133 
110° . 98088 89467 | 1.43243 .8471 | 1.23838 | 1.37493 | 1.69862 1.1932 
120° .95977 94015 | 1.49259 .8227 | 1.19570 | 1.46677 | 1.78411 1.1799 
130° . 93522 .97669 | 1.54802 .8051 | 1.14421 | 1.54358 | 1.86434 1.1743 
140° .90498 | 1.00602 | 1.59928 .7945 | 1.07901 | 1.60653 | 1.93967 1.1776 
150° .87325 | 1.02833 | 1.64839 .7909 | 1.00815 | 1.65630 | 2.01304 1.1906 
160° .83929 | 1.04435 | 1.69616 . 7946 .92987 | 1.69316 | 2.08549 | 1.2144 
170° .80283 | 1.05397 | 1.74305 . 8060 .84355 | 1.71590 | 2.15770 | 1.2499 
180° .76431 | 1.05731 | 1.79054 . 8255 .74951 | 1.72402 | 2.23190 | 1.2991 








where ¢ is defined in Eq. (9). The position coordi- 


nates are obtained by integrating 
dx = v cos 0 dt = vu" cos 6 dé/g(8 — cos A) 
and dy = w sin 6 dt = v) sin 0 d0/g(8 — cos @) 


Introducing dimensionless variables and integrating 
leads to 


(33) 
(34) 


g=o-0 
n = log (h/ho) 


(IV) Maneuvers With Constant Angular 
Velocity and Net Thrust Factor 


If the airplane follows a path with constant angular 
velocity, then the assumption of constant net thrust 
factor leads to simple formulas for the velocity, coordi- 
nates, and time along the flight path. Let d6@/dt = w, 
a constant, then @ = wt and 


(35) 
(36) 


dv/dt = g(a — sin wt) 


vw = g(B — cos wt) 
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The first equation can be integrated to give 


v= %t+ ¢ af _ | (37) 


or, in dimensionless form, 


cos wt — 1 


® 


v/t = 1+ 060+ aod — acosé (38) 
where o = g/vqw is the ratio of the acceleration of 
gravity to the initial angular acceleration. With 


Eq. (37), Eq. (36) gives a relation for the load factor 
8 as a function of time, net thrust factor, angular 


velocity, and initial speed. The coordinates are 


TABLE 


IN 


x 


PLANE 
y= [ 7 sin wt dt 
JV 


(1/w) (vo + g/w + gat) sin wt + 


VERTICAL 


v cos wt di, 
v7 WV 


which integrate to 


THE 


Xx 


(ga/w*)(cos wt — 1) — (g/2w)t — (g/4w*) sin 2ut 
y = (1/w)(% + g/w + gat)(1 — cos wt) + 
(ga/w") sin at — (g/2w?) sin? wf 


When the net thrust factor is zero, then v and y are 
has, in addition to its 


In 


periodic. The coordinate x 
periodic component, a translation term —gt/2w. 
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Flight paths for constant load factor 6 and constant net thrust factor a 







































































(8 = 3; a = 0, 0.25, 0.50, 0.75) 
B=3, a=0 B=3, a= 0.25 

g € 7 r WY € vj r WN% 

0° 0.000 0. 000 0. 000 1.000 .00000 | .00000 | .00000 | 1.0000 
10° 0. 082 0. 007 0. 073 . 996 .08807 | .00773 | .08778 | 1.0143 
20° 0.156 0.029 0.145 .985 .17333 | .03055 | .17483 | 1.0135 
30° 0. 221 0. 060 0.213 .970 . 24992 . 06617 . 25878 .9992 
40° 0.272 0.099 0.276 . 949 .31379 | .11077 . 33778 .9737 
50° 0.311 0.140 0. 336 .927 .36320 | .15999 .41067 . 9406 
60° 0. 337 0.180 0.391 . 906 . 39835 . 20997 .47700 . 9029 
70° 0.354 0.217 0.442 . 885 .42077 | .25779 | .53684 . 8638 
80° 0.362 0.250 0.487 . 866 .43261 .30160 | .59063 . 8254 
90° 0. 364 0.278 0.528 . 850 .43610 | .34045 | .63901 . 7893 
100° 0.364 0.301 0.571 . 836 .43323 | .37407 | .68272 7567 
110° 0.360 0.321 0. 605 . 824 .42565 | .40257 | .72249 . 7283 
120° 0.355 0.337 0. 642 . 814 .41465 | .42626] .75902 7044 
130° 0. 349 0.350 0. 675 . 806 .40143 | .44523] .79297 . 6851 
140° 0.345 0.359 0.720 . 801 .38598 | .46084 | .82485 . 6705 
150° 0. 335 0.366 0.765 .796 .36939 | .47244] .85526 . 6607 
160° 0.328 0.371 0.772 .793 .35189 | .48061 . 88465 . 6556 
170° 0.320 0.374 0.801 .791 .33367 | .48550 | .91349 . 6554 
180° 0.313 0.375 0. 833 .791 .31489 | .48714]| .94221 . 6601 

B=3, a=0.50 B=3, a=0.75 

] é 7 r WM € ” tr WM 

0° .00000 | .00000 | .00000 | 1.0000 . 00000 00000 | .00000 | 1.0000 
10° . 09001 .00796 | .08874 | 1.0366 . 09202 00820 | .08972 | 1.0594 
20° .18099 | .03237 | .17867 | 1.0583 .18911 03431 .18263 | 1.1051 
30° . 26622 | .07208 | .26724 | 1.0653 . 28397 07857 | .27608 | 1.1357 
40° .34026 | .12382 | .35231 | 1.0592 . 36979 13863 | .36768 | 1.1519 
50° .39978 | .18320 | .43235 | 1.0427 .44152 21026 | .45557 | 1.1559 
60° .44372 | .24574 | .50654 | 1.0191 .49644 28854 | .53856 | 1.1502 
70° .47274 | .30770 | .57465 .9915 . 53400 36882 | .61608 | 1.1382 
80° .48856 | .36635 | .63689 . 9626 . 55515 44735 | .68809 | 1.1228 
90° .49335 | .42000 | .69374 . 9346 . 56173 52144 | .75490 | 1.1064 
100° .48925 | .46779 | .74586 . 9087 . 55587 58938 | .81704 | 1.0912 
110° .47816 | .50943 | .79393 . 8864 . 53964 65021 .87516 | 1.0787 
120° .46164 | .54496 | .83867 . 8682 . 51485 70351 .92995 | 1.0702 
130° .44139 | .57420 | .88079 . 8548 .48377 74860 | .98221 | 1.0666 
140° .41689 | .59869 | .92078 . 8466 . 44509 78702 | 1.03239 | 1.0689 
150° .39021 | .61739 | .95940 . 8438 . 40206 81715 | 1.08142 | 1.0775 
160° .36135 | .63085 | .99715 . 8467 . 35445 83934 | 1.12991 | 1.0934 
170° . 33062 | .63908 /|1.03460 . 8557 . 30262 85321 | 1.17854 | 1.1173 
180° . 29825 | .64190 /1.07230 . 8713 . 24682 85804 | 1.22803 | 1.1501 
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dimensionless variables £, », o, and 6, the coordinates 
become 


— = o(1 + o6 + aad) sin 6 + ao*(cos @ — 1) — 


Ss 


(o?/2)0 — (07/2) sin 26 (39) 
n = o(1 + 0+ aad)(1 — cos @) + 
ao” sin @ — o* sin? @ (40) 


(V) Discussion of the Numerical Calculations 
and the Application of the Tables for Calculating 
More General Flight Paths 


A comparison of the loop flown at constant load 
factor and net thrust factor with a loop flown at con- 
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stant centrifugal acceleration and net thrust factor is 
provided in Fig. 2. The values of the parameters 
chosen for this figure are 8 = K = 2anda='/s. The 
path of constant centrifugal acceleration has greater 
lift in the flight path angles of 0° to 90° and 270° to 
360° than the flight path of constant load factor, and 
less lift in the angles 90° to 270°. The evidence of this 
difference is easily seen in the two loops of Fig. 2. The 
loop with constant centrifugal acceleration covers less 
altitude and range than the loop with constant load 
factor. Level flight is attained at & = 7.0 and 7 = 
—5.4 for the constant centrifugal acceleration loop 
and at € = 17.6 and » = —7.0 for the constant load 
factor loop. 


TABLE 4 


Flight paths for constant load factor 8 and constant net thrust factor @ 


(8 = 4; a = 0, 0.25, 0.50, 0.75) 









































B=4, a=0 B=4, a=0.25 
e € 7 r WY € 7 r WN% 
0° 0.000 0. 000 0.000 1.000 .00000 | .00000 | .00000 | 1.0000 
10° 0. 058 0. 005 0. 048 .997 .04276 | .00565 .05840 | 1.0095 
20° 0.112 0.019 0.115 .990 . 11503 . 02028 .11649 | 1.0090 
30° 0.160 0. 042 0.170 .979 .16660 | .04428 . 17316 . 9993 
i 40° 0. 200 0.070 0. 222 . 964 . 21072 .07511 . 22753 .9817 
50° 0.231 0.101 0.270 . 948 . 24602 . 11032 . 27897 . 9584 
60° 0.254 0.133 0.315 .931 . 27213 .14750 | .32710 . 9310 
70° 0. 268 0.166 0.356 .913 . 28949 . 18458 . 37182 .9017 
80° 0.276 0.194 0. 393 . 898 . 29904 .21998 | .41321 . 8722 
90° 0.278 0. 219 0. 428 . 884 . 30195 . 25265 .45148 . 8437 
100° 0.275 0. 242 0.457 . 871 . 29943 .28195 | .48694 . 8173 
110° 0.271 0.261 0. 487 . 860 .29260 | .30761 .51995 .7937 
120° 0. 261 0.278 0.512 . 850 . 28241 . 32955 . 55087 .7734 
130° 0.254 0.291 0.547 . 842 26979 .34750 | .58006 . 7568 
140° 0. 243 0. 302 0.563 . 835 . 25491 . 36258 . 60788 .7441 
150° 0.231 0.310 0.593 . 831 . 23871 . 37393 . 63467 . 7355 
160° 0.220 0. 316 0. 608 . 827 . 22142 . 38201 . 66075 .7310 
170° 0. 207 0.319 0. 629 . 825 . 20332 .38686 | .68643 . 7307 
180° 0.194 0.320 0. 650 . 825 . 18467 . 38849 .71202 7349 
B=4, a=0.50 B=4,a=0.75 

8 € 7 r Wo € ? r Wo 
0° .00000 | .00000 | .00000 | 1.0000 . 00000 00000} .00000 | 1.0000 
10° . 05929 .00522 | .05883 1.0242 . 06016 00533 | .05926 | 1.0393 
20° .11840 | .02107 .11819 | 1.0386 . 12189 02191 .11992 | 1.0691 
30° .17380 | .04689 .17694 | 1.0432 . 18142 04968 .18082 | 1.0891 
40° . 22253 . 08097 . 23409 | 1.0391 . 23524 08735 .24089 | 1.0996 
50° . 26257 . 12095 .28888 | 1.0278 . 28068 13274 .29927 | 1.1023 
60° . 29298 . 16427 . 34083 1.0113 . 31607 18322 .35534 | 1.0985 
70° .31369 | .20855 . 38970 .9915 . 34080 23611 .40875 | 1.0903 
80° .32535 | .25186 | .43547 . 9702 . 35504 28908 .45937 | 1.0794 
90° . 32898 . 29273 .47829 .9491 . 35957 34023 | .50726 | 1.0676 
100° . 32575 . 33021 . 51839 . 9293 . 35542 38816 | .55262 | 1.0565 
110° .31682 | .36371 . 55612 .9117 . 34374 43191 . 59573 1.0472 
120° . 30322 | .39294 .59181 . 8972 . 32561 47085 .63693 | 1.0409 
130° . 28614 .41735 . 62584 . 8864 . 30246 50408 . 67661 1.0381 
140° . 26547 .43819 . 65857 . 8796 . 27378 53284 .71511 1.0398 
150° . 24263 . 45419 . 69038 . 8773 . 24156 55541 .75289 | 1.0465 
160° .21780 | .46578 | .72162 . 8797 . 20591 57203 .79032 | 1.0587 
170° .19135 | .47286 | .75267 . 8872 .16726 | .58237 .82785 | 1.0773 
180° .16362 | .47528 | .78388 .9001 . 12601 58595 . 86591 1.1025 
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(K = 1; a = 0, 0.25, 0.50, 0.75) 
K=1], a@=0O K =1, a=0.25 
8 t € ” WN r € 7] WY 
0° . 0000 . 0000 . 0000 1.0000 . 0000 . 0000 - 0000 1. 0000 
10° . 1734 .1718 . 0150 . 9700 .1775 . 1796 . 0158 1.0288 
20° . 3418 . 3289 . 0568 . 8864 . 3573 . 3515 . 0618 1.0272 
30° . 5003 . 4599 . 1176 . 7649 . 5344 . 5212 . 1394 0.9970 
40° . 6460 . 5595 . 1869 . 6263 . 7040 . 6560 . 2337 0.9423 
50° . 7762 . 6284 . 2552 .4895 . 8624 . 7579 . 3349 0.8703 
60° . 8902 . 6713 . 3161 . 3679 1.0075 . 8268 . 4330 0.7880 
70° . 9883 . 6948 . 3657 . 2682 1.1373 . 8684 . 5206 0.7029 
80° 1.0717 . 7053 .4043 .1915 1. 2528 . 8883 . 5943 0. 6203 
90° 1.1420 . 7078 . 4324 . 1353 1. 3542 . 8936 . 6533 0.5449 
100° 1.2011 . 7062 . 4522 . 0956 1.4433 . 8900 . 6988 0.4784 
110° 1.2509 . 7027 . 4659 2.0683 1.5218 . 8808 . 7331 0.4224 
120° 1. 2932 . 6984 .4751 . 0498 1.5913 . 8692 . 7581 0.3766 
130° 1.3295 . 6941 . 4813 . 0374 1. 6537 . 8564 . 7764 0.3412 
140° 1.3612 . 6900 . 4854 . 0292 1.7112 . 8434 . 7892 0.3149 
150° 1.3895 . 6864 . 4881 . 0239 1.7642 . 8298 . 7989 0.2978 
160° 1.4152 . 6828 .4897 | 4.0207 1.8151 . 8161 . 8053 0.2889 
170° 1.4397 . 6795 .4906 | 8.0189 1. 8658 . 8021 . 8090 0.2885 
180° 1.4632 . 6763 . 4909 .0183 1.9027 tort . 8104 0.2968 
K=1, @=0:50 K=!1, a@a@=075 
6 r € ” WN r é ? WX 
0° . 0000 . 0000 . 0000 1.0000 . 0000 0000 . 0000 1.0000 
10° . 1813 . 1878 . 0167 1.0747 . 1854 1962 .0179 1.1228 
20° . 3733 . 3916 . 0717 1.1211 . 3906 4289 . 0806 1. 2232 
30° . 5709 . 5937 . 1662 1.1364 . 6109 . 6793 . 1979 1.2953 
40° . 7686 . 7765 . 2941 1.1219 . 8348 9280 . 3725 1.3359 
50° .9611 .9271 . 4425 1. 0824 1.0752] 1.1504 . 5950 1.3463 
60° 1.1452 1.0385 . 6023 1.0239 1.3093 | 1.3298 . 8513 1. 3303 
70° 1.3177 1. tig7 . 7569 0.9539 1.5386] 1.4575 1.1242 1.2948 
80° 1.4778 1.1494 . 8987 0.8797 1.7644} 1.5308 1.3966 1.2470 
90° 1.6249 1.1603 1.0224 0. 8069 1.9733 | 1.5522 1. 6533 1.1949 
100° 1.7597 1.1514 1.1265 0.7401 2.1829 | 1.5338 1. 8926 1.1451 
110° 1. 8838 1.1288 1.2118 0. 6824 2.3739 | 1.4776 2.1045 1.1029 
120° 1.9987 1.0970 1.2801 0. 6358 2.5636 | 1.3912 2.2905 1.0743 
130° 2.1067 1.0594 1. 3338 0.6014 2.7495 | 1.2780 | 2.4529 1.0615 
140° 2.2094 1.0170 | 1.3765 0.5814 2.9351 | 1.1383 | 2.5957 1. 0689 
150° 2.3099 . 9693 1.4096 0.5729] 3.1243 .9674 | 2.7078 1. 0806 
160° 2.4103 .9156 1.4340 0.5806 3.3216 . 7636 2.8072 1.1673 
170° 2.5133 . 8568 1.4501 0.6058 | 3.5338 .5156 | 2.8722 1.2719 
180° 2.6228 . 7897 1.4565 0.6511 3.7685 .2055 | 2.9039 1.4641 
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Graphs of the half loop—the flight-path angle range 
0° to 180°—are shown in Figs. 3 and 4 for constant 
load factor and in Figs. 7 and 8 for constant centrifugal 
acceleration. The dimensionless times for these flight 
paths are shown in Figs. 5 and 6 and 9 and 10. The 
tabulations of the variables ¢, &, 7, 7, v/v» for values of 
the net thrust factor a = 0, '/4, '/2, */, and for 8 = 
2,3,4and K = 1, 2, 3 are provided in Tables 1-7. The 
data for ¢, n, 7 for the constant-velocity climb with 
constant load factors of 8 = 2, 3, 4+ are given in Table 8 
and plotted in Figs. 11 and 12 for 0 < @ < 90°. 

More complicated flight paths may be estimated by 


use of these tables. For example, if the load factor 


must vary throughout a maneuver, an approximation 
to the true flight path may be found by considering the 
load factor to vary in a sequence of step functions ap- 
proximating the true variation. By maintaining the 
velocity and angle continuous at each point of dis- 
continuity of the load factor, we find the correct value 
of v to use in determining the actual velocity and co- 
ordinates from the dimensionless variables v/v, Af, and 
An along the flight path. To illustrate the method of 
computing the coordinates, we assume that the aircraft 
flies at a = '/, and B = 2 for 9 = O° to 20° and at 
a = '/,, 8 = 3 from 20° te 30°. Let the initial ve- 
locity be 100 mph; then we see from Table 2 that 
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TABLE 6 
Flight paths for constant centrifugal acceleration K and net thrust factor a 1 
(K = 2; a = 0, 0.25, 0.50, 0.75) - 
K=2,a=0O K=2, @=025 
] r é 7] WM r é ” WY, : 
. 
0° 0. 0000 0. 0000 0.0000 1.0000 . 0000 0. 0000 . 0000 1.0000 S 
10° 0.0871 0. 0864 0.0076 0.9924 . 0876 0.0875 . 0078 1.0143 b 
20° 0.1728 0.1677 0.0285 0.9702 .1918 0.1738 . 0343 1.0135 
30° 0.2560 0.2396 0.0752 0.9352 . 2640 0. 2543 . 0683 . 9985 la 
40° 0.3357 0.2993 0.1043 0. 8896 , 3501 0.3244 . 1169 .9707 tl 
50° 0.4111 0.3451 0. 1502 0.8364 4335 0.3815 1728 . 9329 (2 
60° 0.4816 0. 3780 0. 1968 0.7788 . 5097 0.4234 .2319 . 8877 
70° 0.5469 0. 3997 0.2411 0.7196 . 5850 0.4510 . 2907 . 8384 
80° 0.6071 0.4096 0.2812 0.6615 . 6599 0.4654 . 3465 . 7876 
90° 0.6625 0.4128 0.3160 0. 6066 . 7268 0.4690 . 3969 . 7382 
100° 0.7131 0.4102 0. 3454 0.5561 . 7888 0.4649 4413 .6917 
110° 0.7597 0.4040 0. 3694 6.5112 . 8473 0.4541 4791 . 6499 
120° 0.8026 0. 3950 0. 3884 0.4724 . 9023 0.4439 . 5106 . 6137 
130° 0.8424 0. 3847 0.4033 0.4398 . 9542 0.4225 . 5361 . 5841 
140° 0.8796 0. S735 0.4145 0.4135 | 1.0031 0.4036 . 5563 . 9612 
150° 0.9147 0. 3620 0.4226 0. 3934 1.0516 0.3824 -otLO . 5457 
160° 0.9484 0. 3503 0. 4282 0.3791 1.0972 0.3595 . 5824 ~ 0310 
170° 0.9811 0. 3384 0.4313 0.3707 1.1458 0.3348 . 5889 Doll 
180° 1.0132 0. 3266 0.4323 0. 3679 1.1935 0.3082 .5911 . 5448 
K=2, a2=0.50 K=2, a=075 
6 r é 7) WV Tt é 7 WY 
0° 0.0000 0.0000 0.0000 1.0000 0. 0000 0.0000 0.0000 1.0000 
10° 0.0889 0.0902 0.0073 1. 0367 0. 0899 0.0923 0.0081 1.0596 
20° 0.1805 0.1783 0. 6311 1.0588 0.1844 0.1912 0. 0349 1.1060 
30° 6: 2733 0.2723 0.0746 1. 0660 0.2824 0.2908 0.0818 1.1381 
40° 0.3660 0. 3532 0.1312 1.0592 0.3826 0.3812 0.1450 1.1558 
50° 0.4578 0.4212 0.1992 1.0404 0.4839 0.4679 0. 2301 1.1603 
60° 0.5474 0.4740 0.2743 1.0119 0.5850 0.5347 0.3289 1.1534 
70° 0.6343 0.5104 0. 3526 0.9767 0.6851 0.5831 0.4296 1.1379 
80° 0.7179 0.5309 0.4296 0.9379 0.7835 0.6142 0.5382 1.1167 
90° 0.7979 G. 35310 0.5027 0.8983 0.8885 0.6199 0. 6423 1.0931 
100° 0.8746 0.5318 0.5700 0. 8603 0.9743 0.6153 0.7460 1.0701 
110° 0.9482 0.5160 0.6297 0.8261 1.0669 0.5876 0. 8387 1.0502 
120° 1.0189 0.4922 0.6811 0.7974 1. 1577 0.5470 0.9242 1.0365 
130° 1.0873 0.4611 0.7257 O. Tita 1.2478 0.4943 1.0004 1.0303 
140° 1.1547 0.4246 0.7612 0.7625 1.3544 0.4288 1. 0664 1.0339 
150° 1.2208 0.3833 0.7907 0.7569 1.4282 0. 3516 1.1315 1.0395 
160° 1. 2868 0.3376 0.8121 0.7620 1.5215 0. 2628 1.1626 1.0804 
170° 1.3542 0. 2879 0.8252 0.7783 1.6176 0.1594 1.1891 1.1278 
180° 1.4238 0. 2338 0. 8304 0. 8069 1.7188 0. 0426 1.1924 1.2100 
for 6 = 20° (146.7) ?(0.35141) -- 
x= =e = 235.2 ft 
v/Vo = | 0274 (32.16) 
or v = 102.74 mph. For 8 = 3, a = '/4, and 6 = y = 41.4 it 
20° (Table 3), From Table 3 for 8 = 3 we find that the increments of 
v/v = 1.0135 the coordinates for 6 = 20° to @ = 30° are 
Since these two values of v must agree, we choose v = AE = 0.24992 — 0.17333 = 0.07659 
101.37 for the 8 = 3 portion of the flight path. The An = 0.06617 — 0.03055 = 0.03562 
position coordinates for @ = 20° from Table 2 are From this, we obtain 
Since vy = 100 mph or 146.7 ft/sec, then os (32.16) ——— 
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Ay = 24.5 ft 
The coordinates at 6 = 30° then become 


x = 287.9 ft y = 65.9 ft 


For ease in calculating the example, values of 8 were 
chosen such that the tables could be used directly. 
Solutions for intermediate values of 8 could be found 
by graphical or polynomial interpolation of the tabu- 
lated values. Lagrange’s interpolation formula for 
the polynomial which passes through the three points 
(x1, V1), (X2, Yo), and (x3, V3) is 


(x) = yiLi(x) + Volo(x) + V3L3(x) 


where 
Li = (x — X2)(x — X3)/(x%1 — X2)(x%1 — X3) 
Lo = (x — %)(X — X3)/(X2 — %1)(X2 — X3) 
L3; = (x — x1)(x — xX2)/(%3 — 1) (x3 — Xe) 
For interpolating the three tables of n(@) for 8 = 2, 3, 


and 4 we obtain 


n(8,8) = n(@)\g2 (1/2)(8 — 3)(8 — 4) — 
n(@) 3-3(0 —_ 2)(p8 — 4) + 

n(6) 8-4 (1/2)(8 — 2)(8 — 3) 

The other variables can be interpolated by similar 


relations. 












































TABLE 7 
Flight paths for constant centrifugal acceleration AK and net thrust factor a 
(K = 3; a = 0, 0.25, 0.50, 0.75) 
K=3, a =O K=3, @ =0.25 

6 t é ” WV, : é ” WV, 

0° . 0000 . 0000 . 0000 1.0000 . 0000 . 0000 . 0000 1.0000 
10° . 0582 . 0577 .0051 . 9949 . 0592 . 0585 . 0051 1.0095 
20° - 1197 .1126 . 0197 . 9801 .1181 .1160 . 0204 1.0091 
30° . £721 . 1621 . 0428 . 9563 1775 . 1691 . 0454 . 9989 
40° . 2269 . 2044 . 0722 . 9250 . 2352 . 2160 . 0778 . 9804 
50° . 2798 . 2383 . 1060 . 8877 . 2910 . 2546 .1165 . 9547 
60° . 3303 . 2634 .1417 . 8465 . 3453 . 2839 . 1586 . 9237 
70° . 3781 . 2800 iTS . 8031 . 3976 . 3042 . 2010 . 8891 
80° .4237 . 2893 .2118 . 7592 .4479 . 3156 . 2444 . 8529 
90° . 4666 . 2922 . 2433 . 7165 . 4966 . 3207 . 2846 . 8168 
100° .5071 . 2898 . 2714 . 6762 . 5416 . 3160 . 3218 . 7821 
110° . 5452 . 2833 . 2957 . 6393 . 5883 . 3072 . 3544 . 7503 
120° . 5815 2137 . 3160 . 6066 . 6321 . 2940 . 3830 . 7222 
130° . 6160 . 2622 . 3328 . 5783 . 6745 . 2685 . 4070 . 6987 
140° . 6489 . 2489 . 3460 . 5550 T1511 . 2578 . 4266 . 6804 
150° . 6806 . 2348 . 3559 . 5369 . 7547 ~ 2377 . 4417 . 6678 
160° tano . 2199 . 3628 . 5238 . 7937 . 2128 . 4526 . 6611 
170° . 7418 . 2048 . 3669 . 5160 . 8310 . 1883 . 4593 . 6607 
180° .7718 . 1895 . 3682 . 5134 . 8694 . 1633 .4616 . 6670 

K=3, a =0.50 K=3, a = 0.75 

8 r é 7 | “%| é ” WV, 

0° . 0000 . 0000 . 0000 1.0000 . 0000 . 0000 . 0000 1.0000 
10° . 0587 . 0594 . 0053 1.0242 . 0615 . 0603 . 0054 1.0393 
20° .1191 .1194 . 0215 1.0386 . 1247 . 1228 . 0222 1.0695 
30° .1798 . 1765 . 0480 1.0435 .1971 . 1845 .0511 1.0900 
40° . 2405 . 2282 . 0862 1.0371 . 2552 . 2419 .0912 1.1013 
50° . 3007 . 2722 . 1280 1.0268 . 3210 . 2916 . 1412 1.1042 
60° . 3598 . 3069 .1774 1.0079 . 3863 soar . 1991 1.0998 
70° S177 . 3319 . 2297 . 9844 . 4509 . 3617 . 2623 1. 0899 
80° . 4743 . 3454 . 2827 . 9581 .5145 . 3793 . 3280 1.0764 
90° . 5294 . 3499 . 3343 . 9310 . 5769 . 3852 . 3941 1.0612 
100° . 5826 . 3457 . 3831 . 9046 . 6379 . 3796 4585 1.0461 
110° . 6372 . 3338 .4279 . 8804 . 6980 . 3636 .5188 1.0332 
120° . 6853 3153 . 4675 . 8599 . 7569 . 3376 . 5742 1.0239 
130° . 7350 . 2912 . 5022 . 8441 . 8153 . 3027 . 6249 1.0195 
140° . 7838 . 2624 . 5312 . 8341 . 8734 . 2598 . 6672 1.0225 
150° . 8321 . 2295 .5550 . 8297 .9315 . 2094 . 7026 1.0330 
160° . 8805 . 1930 5711 . 8345 . 9904 . 1520 . 7318 1.0505 
170° . 9294 . 1532 . 5815 . 8463 1.0509 . 0881 . 7463 1.0835 
180° . 9793 .0951 . 5853 . 8667 1.1132 . 0178 7535 1.1260 
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TABLE 8 


Flight paths for constant velocity and constant load factor 


(8 = 2, 3, 4) 


















































B=2 B=3 B=4 
6 é ”? r € 7] T é ” T 
0° 0.0000 0.0000 0.0000 0. 0000 0.0000 0.0000 0.0000 0.0000 0.000 
10° 0.1745 0.0151 0.3490 0.0885 0. 0086 0. 263 0.0595 0.0050 0.234 
20° 0.3369 0.0586 0. 686 0.1699 0.0298 0.519 0.1169 0.0199 0. 466 
30° 0.4814 0.1258 1.005 0. 2464 0.0649 0.770 0.1674 0.0437 0.691 
40° 0.6029 0.2103 1.301 0.3109 0.1107 1.009 0.2129 0.0751 0.911 
50° 0.6985 0.3054 1.572 0. 3633 0.1643 1.236 0. 2493 0.1125 1.122 
60° 0.7778 0.4055 1.825 0.4088 0.2231 1.456 0.2798 0.1542 1.327 
70° 0.8163 0.5056 2.038 0.4313 0.2844 1. 653 0. 3023 0.1983 1.524 
80° 0.8447 0.6024 2.241 0.4497 0.3459 1. 846 0. 3127 0. 2433 1.709 
90° 0. 8472 0. 6932 2.418 0.4552 0.4055 2.026 0.3152 0.2877 1.886 
| 
cS) | 
— 
| / 
/ y | 
f F 4 + 
p=4 / B=3 B=2 
a 
7 OB Lt . 
‘IG ‘li -p% é - -) i -nsi ass ti 4 oO ig S 
Fic. 10. Flight-path angle @ vs. dimensionless time 7 for flight eo, 12. Piletit pathe ot content velocity ® and constant load 


paths at constant centrifugal acceleration AK and constant net 
thrust factora. K = 3; a = 0, 0.25, 0.50, 0.75. 
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Fic. 11. Flight-path angle 6 vs. dimensionless time 7 for flight 
paths at constant velocity v and constant load factor 8. B = 2, 
3, 4. 


factor B. B = 2, 3, 4. 
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The Newtonian Approximation in Magnetic 
Hypersonic Stagnation-Point Flow' 


PAUL S. LYKOUDIS* 


Purdue Unwersity 


Summary 


The hypersonic flow of an electrically conducting fluid around 
the stagnation region of a sphere carrying a radial magnetic field 
is examined. 

By assuming a Newtonian pressure distribution and constant 
density, the differential equation of the inviscid flow is integrated 
and a simple closed-form solution is obtained. 

It is found that the ratio of the stand-off distances of the shock 
wave for the magnetic and nonmagnetic cases does not depend ex- 
plicitly on the magnetic parameter S, (ratio of the ponderomotive 
force to the free-stream inertia force) nor on the density ratio e (the 
value at the free stream divided by the value behind the shock 
wave) but on the product S,v ‘e, at least for values of « between 
1/5 and 1/20. 

The velocity gradient on the body is also calculated and the 
ratio of the magnetic to the nonmagnetic case is shown to depend 


on the parameter S;V e. 


The case of cylindrical shocks is also examined; the same 
general conclusions are drawn. 
Symbols 
B = magnetic inductance 
C = coefficient of proportionality 
C, = pressure coefficient defined in Eq. (3) 
J = current density 
m = parameter defined in Eqs. (16) and (23) 
P = positive root of the quadratic (15) or (22); physically, 
the inviscid nondimensional velocity gradient 
p = pressure 
Q = negative root of the quadratic (15) or (22) 
R= radius of curvature 
S, = ratio of magnetic force over inertia force defined in 
Eq. (8) 
U = free-stream velocity ahead of the shock wave 
u = velocity in the x direction inside the shock wave 
x = coordinate along the shock wave (see Fig. 1) 
XY; = coordinate inside the shock wave (see Fig. 1) 
4 = coordinate in the radial direction (see Fig. 1) 
P4 defined in Eq. (33) 
z defined in Eq. (12) 
A = nondimensional stand-off distance [= (Rs — Rp)/R] 
€ = density ratio across the shock wave ( = p,,/p) 
c defined in Eq. (34) 
p = mass density 
o = coefficient of electrical conductivity 
Subscripts 
b refers to the radius of curvature of the body 
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e refers to the edge of the boundary layer 
0 refers to the nonmagnetic case 
§ refers to the radius of curvature of the shock 
refers to conditions prevailing ahead of the shock wave 
cyl refers to the cylindrical case 
sph. refers to the spherical case 


A bar (~) denotes a nondimensional parameter 


(1) Introduction 


Ww" A vehicle with a blunt nose moves at hyper- 
sonic velocity, part of the kinetic energy asso- 
ciated with its motion is transformed, in the stagnation 
region, into thermal motion, a substantial fraction being 
used for the dissociation and to some extent ionization 
of the different species present. 

An ionized gas is of course a good electrical conductor, 
and hence it is possible in the presence of electromag- 
netic fields to impart forces whose relative magnitude 
may be considerable in comparison with the inertia or 
viscous forces. It is well known that in the presence of 
a magnetic field fixed relative to the moving body, and 
acting perpendicularly to the oncoming flow, these forces 
are decelerating. 

Qualitatively the following things are expected to 
happen: 

(a) Since the flow is retarded, a wider area is needed 
for the passage of the hot gas between the shock wave 
and the nose, and hence the stand-off distance of the 
shock wave will be increased. 

(b) The coefficients of skin friction and consequently 
heat transfer will be lowered since the free-stream ve- 
locity is diminished. 

In the last two years a large number of papers have 
appeared in the literature dealing with different aspects 
of the above problem (see, for instance, references 1-9). 
The present work is concerned with the phenomenon 
(a) above. In reference 8 a numerical solution is given 
as a result of the digital integration of the basic con- 
servation equations for the determination of the stand- 
off distance of a spherical shock wave. Numerical solu- 
tions in general give a specific answer for a particular 
set of flight conditions whereas they fail to reveal the 
meaning and the number of independent parameters 
describing a physical situation. 

An attempt will be made here to give a closed-form 
answer to this problem by making use of the Newtonian 
approximation. It will be seen that although the shock 
stand-off distance depends on the ratio S, of the pon- 
deromotive force to the free-stream inertia force and the 
ratio e« of the mass densities across the shock wave, the 
important parameter is the product S.ve. The same 
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Fic. 1. Nomenclature for the basic geometry 


is true for the velocity gradient at the wall. The closed 
form solution obtained will be compared with the 
numerical, exact findings of reference 8 given only for a 
specific density ratio across a spherical shock, and it 
will be seen that the agreement is good. The case of 
cylindrical shocks will also be considered. 


(2) Basic Equations and Approximations 


The following approximations will be made in this 
work: 

(a) The magnetic field is constant throughout the 
region between the shock wave and the surface of the 
vehicle; it remains also always perpendicular to the 
oncoming stream lines. 

(b) The induced electric current does not distort the 
applied magnetic field. 

(c) The shock waves are always spherical or cylindri- 
cal depending on the geometry considered. 

(d) The pressure distribution in the stagnation 
region is not altered by the presence of the magnetic 
field. A Newtonian pressure distribution with constant 
density will be assumed. 

(e) The coefficient of electric conductivity is a scalar 
and remains the same everywhere. 

(f) Viscous effects are negligible. 

Assumption (a) is difficult to realize in practice since 
the magnetic field will originate from a coil located co- 
axially with the vehicle in the neighborhood of the 
nose. One expects it to diminish rapidly in the outward 
direction. For instance the rate of reduction far away 
from a magnetic dipole is inversely proportional to the 
third power of the distance, whereas for the experimen- 
tal apparatus described in reference 9 the power was 
found to be approximately 3.61. It was shown in 
reference 4 that if the assumption is made that the 
constant magnetic field assumed in the calculations is 
an average one calculated inside the shock region for the 
specific law of variation of the actual magnetic field, 
then the final result is equally acceptable. 

Assumption (b) is true as long as the magnetic 


Reynolds number based on the body radius is much 
smaller than one; for most of the flight conditions in 
which the present theory applies this is always the case 

In the neighborhood of the stagnation point, assump- 
tion (c) is expected to be valid.* The theoretical and 
numerical solutions described in references 3, 7, and 8 
have shown that the pressure distribution changes very 
little with the imposed magnetic field at least for a 
moderate variation of the appropriate magnetic non 
dimensional parameter. The success and simplicity of 
the Newtonian approximations in hypersonic flow are 
well known. The basic premises of this theory and its 
consequences are described very satisfactorily in ref- 
erence 1(). For the calculation of the nonmagnetic 
shock stand-off distance this assumption gives lower re- 
sults than more exact numerical solutions. Neverthe 
less we expect a priori to find the ratio of this distance 
for the magnetic and nonmagnetic cases very close to 
the true value, since the same approximations are 
made for the calculation in the numerator and de- 
nominator. 

Finally, assumption (e) should be understood to be 
valid as long as the electrical conductivity is calculated 
at an average temperature. The present calculations 
are confined to the inviscid region near the stagnation 
point, where this assumption becomes more plausible 
On the other hand the assumption is well justified in 
the case of a seeded fluid within the temperature range 
where the electrical conductivity is insensitive to tem- 
perature variations. Its scalar character will be con- 
served as long as the collision frequency of the particles 
is much higher than the cyclotron frequency. In this 
case the charged particles are not free to spiral around 
the magnetic lines and as a result they lose the sense of 
a preferential direction. 

In order to facilitate the reading of the present paper, 
use will be made, in so far as possible, of the nomencla- 
ture of reference 10. As shown in Fig. 1 the coordinate 
x denotes the distance of a streamline measured on the 
shock wave from the stagnation axis. The quantity x; 
measures in the direction « the coordinate of any point 
inside the shock layer lying on the same streamline. 

The ‘‘Newton-Busemann”’ pressure laws for the cases 
of a sphere and a cylinder are, respectively: 


l C | x," l ( ™ 2) ( | ) 
»=1- — re — xX 
:™ R?2 3R2°- 
] - x1? l , a ) 
5 oe = l ~ R2 - oRe — 29 (2) 
where Cp = (pb — po)/(1/2)p.U? (3) 





The above distributions contain the Busemann cen- 


* The Newtonian approximation is successful in hypersonic flow 
because the velocity after incidence remains tangent to the wall, 
since it is bounded by the thin shock layer. In this manner the 
momentum associated with the normal component of the inci 
dent free-stream velocity is balanced by the wall pressure. Now 
this last component does not interact with the radical magnetic 
field, since it is essentially parallel to it and therefore the pressure 
remains unaltered at the wall as long as the shock stand-off dis 
tance does not become too high. 
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trifugal pressure correction; reference 10) gives their 
derivation in detail. 

Behind the shock wave the equation of conservation 
of momentum is given as follows: 


pu(du/dx;) = —(Op/0x) + (J XK B)x, (4) 
This is a generalized Bernoulli equation written along 
the streamline x,. By virtue of the assumptions (a), 
(b) and (e) the current induced in a direction perpen- 
dicular to the plane formed by the coordinate x and 
the magnetic field B is given from Ohm's law: 

J = ouB (5) 

Substituting this expression into Eq. (4) we obtain: 
pu(du/dx,) = —(Op/0x) — oB*u (6) 
Nondimensionalizing the velocity u by using the free- 
stream velocity U and the distances by using the radius 
of curvature of the shock, R,, we can write the above as 


follows: 
pu(di/d%,) = —(Op/U?0%) — (cB*ui/UR,) (7) 
At this point we define the ratio of the ponderomotive 
forces over the inertia forces as follows: 
S; = cB?R,/p.U (S) 
By using the above definition and definition (3) in Eq. 
(7) we obtain: 
adi/di%, = «} —(1/2)(0C,/d0%) — S.a} (9) 
where ¢ is the density ratio across the shock wave: 
€ = ?p./p (10) 
Integration of the basic Eq. (9) will yield the velocity @ 
as a function of x, 4, S;, and e. A simple application of 
the equation of conservation of mass will then yield the 
shock stand-off distance. A detailed analysis of this 
method for the nonmagnetic case, along with a discus- 
sion of the different assumptions involved, is given in 


reference 10. 


(3) The Sphere 


Substitution of the pressure distribution given by Eq. 
(1) into Eq. (9) gives: 


~ 
= 


u(di/dk,) = €}(8/3)%, — S,a} (11) 


This is a nonlinear homogeneous differential equation of 
the first order that may be solved by introducing the 
variable: 
2(21) = u ry (12) 
After separation of the variables Eq. (11) becomes 
dz di, : 


e| (8/32) — S,} — 2 zy 


Integration starting from the limit 7; = Z results in 


=\9 ae i—m 7, __ 1+m 
(ey (=2 "(2g)" os 
Zi 1—P i-@ 


where P and Q are the positive and negative roots of the 


quadratic 
s+ S.es — (8/3)e = 0 (15) 
P = —S.e + V (S.€)? + (32/3)e (15a) 
>) 
— Sie — V (Sie)? + (32/3) , 
O=- ~ _ = (15b) 
x 9 


and m is defined as follows: 


~f l . 
m = P-O a + (32/35. (16) 
c Oa / Ws €) 

It is not difficult to see by inspection of Eq. (15) that 

the absolute value of the root Q is always greater than 

P so that the quantity m is always a positive number. 

Furthermore, for zero magnetic field S, is zero and m is 

equal to zero also, whereas for an infinite S,, mm asymp- 

totically approaches the value one. It will be shown 

later that this parameter plays an important role in the 
final solution. 

From Fig. | it is seen that the following equation of 

mass continuity for the case of sphere may be written: 

p..(2ax dx) U p(2mx, dy)u 


Simplifying, using nondimensional quantities and inte- 


y *  Fdi 2 
: = ¢€ (17) 
R, x=0 FU 


For the calculation of the above integral it is necessary 
It is evident that 


grating, we find: 


to solve Eq. (14) for the velocity @. 
such a solution in closed form cannot be found. Instead, 
we differentiate Eq. (14) and by substitution in Eq. (17) 
it becomes possible to use the quantity z as the variable 
of integration. The result is surprisingly simple and 


is given as follows: 


’ 3- ) m 
—_ ~ f ( 5) dCs) 
R, (1 — P)'!-"(1 — Q)'** Js=P \z — P 


* The lower limit of integration can be understood by looking at 

















Eq. (14); for positive values of z, x becomes zero when zs = P 
3.0 
— 
“s° 2.0 
“| Dato of Ret. (7 
7) 10 20 
Ss 


Fic. 2. The stand-off distance for a spherical shock wave 
versus the magnetic parameter S, for various values of the den- 
sity ratio e. 
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Fic. 3. The stand-off distance for a spherical shock wave versus 
the parameter m 


The coordinate y is measured from the stagnation point 
at the body; integration between the distance R, — R, 


yields 


A, = (R, — R»)/R, = 


1 - 
€ se Q m ; 
(1 _— P)} m(] — Q)i+m {. (: ie S) dz (19) 


Closed form calculation of the integral (19) is possible 
only if the parameter m is a rational number. This is 
not very much of a limitation since for all practical pur- 
poses m may be more than adequately approximated by 
a ratio of two natural numbers. 

Values of the stand-off distance A, have been calcu- 


lated for four values of e('/s5, '/1, 1/15, 1/2) and seven 
values of m ("/s, 1/4, 1/3, 1/2, 5/5, 3/4, */g). For the case 
m = ( corresponding to zero magnetic field we obtain 


the well-known result!’: 
Aw = ¢/(1 + V8e/3) (20) 


Fig. 2 shows the ratio A,/A,o of the stand-off distances 
for the magnetic and nonmagnetic cases plotted versus 
S, for different density ratios «. In Fig. 3 the same 
data are plotted, this time versus the parameter m. It 
is seen that the dependence on ¢ has completely dis- 
appeared at least up to the highest calculated value of 
m (= 7/8). We draw the conclusion that the ratio 
A,/A. depends only on the quantity S,/e rather than 
on S, and ¢ explicitly [see Eq. (16)]. A physical argu- 
ment as to the meaning of this finding will be given in 
Section 5. 

The data presented in reference 8 for « = 1/10 ac- 
cording to the findings of reference 7 have been recom- 
puted in terms of the parameters S, and A,/A, rather 
than S, and A,/Aj». The agreement with the New- 
tonian approximation is perfect, as shown in Fig. 2. 
One can also verify by comparison with reference 8 that 
Fig. 2 is equally valid, by merely changing the coordi- 
nate subscripts s to 0. 

It should be pointed out that in the work of references 
6 and 8 the intensity of the magnetic field in the radial 
direction was varied inversely with the square of the 


distance away from the stagnation point; however, the 
solution obtained there was shown to be insensitive to 
different numerical values involving the coefficient of 
proportionality C of the inverse square law B, = C/r’. 


(4) The Cylinder 


For the cylindrical case the equation of conservation 
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Fic. 4. The stand-off distance for a cylindrical shock wave 
versus the magnetic parameter S, for various values of the density 
ratio e 


of momentum in its nondimensional form is the follow- 
ing: 

idi/d, = €)3%, — Sai} (21) 
Use was made of the pressure distribution given by Eq. 
(2). This differential equation is exactly the same as 
Eq. (11) except for the numerical coefficient character- 
izing the pressure distribution. As a result Eq. (14) is 
still valid, with the roots P and Q appearing there found 
from the following quadratic: 


2? + S,ez — 3e = 0 (22) 


—S.,e + V (S.€)? + 1Qe 


9 


99 
Q 9 (22b) 


The quantity m is now defined as 
on.’ | 
P-Q V1 + (12/S,%e) 


m= (23) 


The equation of conservation of mass for the cylindrical 
ase is written as follows: 


po. dx U = pdyu (24) 


A unit length was assumed in the direction parallel to 
the axis of the cylinder. 
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Using again nondimensional quantities and integrat- 


y dz n 
= ¢€ - (25) 
R, r=0 W 


Using the same substitutions and general method ap- 
plied in the case of sphere, we obtain the cylindrical 
shock stand-off distance 


A, = (R; — Ro)/R; 


ing, we find: 


rl 
J 
(] — P) 1—m 2(] — Q)Utm/2 P 
dz 
- ———-_ (26) 
(g — P)G+mi2(g — Q)U-m/2 
For the case m = () this integral reduces to the known 


expression 
Aw = (« Wi- 3e) cosh—!(1 V 36) (27) 


The integral (26) may always be calculated in a closed 
form provided m is a rational number. Such calcula- 
tions have been undertaken for two values of ¢« ('/; and 
'/.9) and four values of m (1/3, 1/s, 2/3, 3/4). 

The ratio A,/A,. has been plotted versus S, for two 
values of ein Fig. 4. Thesame data plotted versus m are 


given in Fig. 5. It is seen that the dependence on the 














2.4 
2.0) 

| | 

€ }_ —_ 

Os 2 ° 

dso 
5 
Lo 
0 05 10 


m 


Fic. 5. The stand-off distance for a cylindrical shock wave versus 
the parameter m. 


density ratio « is again small. For instance, the dif- 
ference between the values of A,/A,y at « = 1/5 and 1/20 
is about 2.5 percent for m = 0.5, whereas for m = 0.7 
(almost a limiting value for an actual experiment), it 
becomes 5.0 percent. These seem to be rather small 
differences compared to the 400 percent change in e. 
It is therefore in order to conclude that for all practical 
purposes in the cylindrical case also the ratio A,/ A, de- 
pends on the parameter m only. 


(5) The Velocity Gradients 


The positive roots of the quadratics of Eqs. (15) and 
(22) have a definite physical meaning. It can be seen 
by inspection of Eq. (14) that P is the value of z for 
which x = 0). 


= t/Zi, = uR,/Ux, = P (28) 


From this expression we conclude that the velocity 
gradient at the body is given as follows: 


0u/Ox, = (U/R,)P (29) 


In other words, P is the nondimensional velocity 
gradient at the body and is given by Eqs. (15a) and 
(22a) for the spherical and cylindrical cases respec- 
tively. 

The ratio P/ Py comparing the velocity gradients for 
the magnetic and nonmagnetic cases is given as follows: 


S, /2 7s IaX2 
P/Tiak. * — = '; 4. (> ¥) +1 (30) 
“ 2 8 
= . » 
(P/Po)cyi. = — 3 \. |. y(° \") +1 (31) 


For both geometries the quantity P /, is given in terms 
of m as follows: 


P/Py = VQ — m)/(1 + m) (32) 


Eqs. (30), (31), and (32) have been plotted in Figs. 6, 
7, and 8. 

The above results are not new. In references 1, 2, 
1, and 5 the same basic formula is given, although writ- 
ten in terms of another parameter which will be dis- 
cussed in the following section. 

For the spherical case and for e« = 1/10 we have the 
results for the velocity gradient of reference 8 up to a 
value S, ~ 5. Comparison of these data (after their 
transformation to values referred to the radius of the 
shock wave) shows that the data of Fig. 6 are higher by 
about 8.5 percent at S, = 5. The data also show that 
Fig. 6 coincides completely with the digital computa- 
tions of reference 8 if the coordinate subscripts are 
changed from s to 6. This substitution apparently al- 
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Fic. 6. The velocity gradient for a spherical shock wave versus 
the magnetic parameter 5S, and various density ratios e. 
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Fic. 7. The velocity gradient for a cylindrical shock wave ver- 
sus the magnetic parameter S, and various density ratios e. 


lows for the change of the pressure distribution in the 


presence of a magnetic field. 


(6) Discussion 


In reference 5 an approximate analysis was made of 
the shock stand-off distance for both steady and non- 
steady hypersonic stagnation flows. An application of 
Bernoulli's equation at the edge of the boundary layer 
had shown (see references 1, 2, or +) that the ratio of 
the two velocities for the magnetic and nonmagnetic 


cases is equal to: 
u, ¢ ‘Se , - 
-— + 1+ ( = = 7 (.3:3 ) 
Ugo ys ; 


where ¢ = oB?/p( UP. /R;) (34) 





The quantity in parentheses is the velocity gradient at 
the body calculated for the nonmagnetic case and is 
always proportional to the square root of « in the New- 
tonian approximation; it is not difficult to deduce that 
¢ is given as follows: 

Soph. = V3/8S.V (34a) 

. mace ‘ ‘ 

ta «BV s/¥2 (34b) 

An equation for mass conservation within the thick- 

ness of the shock layer yields: 

p(R »— Ro)ttay = p(Rs — Ro)ottay (35) 
Assuming that Eq. (33) is valid also for the average 
velocities and since A, and A,» are negligibly small com- 
pared to unity we may write: 

y ge /é / se /)\9 or 
A./Aw = 1/Z = (§/2) + V1 + (§/2)? = (36) 


The above order-of-magnitude argument shows that the 
ratio of the stand-off distances depends solely on ¢ 
which, in turn, has been shown to depend on S,+/e 
only. 

It is thus seen that use of the parameter ¢ (the ratio 
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of the ponderomotive forces to the inertia forces calcu- 
lated behind the shock) is a more natural parameter 
than S,, which contains the inertia forces calculated at 
the free stream. The approximation (36) predicts very 
high values as compared with the better solutions in the 
present work. A power-law approximation for the 
stand-off ratio A,/A,j versus 1/Z, in good agreement 
with Eqs. (19) and (26) for values of m <0.5, is given as 


follows: 
(Ay/Aso)spn. © (1/Z)'/? (37) 
(A,/Ago)ey1. © (1/Z)3/4 (38 
(7) Conclusions 
The following conclusions can now be drawn. A 


theoretical, closed-form calculation of the shock stand- 
off ratio A,/ Axo for spherical and cylindrical shocks is 
feasible within the limitations of a Newtonian nonmag- 
netic pressure distribution and constant electrical 
conductivity and density. 

It is shown that this ratio depends essentially on the 
parameter Z which expresses the ratio of the velocity 
gradients at the stagnation point with and without a 
magnetic field. The dependence of Z on the parameter 
¢, which represents the ratio of the ponderomotive force 
to the inertia force calculated in the shock layer, is 
given by Eqs. (33), (34), and (8). The above results 
are in good agreement with the limited results of exact 
numerical solutions in the literature. 

The following list of ratios for the magnetic and non- 
magnetic cases summarizes approximately the theoreti- 
cal results: 


(As/Aw)epn. & (1/Z)'? 

(A,./As)ey1. = (1/Z)*/* 

P/Po = e/Ua = Z 
G/ag = 2" 


9 


(Continued on page 562 
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Fic. 8. The velocity gradient for both spherical and cylindrical 
shock waves versus the parameter m. 
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Body Under Lifting Wing' 


C. F. CHEN* anp JOSEPH H. CLARKE** 


Brown Unwersity 


Summary 


An investigation is made of supersonic-aircraft configurations 
composed of a cambered body positioned a certain distance 
beneath an arbitrary lifting wing. The geometry of the wing is 
regarded as given and the geometry of the body may be given or 
Expressions for the drag and lift are obtained from 
such a 


optimum. 
reverse-flow considerations; these greatly 
study when interference cross flows must be cancelled. The 
drag advantage to be gained when a given body and wing as- 
Treated more extensively 


implement 


sume a given orientation is studied. 
is the variational problem of determining the optimum wing 
incidence and optimum body shape, for the given volume and 
length, to yield the minimum drag for prescribed lift. |. Numeri- 
cal results are provided to indicate the significance of the large 
number of parameters appearing in the problem. Of these, the 
gap between the wing and the body is found to be particularly 
It is found that at low gap moderate body distor- 
Drag reductions 


important. 
tions have a significant influence on the drag. 
of up to 44% relative to the case of no interference have been 
found at a Mach number of 2.24 in a configuration having a 
gap approximately equal to the maximum diameter of the body, 
and a wing chord of about three eighths of the length of the body 
Comparison is made with the conventional wing-body combina- 
tion including the effects of skin friction, and it is concluded that 
the advantage suggested by the preceding considerations is not 
appreciably diminished. Finally, it is shown that the con- 
figurations studied lead to bodies of fineness ratios much lower 
than are appropriate to conventional wing-body combinations. 
Tests were made on an arrangement consisting of a Sears- 
Haack body located under a lifting rectangular diamond-profile 
The Mach number was 1.6 and the Reynolds number 
It was found that the 


wing. 
was 9.17 X 10° based on the body length. 
measured lift developed on the wing due to the flow field of the 
body agrees very well with the theoretical value. Downstream 
of the impinging shock from the wing, flow separation was ob- 
served on-the exterior of the body but not in the interior. The 
separation is attributed not to the pressure rise across the shock 
but to the pressure field arising from the reflection from the 
body of the shock-induced cross flow. Further observations 
suggest that the separation can be avoided by pitching the body 
or by kinking the body at the shock wave to accommodate the 


shock-induced cross flow. 


~ 
Symbols 

A = aspect ratio, or surface area of bodies 

Qn = Fourier coefficients of source distribution 
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B = VM? -—1 
bn = Fourier coefficients of vorticity distribution 
Cop = drag coefficient based on wing area 

Cr = lift coefficient based on wing area 
ri = chord 

D = drag 

' = lineal source intensity 

f = volumetric source intensity 
h = gap distance between wing and body axis 
i, E. i” P : . , ‘ 

7,0, J, J,* (= functions of n defined by Eqs. (43), (44), 

and (47 
T,, J,@ = functions of n, a,, and b, defined by Eq. (45 
i, j,k = unit vectors in the direction of x, y, 2 
L = lift 
l = body length 
M = Mach number of the free stream 
m = lift-curve slope 
N = unit outward normal vector to the section 
of A made by a plane x = const 
n = unit outward normal vector to A 
p = pressure 
gq = dynamic pressure of free stream 1/2pl” 
R = body radius 
= surface area or cross-sectional area 
= volume enclosed by A 

t = thickness ratio 

U = free-stream velocity 

V = perturbation velocity 

WwW = @V 

U, Vv, W = components of perturbation velocity 

Ve = Cartesian coordinates 

a = angle of attack 

r = lineal vorticity intensity 

0 = angular variable defined by Eq. (38 

= Lagrangian multiplier 

p _ density 

oy = defined by Eq. (3 

Q = volumetric distribution of vorticity 

re) = bound vorticity 

Subscripts 

0 = free-stream conditions 

l = body 

2,w = wing 

F = forward flow 

R = reverse flow 

comb = conventional wing-body combination 

iso = isolated system 

S-H = configuration with Sears-Haack body 

opt = configuration with optimum body 


Primes denote interference flow quantities 


(I) Introduction 


UCH HAS BEEN WRITTEN!~ in recent years on the 
M subject of favorable interference in supersonic 
flow. Attention has been given both to a number of 
specific arrangements which appear to offer possibilities 
for reducing the inviscid drag due to volume or lift, 
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and to the underlying theory concerning the nature of 
the forces arising from a general field of disturbance 
to a uniform incident stream. In reference 1, a con- 
figuration composed of a mutually interfering lifting 
wing and slender body located underneath was pro- 
posed and analyzed, for one special case. The inter- 
ference considered has three principal favorable fea- 
tures. The important upwash field of the body en- 
closing the required volume results in the appearance 
of an interference lift on the wing above, this lift appear- 
ing without wing incidence according to the mechanism 
of references 1, 2, and 3. Further, the overlapping flow 
fields can produce drag reductions in the manner 
suggested by the well-known two-dimensional biplane 
of reference 4. Finally, unlike the analogous interfer- 
ence in the biplane, the interference lift on the wing 
above cannot be counterbalanced by a negative lift 
on the body below, since the total lift on the body must 
vanish when it is closed. 


Fic. 1. Interfering configurations: (a) two-reflection case; (b) 
single-reflection case. 


Later, as an illustrative application of a general re- 
verse-flow relation,’ one of the present writers pre- 
sented certain expressions for the lift and drag of 
arbitrary arrangements of the wing and the fusiform 
body*; these expressions may require only limited 
knowledge, or no knowledge, of the interference flows 
needed to preserve the boundary shapes. Using 
these expressions, in the present paper we undertake to 
study the drag advantages to be gained both when a 
given body and wing assume a given orientation and, 
more extensively, when the said body is recontoured 
to its optimum shape for the given wing and orienta- 
tion. The body may be interpreted as a fuselage, a 
store, etc. 

The configuration analyzed then consists of a wing 
located above a fusiform body. Three possible formu- 
lations of the problem are of interest here: (1) wing 
represented by given singularities and body represented 
by given or optimum singularities, (2) wing with given 
geometry and body with given or optimum geometry, 

* By fusiform body is meant a body of fuselage or projectile 


shape which does not necessarily satisfy the smoothness property 
of a slender body as specified by Ward.® 


and (3) optimum wing and optimum body. In the 
first problem the required operations, such as opti- 
mization, are performed on the singularities represent- 
ing the configuration. The geometry of the con- 
figuration is then obtained by constructing the stream 
surface induced by the resulting singularity distribu- 
tion. This may be termed the indirect method. 
The optimum problem indicated has been treated by 
Lomax and Heaslet.’ In the second formulation, in 
which both the wing and the body are given geometri- 
cally, interference fields must be added to the isolated 
disturbance fields to preserve the boundary shapes. 
Operations must be performed on the summation of 
these fields. This may be referred to as the direct 
method. (Two numerical examples are presented 
in this paper to compare the results given by these 
two different formulations.) One may, however, 
choose to find the optimum body in the second formu- 
lation using the indirect method, thus shifting the 
difficulty of determining the interference field of the 
body to that of subsequently determining the body 
geometry. But if the reverse-flow relations are used 
to obtain the forces, there are interesting instances in 
which knowledge of the interference field of the body 
is not needed in the direct method. Evidently the 
third problem can be approached by either the indirect 
or the direct method. In the process of optimization, 
one may regard the horseshoe vorticity representing 
the lift distribution of the wing as seeking to develop 
the favorable interference lift by tending to cancel the 
upwash field of the body and, in addition, warping the 
wing in an optimal way. At the same time, the body 
geometry is disposed to obtain maximum advantage 
from the wing flow field and to affect the wing favor- 
ably, much as in the other problems. 

In this paper the second formulation is chosen, and 
interference flows are associated with both wing and 
body. It is believed that this formulation, in which 
the wing geometry is given, gives in addition an indica- 
tion of what results may be expected from the most 
difficult third formulation. Whereas it takes no ac- 
count of any advantages to be gained by warping the 
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Fic. 2. Drag polar for configuration Rl; A = 1.40, c/l = 
0.293. 
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wing and redistributing its thickness (advantages which 
are well understood for isolated wings), it does account 
for the additional wing vorticity connected with inter- 
ference lift, the effect believed to be the more important 
one. On the other hand, the first formulation, what- 
ever its intrinsic interest, cannot give any indication 
of the results to be expected from the third, since it 
excludes interference lift (and the decrease in the drag 
of the given wing vorticity due to body upwash turns 
out to be much less important). The choice made 
also provides an ancillary advantage in that it is the 
one appropriate for assessment of off-design char- 
acteristics of optimal arrangements. 

The problem treated is posed as follows: for a wing 
of given geometry, and for given gap, stagger, free- 
stream Mach number, and lift, (1) find the drag of the 
arrangement when the closed body is given, and (2) 
find the minimum drag when the body geometry is 
modified to its optimal shape for unchanged length and 
volume, and when the wing incidence is optimum. 

First, convenient expressions for the lift and drag 
of the configurations are derived. The role of an ex- 
tended slender-body theory in these results is discussed. 
After this analytic framework has been established, 
a physical and interpretive discussion of the interfer- 
ence effects is given. The variational problem is then 
formulated and solved with the Rayleigh-Ritz method. 
Numerical examples are given and results are com- 
pared to the case of no interference; these are intended 
to show the influence of the considerable number of 
geometric parameters on the phenomena. In con- 
clusion, comparison is made with the conventional 
wing-body combination when the effects of skin friction 
are included. 


(II) Expressions for Drag and Lift 


We consider a steady, supersonic flow of moderate 
Mach number M over a wing-body configuration with 
inclinations to the incident stream such that lineariza- 
tion of the equations of motion is permitted. We 
introduce a Cartesian coordinate system with the x-axis 
in the direction of the incident stream U = Ui, where 
U > 0 in forward flow and U < 0 in reverse flow, and 
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the y-axis along the span of the wing. Then the per 
turbation velocity V = ui + vj + wk is governed by 


VeWe=f (1) 
and VYxVeao (2) 
together with the definitions 
—B? 0 O 
W = ®-V, @ = 0 1 @ (3) 
0 @ 3 


Where B? = /*? — 1. Exterior to source and vor- 
ticity distributions, f and Q vanish. The vorticity 


can be decomposed as 
Q= 2i+ © (4) 


where w, being perpendicular to the x-axis, is referred 
to as the bound vorticity. If n is the unit outward 
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normal to the surfaces A of the bodies concerned, then 
the tangency boundary condition in terms of V is 


U-n = —V-n (5) 


It is understood that Eq. (5) is conveniently satisfied 
on the reference surface of the wing, whereas it must 
be satisfied on the actual surface of the body. If N is 
the unit outward normal to the section of A made by a 
plane x = const., then V-n in Eq. (5) can be written 
interchangeably with any of the terms V-N, W-n, 
W-N, essentially because the x-component of n is small. 
The surface pressure p in excess of that of the undis- 
turbed stream Po, is given by 


(p — po)/po = —U-V — 1/2V-W (6) 


When the pressure is evaluated on the wing, for which 
the approximation is fg < 1, t2 being the wing thickness 
ratio, the last term in Eq. (6) is negligibly small. 
When the pressure is evaluated on the fusiform body, 
for which the approximation is ¢,? In t; < 1, 4 being 
the body thickness ratio, only the product of the x- 
components of V and W in the last term of Eq. (6) is 
negligible. Here we are considering a configuration 
consisting of two basically different bodies. We shall 
assume that 4,2 In ¢; = O(t) <1: This establishes the 
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order of magnitude of each component of the perturba- 
tion velocity that arises in the configuration. 

Convenient expressions for the drag and lift forces 
acting on the wing-body configuration have been de- 
rived in reference 5. An integral relation, the general- 
ized reverse-flow relation, involving the forward and re- 
verse flow over related configurations of wings and 
fusiform bodies is obtained therein. This reverse- 
flow relation is subsequently applied to the forward 
flow over the present configuration and the reverse 
flow over a related configuration chosen so as to give 
convenient expressions for the forces acting in the 
forward flow. Further applications of reverse-flow 
methods to interference problems are given in reference 
8 for wing-body combinations with quasi-cylindrical 
surfaces. 

Thus the perturbation velocity induced by the con- 
figuration considered in forward flow, V;, is the sum 
of the velocities induced by the isolated body and 
wing, Viz and Vo,, and the interference velocities 
associated with each member, V’;, and V’s,, which 
are needed to preserve the boundary shapes. Hence 


Ve => Vip + Vor ++ Vip + V' op (7) 


We wish to compute the forces in forward flow but to 
avoid or reduce the labor of determining the inter- 
ference flows. With this objective, we choose for the 
reverse flow the same arrangement of wing and body 
with the same isolated geometry of each as in forward 
flow; we do not, however, introduce interference flow 
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fields, so that mutual distortions of the boundaries are 
permitted. Therefore 


Ve = Vir + Ver (8) 
and where Vir and Voz satisfy the boundary conditions 
Urn = —Vig-n on A; 

and Ur-n = —Vor-n on A» 


Here n is the surface normal in forward flow and Ur = 
—U,. By use of Eqs. (7) and (8) in the reverse-flow 
relation, the drag of the configuration in forward flow 
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Fic. 6. Drag polar for configuration Rl — A,h/l; A = 2.85, 
h/l = 0.177. 


appears expressed as a sum of the following surface 
integrals on the body and the wing: 


Dy = ri ' pen -idS + g Prn: idS — 
Ai Az 
l l 
g ( = poU r * Vy) z Vor : NdS eed g Pr - Vir . NdS 
A\ l R Az l R 


Evidently, the forward flow (7) is not necessarily re- 
quired over the entire surface of each body. Further- 
more, by re-expressing in the above the disturbance 
induced on the wing due to the body in reverse flow 
as the disturbance induced on the body due to the wing 
in forward flow, the surface integrals on the wing can be 
written® as a volume integral over the body 7}. Then 
the drag expression involves only integrations over the 
body plus a term expressing the drag of isolated wing 
in reverse flow: 


D, = g _ pan-idS + 


psf (i X oir + firi)- (Ver + V'er) dT + 
Ti 


l 
pen-idS = g, (— poU r- Vy) U Vor: NdS (9) 
433 R 


The first term of the right-hand member of Eq. (9) 
is the drag of the isolated wing in reverse flow. The 
volume integral is the drag due to cross flow and pres- 
sure of the wing, Vor + V's, induced on the singulari- 
ties representing the body in reverse flow. We antici- 
pate consideration of bodies that can be represented 
by a lineal distribution of singularities, for which this 
volume integral will, in the limit, become a line integral. 
The third term is the drag of the body in reverse flow 
due to the complete pressure field, which is quadratic 
in the complete perturbation velocity Vir + Ver. The 
last term is the drag arising from the linear part of the 
pressure field in forward flow supported on a fictitious 
slope due to the cross-flow of the wing in reverse flow. 
The last two integrals extend over the actual surface. 

Finally, it is desirable to transform the last two sur- 
face integrals on the body in Eq. (9) into volume in- 
tegrals over the body. Then, in the limit, a line 
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an 
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integral results. Furthermore, consideration of the full 
quadratic pressure and the actual surface of the body 
is avoided. As a consequence of this transformation, 
the interference effects are clearly exhibited. The 


surface-volume identity 
$ (VeWer-n + VeWer-n = Ve: Wen) dS = 
4; 


| [VeV-Wr + VeV-We — We X (V X Ve) — 
Wr X (V X Ver) |JdT (10) 


given by Ursell and Ward® and subsequently used by 
Ward® and Clarke® will be used to effect the desired 
transformation. 

We begin by expanding pz on A; according to Eq. (6), 
using the flow field given by Eq. (8), then the last two 


integrals of (9) become 


g pirn-idS — po g Vir: Ween-idS + 
Ai Ai 
¢ porn: idS — 
A, 


l 
ri (— poU,- Ve) U Vor: NdS (11) 
Ai R 


The last two terms of (11) can be written in the form 
of the first two terms of the surface integral of Eq. (10) 
by applying the boundary condition on A and noting 
the relation Vor-N = Wor-n. We obtain 


oer 5 g (U--VeWor'n —_ Ur: VerW,e-n)dS (12) 
F Ay 
which is transformed according to Eq. (10) as 
Po r . . 
= U { [Up-Vefor — Ur-Verfe — Ur-Wer X 
pT; 
wer + Ur- Wor X wr] dT — 


Po g Vie: Wori-ndS ( 13) 
A, 


in which Eqs. (1), (2), and (4) have been used. When 
(13) is simplified, by noting that f. and w. vanish in 7}, 


and substituted into (11), we obtain 
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¢ Pipn-idS —_ pf (i x @Wir Tt fipi)- VordT — 
Ai 7; 
Po ri (V; a Vir): Woern-idS (14) 
A, 


When the term (Vy + Vir)-Wee is evaluated by the 
respective components of the vectors in a local orthog- 
onal coordinate system as in reference 5, it can be 
shown that the last surface integral in (14) is negli- 
gibly small. The drag formula finally becomes 





= (0.155 


Fic. 8. Drag polar for configuration R2 — h//; h/l 


Dy = ri Pipn-idS + ¢ Popn-idS + 
A 1 


m | ) (i x QR i ifir)- (Vox + V'sp) a 
Ti 
[i X (oir + oir) + i(fir + fie) |-VerjdT (15) 


Knowledge only of the singularities representing the 
body and of the flow fields of the wing both in forward 
and reserve flow is explicitly required for the evaluation 
of Dy. Implicitly, however, the flow field of the body 
in forward flow is needed in order to evaluate V’o,. 
Let the first two terms of Eq. (15), which are the self- 
induced drags of the isolated body and isolated wing 
respectively in reverse flow, be denoted by Da,y 
and Dr». It has been shown that the drag of the 
isolated body® or wing’ in either direction is the same, 
i.e., Deyn = Dey and Dry» = Dy,22, aside from possible 
edge-force contributions in the latter case.!” 

By a similar application of the reverse-flow relation, 
the interference lift, L’2,, produced on the wing is de- 
termined as° 
Lop = ~~“ f li x (@ir + wip) 2 

me V* oR 
(fie + fir) |- 


a 


dT (16) 


where V*., is the perturbation velocity induced by a 
flat plate in reverse flow which is placed at an angle of 
incidence of a*, and which has the same planform as 
that of the given wing. As previously noted, the lift 
on the closed body is zero. This fact: can be demon- 
strated as follows. If we apply the momentum theo- 
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rem to a surface enclosing the configuration, the lift 


force is® 
L= pux ¢ n X VdS (17) 
A 


Now we shrink A to include the body surface A,, 
the wing surface A», and a tubular surface connecting 
the two. The integral over the tubular surface 
vanishes in the limit. Since the velocity V can be de- 
rived from a potential ¢, the surface integrals on A; and 
A, can be transformed into line integrals along the trail- 
ing edges where ¢ may be discontinuous though finite. 


Then 
L = poU x | f ods +f ois | (18) 


where ds is the elemental tangent vector along the path. 
We see then that the integral around ¢, which is pro- 
portional to the lift on the body, vanishes since the 
body is closed. 

We shall assume throughout that the gap between 
the wing and the body is such that on the body the 
domain of influence and the domain of dependence 
of the wing in forward flow do not overlap (Figs. la and 
lb depict two such situations) ; then the terms contain- 
ing w’iy and f’;, vanish within the region of integra- 
tion in both the drag and the interference lift expres- 
sions and so need never be determined. If the gap is 
such that the domains of influence and of dependence 
of the body do not overlap on the wing (Fig. 1b depicts 
the situation hereafter referred to as the single- 
reflection case), the drag formula is further simpli- 
fied by the vanishing of V’2, in the region of integra- 
tion. In this case, then, the drag of the configuration 
is evaluated solely from a knowledge of the flow fields 
of the isolated wing and the singularities representing 
the isolated body both in forward and reverse flow. 
If the gap is smaller than that giving the single- 
reflection case (Fig. la depicts such a situation, here- 
after referred to as the two-reflection case), V’2, must 
be determined on the rear portion of the body axis. 
When wing edges are suitable, V’2, is conveniently 
evaluated by the use of an image. More complicated 


wave reflection patterns will not be treated, since these 
seem to be adequate to this study for the Mach num- 
bers and geometries contemplated. 

Further, we shall take the bound vorticity in the 
direction of the y-axis because of vertical symmetry, 
and assume that the spatial distributions of sources and 
vorticity are zero everywhere except along the x-axis 

i.e., we take only the first two terms of the multipole 
distribution which replaces the volumetric distribution 
when the limiting process is used. Let the lineal in- 


tensity of the sources and vorticity, F and jl, be de- 
fined by 


fipdT = Fi pdx (19) 
jo, pdT => jlir dx (20) 


0.08 





004 


1, 


© 0005 0010 0015 0020 0025 0030 0035 
0 10 20 30 40 RED Cp 





Fic. 10. Drag polar for configuration R38 — h/l; h/l = 0.177 


on the x-axis. If the length of the body is /, then Egs. 
(15) and (16) become 


Dr = Dyyul Fie, lj) i D p,22( a2) = 
1 
po f ) (kT r(x) + iFir(«)]-[Ver(x) + 
V'er(x)] — [kT ir(x) + iFir(x)]-Ver(x) fdx (21) 


V*.R (x) 
, * 


Qa 


l 
L'sp = | [kT ip(x) + iFip(x)] dx (22) 


(III) The Role of Slender-Body Theory 


Since the singularities associated with only the iso- 
lated body appear in the drag expression, we now dis- 
cuss briefly the flow over such a body. The singulari- 
ties representing a body in forward flow are related 
to the geometry by the tangency boundary condition 
(5). When the same body appears in reverse flow, the 
associated singularities are implicitly related to those 
in forward flow by the expression 


(Vir + Vir) -n = (0) (23) 


Consider in particular a lineal distribution of sources 
and horseshoe vorticity representing a slender body of 
circular cross-section S;(x) having local angle of attack 
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or camber a(x). If F(x) together with its derivatives 
are of uniform order, and similarly for T\(«), we find 


F(x) = U(d dx)S (x) [1 + O(t;" In ty) | (24) 
T(x) = U(d/dx)[Si(x)ay(x)|[1 + O(4? In 4)] (25) 


when LU’ § 0. The analysis is similar to that for a 
slender body of revolution at incidence and need not 
be given here. It is seen that the singularities repre- 
senting a slender body in forward and reverse flow differ 
in sign only, when this result is substituted into Eq. 
(21), the expression for the interference drag is further 
simplified. 

However, according to this slende1-body approxi- 
mation, the self-induced drag of the vorticity distribu- 
tion is of order f,°, which is negligible compared to that 
due to the source distribution of order ¢,;'. The drag 
expression then becomes linear in I'\,; therefore no 
relative minimum exists for any Ij. Uninhibited 
by a self-induced drag, the vorticity distribution will 
grow to a point where Eq. (25) fails. Further, results 
indicate that the optimum vorticity distribution is 
sometimes insufficiently smooth to justify the use of 
Eq. (25) in the first place. By an extension of the 
slender-body theory, we can admit a less smooth vor- 
ticity distribution whose self-induced drag is not of 
negligible order, and, at the same time, obtain the same 
simple expression for the interference drag that results 
from use of slender-body theory; the latter is possible 
because the vorticity term in Eq. (21) is smaller than 
the source term. In anticipation of a rapidly varying 
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ai(x), we suppose the order of magnitude of T'\(x) be- 
comes larger by order ¢,'/? for each successive deriva- 
tive—i.e., 

T(x) /Ty'(x); Ty’(x)/T1"(x); ete. = O(4/*) (26) 
where the prime denotes the derivative with respect 
to x. By calculating the x-component of the velocity 
u induced by such a distribution of vorticity and apply- 
ing the assumption « < U used to derive the basic 
equations, one finds that T';(x) is of order ¢,° In h, in- 
stead of the #,° in slender-body theory. The self- 
induced drag becomes of order 45 In? 4. By applying 
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the boundary condition in the cross-flow plane and 
noting the relation (26), one obtains an expression for 
I'\(x) like Eq. (25) but with error [1 + O(f,; In 4) |; there- 
fore 

Tir(x) = Tye(x) [1 + O(4 In 4)] (27) 
Substituting Eq. (27) into the drag expression (21), and 
noting that V> is of order ¢,* In f;, we obtain, just as for 
slender-body theory, 
Dy = DewlFir, Tir) + Dr,22(a2) — 

Po [kT (x) Tr iF p(x) |- [Voy(x) + 

7 V 
Vor(x) + V'or(x) |dx + O(t)8 In ¢,) (28) 

Since the second term in Eq. (28) is of order t;* In? 4, 
the last term is clearly negligible. Therefore, in evalu- 
ating the drag, [')(x) may be related to the geometry 
by Eq. (25) within the approximation made. On the 
other hand, to obtain a;(x«) from a computed I';(x), the 
term of order f; In ¢; in Eq. (25) is to be accounted for as 
follows: 


I ‘ : 
a(x) = 2 f r(é)dé + 
aR? Jo 


B? \ — 4 x—& l J 
read A Ir" (&) | cosh BR 2 dé + 


x— ] 
r’(0) | cost - 


9 


Ser 


BR _ ] + ows imty (29) 


2Qr 


where & is the radius of the body. 


(1V) Examination of the Lift and Drag 
Expressions 


The nature of the interference is more easily appre- 
ciated by examining Eqs. (22) and (28) than by con- 
sidering the forward flow directly. A positive inter- 
ference lift generally results if the configuration is such 
that within the region of influence of the wing in re- 
verse flow, the source strength is positive and the vor- 
ticity strength is negative. If the wing has neither 
thickness, camber, nor incidence, then lift can be de- 
veloped on the wing without incurring an attendant 
drag. The interference drag in Eq. (28) may be de- 
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composed into (1) that due to wing angle of attack, 
Vege + Voor, (2) that due to wing thickness and 
warp, Var + Vor, and (3) that due to the inter- 
ference flow, V’2,. If we consider, for purposes of illus- 
tration, a configuration with body and wing planform 
symmetrical about the plane x = //2, then, in slender- 
body theory, Fix(x) is odd and T),(x) is even with 
respect to the plane of symmetry. The sum of per- 
turbation velocities toy” + Mer is even and Weyer + 
Wr is odd with respect to the same symmetry plane; 
the interference drag due to wing angle of attack there- 
fore vanishes. By a similar consideration, the inter- 
ference drag due to wing thickness and warp can be 


written as 


l 1 
—2p f Fy plto4p dx — 2m f Ti pW p dX 
1/2 1/2 


As is clear from physical considerations, a compression 
wave from the wing will produce a favorable effect, as 
willa positive a;. Since the portion of the body affected 
by the interference flow from the wing has in general 
negative surface slope, a compressive V's, will con- 
tribute to the drag reduction. A consideration of these 
various effects is helpful for recognizing favorable 
orientations of the configuration considered. 

When the body geometry is given, the aerodynamic 
forces acting on the configuration can be determined 
with comparative ease from Eqs. (22) and (28). (In 
the subsequent numerical calculations, the Sears- 
Haack body is taken as the given body in each of the 
configurations considered.) In order to exhibit the 
form of the polar drag curve for a configuration con- 
sisting of a given wing and a given body, which does 
not pitch with the wing, we write these drag and lift 
expressions as 


Dy = Dp, ur + D, + Dya2 + D,a2? + D, + D, ay (30) 
Le = Maz, — Qz,) sae Ager f aa 
where 0,, D,, D,, IN, and the angle of zero lift a,, 
characterize the isolated wing for given Mach number, 
and ), and , characterize the interference drag. If 
we ignore any leading-edge singularity of the wing, 
then we can rewrite the above expressions in coefficient 


form as 
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Fic. 13. Drag polar for configuration T = c//]; c/l = 0.4. 
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l Ca + Ca? 
Ca. = | mas a by + =| 
1 


n 2 
(Cp, + Co,)?} (31) 


tm 


Co, + Co, + Cry — 


Cr = mMaz+ C",, — Me, 


where m is the lift-curve slope of the isolated wing. 
It is seen that the effect of the body and the interfer- 
ence is like the effect of thickness and camber on flat- 
plate wing characteristics, with Cp,, and Cp, formally 
corresponding to the thickness and the warp part of 
Cp,, respectively, Cp, corresponding to Cp,, and C’, 
corresponding to —maz,; but there are possible dif- 
ferences in signs and magnitudes. Evidently, the 
effect of favorable interference is to shift the drag polar 
of the interfering system upward and toward the 
left. In the case of configurations without wing warp 
and with the fore-and-aft symmetry noted, the mini- 
mum Cy occurs when a, = (). The performance of 
such a configuration depends then on C’;, and Cp,. 
When the latter vanishes, the total effect of the inter- 
ference is to shift the polar curve upward. When the 
body is contoured optimally, the further reductions in 
drag at given C, are due partly to a slight increase 
in the focal length, and partly to an additional shifting 


of the drag parabola. 
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Fic. 14. Shape of optimum body; (a) configuration Rl — 7}, 
(b) configuration R2 — h/l, (c) configuration T — 7}. 


(V) Variational Problem and Method of Solution 


The variational problem requires determination of 
functions F(x) and T\,(x) which satisfy respectively 
the end conditions 

F,r(0) = Fip(l) = 0, Tir(0) = Tir() = O (32) 
and of a, such that the drag, Eq. (28), be a minimum 
subject to the following subsidiary conditions respec- 
tively on (1) the lift of the configuration, (2) the vol- 
ume, (3) the closure, and (4) the absence of total lift 


on the body: 


L = Maz — agi) — mf [kT 1 p(x) a 
0 
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V*ip 
iFip(x)]-—— dx (33) 
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. : us ] ”y 
Kw f ‘i {  Fyp(&) dt (34) 
0 0 Us 


Qe l 
F(x) dx = 0, and f Tir(x) dx = 0 (35), (36) 
0 


ed { 


The drag of the isolated body is" 


Dru = — : tf F'ip(x)F'ir(é) X 


Qn 


In | x _ t| dé dx + 2B? 


Lei 
ff l’'ir(x)T'ir(é) In |x — &| dé dt (37) 
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Fic. 15. Drag polar for configuration Rl — A (1) 


The drag of the isolated wing can be written generally 
asin Eq. (30). 

The direct method of Rayleigh-Ritz is used for the 
solution of this variational problem. For discussions 
of the convergence and limitations of this method, see 
reference 12. Then, as suggested by the technique 
employed in the Prandtl wing theory and later used 
by Sears,'* we first transform the portion of the x- 
axis from x = 0 to x = / into an angular variable 6, 
where 


x = (1/2)(1 + cos @) (38) 


We then choose the functions sin 6, which satisfy 
the end conditions trivially, to be the coordinate 
functions. Expanding Fip[x(6)] and T,[x(6)], we 
obtain 


F,p(0) = Urrl >> a, sin nd (39) 
n=1 

and Tir(@) = Urrl DY bd, sin n0 (40) 
n=1 


The subsidiary conditions (35), (36), and (34) become 
a, = 0, bd; = O, and 


a, = —(16/2*)(7T;/1') (41) 


respectively. Substituting Eqs. (39) and (40) to- 
gether with the transformation (38) into the drag ex- 
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Fic. 16. Drag polar for configuration R1 — A,h// (1). 


pression (28) and integrating, we obtain the drag writ- 
ten in coefficient form as 


. 5 

ml? — ~ 

b= — L n(a,2 + 2B%,2) + Coa + 
4S, n=2 

9 N 

. i rl? . 
Cp,a2 + Cp,a2? — 5 x 

Sy n=2 


An[Ln + al, + Tr(a, b)] + 
Bala + arty + Jy (a, b)]} (42 


in which J, and J,“ are defined as 


9 n 
[© = 7, | (tt ¢ + Uxr) sin @sin n6dd (43) 
Fv 


bo 


I, = TT, f (oar + Urer) sin 6 sin n6d@ (44) 
F 0 

and J,, which depends on all the a’s and 0’s, is de- 

fined as 


9 ™ 
I, (a; 6) = = i) u'op sin @sin n0d@ (45) 
Ur Jo 
i= 2,3,..., N, and J,, J,™, and J, are similarly 
defined for wer + wee and wo’. The only remaining 
subsidiary condition, the prescribed lift, becomes 


rl? x 
Cr = m(a2 — agi) — © DL (aaln* + brJn*) (46) 
2 Uor™ ‘ a 
where /,* = — | ; sin 6 sin 6 d9 = =(47) 

L F 0 a 


and J,,* is similarly defined for wor*/a*. Since bodies 
with discontinuous surface slope and camber are not 
considered here, convergence of the sums = na,’ and 
> nb,2 as N > © is assured. Terms like = a,J, 
are convergent because /,,“ can be regarded as the Four- 
ier coefficient of (22,7 + Uxr) sin 6 and behaves like 
1/n. The drag and lift of a configuration with a 
Sears-Haack body are simply evaluated from Eqs. 
(42) and (46) by substituting a, according to (41), and 
by setting a, = O form > 3, and b, = Oforn > 2. 
The variational problem reduces now to a problem 
in the ordinary calculus of maxima and minima, in terms 
of 2(N — 1) variables with one subsidiary condition. 
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Fic. 17. Drag polar for configuration R2 = h// with friction 
included. 


The necessary conditions for the existence of a rela- 
tive minimum give the following equations for a,, b,, 
a, and the Lagrangian multiplier ): 


(x?/2)nay + 1, + al, + 
(0 04») [AnD (ai, b,;) + DJ,“ (dj, b;)| —_ 
7,* = 0 (48) 


1 Bnb,, + DP + a2J ,\) ae 
(0/Obn) [AnD n (aj, b:) + b,J, (aj, b:)| — 
\J,* = 0 (49) 


9 N 


De (aula + dyn) + 


Am = 0 (50) 


Cp» 2 2a2C p, + 


where n > 3 in Eq. (48) and m > 2 in Eq. (49). These 
2(N — 1) equations together with the lift expression 
(46) will suffice to determine the optimum values of 
fie Ga HOE 1 = SF, <s.0-> LW) Oe COE = 2.3, .«.. 54), 
and X. It is noted here that the optimum geometry 
of the body varies with the design lift. The minimum 
drag is obtained when these values are substituted 
back into Eq. (42). The sufficient condition for the 
existence of a relative minimum can be obtained by 
considering the second-order terms in the Taylor’s 
series expansion of the drag expression about the opti- 
mum point. Having obtained the quadratic form, we 


find that if 


2Cp,. al 

9 N 9 

rl? CE 1))2 

: ; >)» > eee an eres 

Se n=3 (n1r", 2) + (0? Od ,,”) ia,1,° + b,J,] 
Sy 
n=2 nw? + (02/da,”) [anln + d,Jn]$ 

in which the derivatives are evaluated at the stationary 

point, then the stationary point is a minimum. 


(VI) Numerical Results and Discussion 


In this section, the results of numerical calculations 
for four basic configurations, three with a rectangular 
wing and one with a triangular wing having supersonic 
leading edges, are presented in order to show (1) the 
magnitude of the drag reduction made possible by the 
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interference considered, (2) the further improvement 
resulting from the substitution of the optimum body 
for the Sears-Haack body, (3) the importance of the 
gap between the wing and the body, and (4) the effect 
of the Mach number and various geometrical param- 
eters on the drag reduction. Both types of wings have 
symmetric diamond profiles of 5% thickness ratio 
with the maximum thickness at the mid-chord. The 
sweep angle of the triangular wing is taken to be tan™! 
0.86. From a qualitative consideration of the drag 
formula, we locate the wing centrally above the body 
in all configurations considered; no studies on the 
effect of stagger are made. The gap h, except in the 
low-gap studies, is taken to be //25, so that the inter- 
ference flow V’s, is not needed for calculating the drag. 
The lower gaps are so chosen for each wing planform 
that V’,, can be determined conveniently by the 
method of images. The aspect ratio A of the wings is 
chosen such that the influence of the body is not felt on 
the upper side of the wing; with this choice, the flow 
field of the rectangular wings at the body axis is con- 
veniently two-dimensional. For all basic configura- 
tions considered, the design Mach number is taken to 
be 2.24 (B = 2) and the given volume of the body is 
taken to be 0.00656/*, corresponding to a Sears-Haack 
body (referred to hereafter as the S-H body) of 12% 
thickness ratio. The values of geometrical parameters 
defining each basic configuration are listed in Table 








TABLE 1 
Wing 
Configuration planform A c/l 
Rl Rectangular 1.40 0.293 
R2 7 1.62 0.380 
R3 a 1.38 0.467 
é Triangular 2.50 0.500 





Modifications in the basic configurations are introduced 
in order to study the effect of varying a parameter; 
these modified configurations are designated by R1-t. 
if the wing thickness alone is changed, and so forth. 
The values of the modified parameters are indicated 
on each graph presenting the results. Extensive 
studies of the effects of different geometric parameters 
have been carried out for configuration R1; the inter- 
relationship among these studies is exhibited in Sketch 
1. The symbol (I) denotes the fact that the calcula- 
tions are made by the indirect method; the results 
are discussed in a subsequent section. 








l 
fk1-A(I)}4r1-A, h/2 (1)} 
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ent 
ody 
the 
‘ect 
4m- 
ave 
itio 
The 
a 
rag 
dy 
the 
the 


ag. 
rm 
the 
3 1s 
on 
OW 
on- 
ra- 
to 
1s 
ick 
y// 


ers 








BODY UNDER 


In each configuration, the procedure of calculation 
is as follows: We first calculate the polar drag Cp,,, of 
the isolated system consisting of the given wing and 
the S-H body in isolated flow, then the wing and the 
body are brought together to form an interfering con- 
figuration and its polar drag Cp, is evaluated. 
Lastly, we calculate the minimum drag Cp,,, of the 
configuration by shaping the body optimally for each 
given. lift. The results, reduced to coefficient form 
with respect to wing area, are presented in the figures 
in the form of these three drag polars. On each of 
these curves, the point of maximum ratio of lift to drag 
is indicated by a circle, the vertices of the Cp, an 
Cp,» Curves are indicated by a horizontal dash. In 
addition, three curves of percentage drag reduction are 
presented in the same graph. The circle on the 
abscissa denotes the drag of the isolated S—H body. 
Two views of the configuration are sketched on the 
graphs for each of the basic configurations considered, 


and for those modified configurations in which the 
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Fic. 18. Schlieren photographs of flow about wing-body configuration (knife-edge horizonta!). 
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geometry of the wing or the gap is changed. The 
equations of the drag polars for each configuration are 
listed in reference 14. 

The integrals involved in the calculations are evalu- 
ated in reference 14. In the single-reflection cases, 
the integrals are elementary for the rectangular wing 
and elliptic for the triangular wing. For the two- 
reflection cases, the value of V’., at the body axis is 
approximated by the asymptotic flow field of the image 
body. (When V‘%, cannot be evaluated by the 
method of images, it might sometimes be possible to 
approximate V’,, by a flat-plate flow with the average 
w,/lU,» determining its angle of attack.) The result- 
ing integrals again are elliptic and are evaluated accord- 
ing to reference 15. All calculations are made by tak- 
ing four terms—1i.e., V = 5—in the series expansion of 
F\,-(x) and T;,(x). The rapid convergence of the series 
terms in the drag expression for the single-reflection 
cases may be indicated by the value of Cp,,,, obtained 


in configuration R1-A: 
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Crypt = 0.0107 + 0.539(C_ — 0.0178)? for N = 20 
Cp, p, = 0.0107 + 0.540(C_ — 0.0176)? for N= 5 


No such comparisons have been made in the case of 
triangular wings or in any two-reflection cases. A 
brief discussion of the results obtained for each of the 
configurations considered is presented below. 

Configuration R1: All the features of favorable inter- 
ference discussed previously are clearly exhibited here 
(see Fig. 2). The drag polar shifts when the inter- 
fering arrangement is assumed by the given wing and 
body. When the body is contoured optimally, the 
drag polar shifts again and opens up. The relative 
drag reduction obtainable is on the order of 25%, how- 
ever, the advantage gained by contouring the body is 
small. 

Configuration R1-1\: The effect of changing the 
body volume to that corresponding to an S—H body 
of 10.4% thickness ratio on the drag reduction is not 
pronounced (see Fig. 3). Compared with the previous 
configuration, the maximum drag reductions from C’p,,, 
vecur at a lower lift coefficient. A useful rule of thumb 
may be deduced here: the maximum relative drag 
reductions occur at a lift for which the drag of the 
isolated wing and body are nearly the same. The 
shape of the optimum body corresponding to C;, = 0.15 
has been calculated and is presented in Fig. 14a. The 
body camber is calculated by the slender-body theory 
approximation, which is justified a posteriori. The 
body has fore-and-aft symmetry; its cross section is 
bigger at the two ends and smaller at the middle than 
that of the S-H body, and its centerline curves up 
like a circular arc with maximum camber of 2%. The 
large rate of increase of the cross section and the upward 
camber in the front result in a higher pressure acting 
on the rear portion of the wing, while the symmetrical 
geometry in the rear portion of the body provides a 
steeper slope for the wing pressure to act on. Of 
course, the distortions are limited by the self-induced 
drag. 

Configuration R1-t.: The effect of the wing thick- 
ness is studied here by setting the thickness ratio of 
the wing equal to zero. Fig. 4 shows that the inter- 
ference drag vanishes by virtue of symmetry and zero 
thickness of wing. The reduction in drag is due solely 
to the interference lift, as is manifested by the absence 
of lateral shifting of the Cp,,, and Cys, curves. An 
interesting example is offered here. Consider the 
vertex of the Cp, curve: it represents, an S-H 
body configuration with a level flat plate placed cen- 
trally above. This configuration has the drag of an 
S-H body, which is minimal for given volume and 
length, while at the same time a lift of Cr, = 0.033 is 
generated on the flat plate by the flow field of the body. 

Configuration R1-A: When the aspect ratio of the 
wing is increased, the effectiveness of the favorable 
interference due to the body is reduced (Fig. 5). The 
reasons are twofold: first, the region of interaction be- 
comes relatively smaller as the aspect ratio is in- 
creased, and second, the wing is so efficient in carrying 
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the lift that the optimum configuration does not re- 
quire body distortions for good lift interference. The 
results shown in Fig. 5 are to be expected. All three 
relative drag reductions are smaller than in the pre- 
ceeding configurations. The large aspect ratio is 
needed also for the off-design case, the low-gap case, 
and the case in which the inverse method is used. 
For all the configurations considered thus far, the most 
favorable drag reduction from that of the isolated 
system is on the order of one-half the drag of the isolated 
S-H body. 

Configuration R1-A,\M: To study the off-design 
conditions, the drag of the optimum configuration cor- 
responding to C, = 0.14 in R1-A is calculated at the 
reduced Mach number of 1.73 (Fig. 5). The polar 
drag of the S-H configuration is also calculated for 
that speed. For a meaningful comparison with the 
results obtained at design speed, the off-design force 
coefficients are normalized with respect to the original 
dynamic pressure and are presented also in Fig. 5, with 
the diamond point representing the off-design optimum 
configuration. Off design, this configuration evidently 
has a higher drag than that of the S-H configuration 
at the same lift because the optimum body is distorted 
in a special manner to give the minimum drag only for 
a particular design condition. When the speed is in- 
creased from that of the design condition, it is con- 
jectured that the favorable interference decreases 
since the regions of interaction both on the wing and 
on the body become smaller. 

Configuration R1-A,h/l: When the gap is reduced, 
the interference becomes much more favorable (see 
Fig. 6). All three relative drag reductions are larger 
than those in configuration R1-A except at very low 
lift; the maximum drag reduction for Cp,,, has in- 
creased from 20% in R1-A to 32%. The importance 
of body distortions is shown by the large relative drag 
reduction of Cp,,, from Cps,, with an average value 
of 7% and a maximum of 9%. Two factors are re- 
sponsible for the more favorable effects: (1) the larger 
interference lift due to the larger region of influence of 
the body on the wing, and (2) the favorable interference 
drag due to V’2y. In this case, interference drag due to 
the wing thickness is not important. 

Configuration R2: This case, given in Fig. 7, is 
studied primarily to provide a basis of comparison for 
the next study, in which the gap is reduced. The 
choice of parameters is such as to give the minimum 
gap and span under the limitations we have specified. 
The main features of these results resemble those of 
configuration R1. 

Configuration R2-h/l: In this case, shown in Fig. 8, 
the gap is reduced to 0.155/, or 1.3 times the maximum 
diameter of the body. The maximum drag reduction 
for the optimum configuration is increased from 26% 
to 44%. The favorable effects obtained by body dis 
tortions are even more pronounced than those of con- 
figuration R1-A,///; the relative drag reduction from 
Cps-y 1S on the order of 20% with a maximum of 
26%. By comparing the Cp,., curve with that of the 
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Fic. 19. Experimental lift on wing. 


configuration R2 (Fig. 7), it is seen that below C, = 
0.08, the drag of the low-gap configuration is greater 
than that of the configuration with large gap. This 
illustrates the fact that reflecting the body flow field 
onto the body again is not necessarily advantageous if 
the body is not optimally contoured. The calculated 
shape of the optimum body corresponding to C, = 
0.12 is presented in Fig. 14b. The body is no longer 
symmetric about the plane x = //2. The slender- 
body approximation is used to obtain the shape of the 
centerline, which is wavy and bends upward in the rear 
portion of the body. This shape may be affected by 
the asymptotic approximations for V’s,. 

Configuration R3: The results of this study (see 
Fig. 9), together with those obtained for R1 and R2, 
indicate that in configurations of a body with a rec- 
tangular wing, when the gap is restricted to give a single 
reflection, the drag reduction from that of the isolated 
system is insensitive to variations of the wing chord- 
body length ratio ¢//. 

Configuration R3-h/l: Since the gap in this con- 
figuration (Fig. 10) is the same as that in configura 
ation R1-A,h//, the results obtained in the two con- 
figurations do not differ greatly even though the plan- 
forms are different. From the three low-gap studies, 
it may be concluded that the gap distance is of great 
importance in obtaining favorable interference, and the 
most favorable drag reduction is on the order of the 
drag of the isolated S-H body. 

Configuration T: The asymmetry of the triangular 
wing introduces no radical changes in the results (see 
Fig. 11). The relative drag reduction is about 25%, 
and the advantages gained by body distortions are 
small. 

Configuration T-7;: It is seen in Fig. 12 that the 
effect of reducing the volume of the body is much the 
same as for the rectangular wing. The rule of thumb 
deduced from configuration R1-7; applies in con- 
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figurations with a triangular wing also. The shape 
of the optimum body at C, = 0.14 is calculated and 
presented in Fig. l4c. The slender-body approxima 
tion is used to obtain the camber. The cross section, 
compared with that of the S-H body, is greater over 
the front half and smaller over the rear half. The 
centerline has almost constant slope up to x = 0.6), 
and then it becomes horizontal. 

Configuration T-c/l: Except for the fact that the 
maximum relative drag occurs at higher lift (see Fig. 
13), which is consistent with the rule of thumb, the 
effect of reducing the wing chord-body length ratio 
is not important. We see again that drag reductions 
are insensitive to c// for single reflections. 


(VII) Conclusions Drawn from Numerical 
Results 


(A) For single-reflection cases involving both types 
of wing planforms considered, the following conclusions 
may be drawn: 

(1) With the lift coefficient for maximum drag reduc- 
tion, the drag reduction obtained through interference 
is about one half the drag of the S-H body in isolated 
flow. 

(2) The advantages gained with body distortions 
are not important. Therefore, calculations should 
be made for the S-H configuration, since these are 
evidently very simple. 

(3) When the aspect ratios of the wings are approxi- 
mately the same, the relative drag reduction from 
Cp,,, 18 not sensitive to variations in c// and the body 
volume. However, the maximum relative drag re- 
duction will occur at a lift for which the drag of the 
isolated wing is nearly the same as that of the isolated 
S-H body. 

(4) When the Mach number of the uniform stream 
is lowered from the design value, the drag of the off- 
design optimum configuration is higher than that of the 
S-H configuration at the same lift. 

(B) For the two-reflection cases considered, the fol- 
lowing conclusions may be drawn: 

(1) The drag reductions obtained were considerably 
higher than with the larger gaps. At the lift coefficient 
for maximum drag reduction, the drag of the optimum 
configuration is less than that of the isolated system 
by approximately the drag of the S-H body in isolated 
flow. Fig. 8 shows a maximum drag reduction of 44%. 

(2) The effect of body distortions is very important 
here, although the distortions required are quite moder- 
ate. Thus at low C, the drag of the S-H configura- 
tion may be larger for the lower gap than for the higher 
gap. 

(3) The principal contributions to the favorable 
interference effects are the high interference lift gener- 
ated on the wing, and the action of the interference 
flow of the wing V’2,» on the body. 

(4) The low-gap case should be studied also for 
sweptback wings. The improvement obtained by in- 
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creasing the region of interaction on the wing suggests 
the possibility of sweeping the wing in such a way that 
the region of interaction includes the entire planform. 
The influence of stagger should also be studied. 


(VIII) Interference With Mutual Distortion 


In this section we shall consider the first formulation 
of the problem, as discussed in the introduction, in order 
to conipare results. The drag of the two singularity 
distributions in forward flow is conveniently written 
according to Eq. (23) of reference 5. In this case 
the same singularities (but with different signs) appear 
in both flow directions. Eq. (28) results again except 
that Vo,’ = 0, and Voz is reinterpreted as above. No 
restrictions on orientation or smoothness now appear. 
The lift of the configuration is given merely by the 
integral of the wing vorticity. Again we expand the 
body singularities into four-term Fourier series and seek 
the optimum coefficients subject to subsidiary condi- 
tions (84), (35), and (36). For the single-reflection 
case R1-A(I), the results are nearly the same as those 
obtained with the direct formulation in Fig. 16. For 
the two-reflection case R1-A,///(1), the two principal 
factors responsible for good drag reduction in the direct 
formulation, namely, the good lift interference and the 
interference flow field of the wing, are absent here. 
The results for this case, shown in Fig. 17, are seen 
to be quite different from those obtained for configura- 
tion R1-A,h/l (Fig. 6). The drag reduction in this 
case, which is due principally to the upwash field of the 
body acting on the vorticity distribution of the wing, 
is not as great as that obtained when the direct formu- 
lation is used. 


(IX) Further Considerations 


Although the isolated system is the logical basis with 
which to compare configurations with S-H and optimum 
bodies, the question of practical interest is: how do 
the interfering configurations compare with the con- 
ventional wing-body combination? Furthermore, how 
does the interfering configuration perform when friction 
is accounted for? We shall now give some estimates 
for these two effects. 

Without going into the details of the interference 
encountered in a wing-body combination, we make the 
simple assumption that the portion of the wing blank- 
eted by the body does not incur thickness drag but 
supports the lift distribution which would normally be 
found there. To estimate the skin friction, we assume 
that separation of the boundary layer owing to inter- 
action with the shock systems does not occur. De- 
tailed calculations have been made for the configura- 
tion R2-h/l. The Reynolds number was assumed to 
be 2.65 X 10® per ft, which corresponds to M = 2.24 
at 50,000 ft, and the length of the body was taken to 
be 120 ft. Data for a compressible fully turbulent 
boundary layer along a flat plate’® ' were used. It 
was found that the coefficient of drag due to friction is 
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0.0044 for the wing-body configuration and 0.0039 for 
the conventional wing-body combination. The results 
in the form of drag polars are shown in Fig. 17, where 
the drag of the wing-body combination Cp,,,,,, instead 
of Cp,,,, is the basis for comparison. The maximum 
drag reduction decreases from 44% to 33% when fric- 
tion is included. Fig. 17 also shows that at low C, 
the S-H configuration has more drag than the wing- 
body combination. 

In the absence of friction, the minimum wave drag 
of the body of given volume and length (S-H body) 
decreases indefinitely as the given length is increased 
indefinitely at constant volume. When friction is 
taken into account, however, the elongation is ac-om- 
panied by an increase in friction drag, and there is an 
optimum fineness ratio, on the order of 20. The same 
is evidently the case when the body is part of the iso- 
lated system and the lift is also given. However, 
when interference is present, the increase in body 
length is accompanied by a change in the interference 
drag which might counterbalance the improvement 
in the isolated (inviscid) drag of the body. It is inter- 
esting to determine whether there would exist then an 
inviscid minimum drag corresponding to an optimum 
ratio of the body length to wing chord, //c. Such a 
calculation can be made using the three low-gap studies 
discussed, and it is found that an optimum //c for 
maximnm lift to drag ratio, (L/D)maz, does in fact 
exist. The calculations made for these low-gap cases 
are modified by changing the wing aspect ratio to 1.85 
while keeping the body volume constant at 0.0644c° 
for all three bodies of different lengths. For each case 
the polar drag of the optimum configuration is calcu- 
lated, from which the maximum L/D is determined. 
The results obtained are shown in Table 2. It is seen 











TABLE 2 
Approximate 
fineness ratio 
l/c h/c (L/D) nes of body 
2.14 0.378 9.4 8 
2.63 0.407 10.6 11 
3.41 0.603 of 17 





that the optimum fineness ratio of the body in an inter- 
fering configuration is very much lower than that found 
in the isolated system or wing-body combination, and 
this is a decided advantage from the point of view of 
structural design. 


(X) Experimental Study 


The experimental investigation, performed in the 
Brown University 9 X 9-in. transonic-supersonic wind 
tunnel, has been described in more detail in reference 
18. The design of the configuration tested was based 
essentially on that of the configuration R1-A,h/l. 
The rectangular wing with symmetrical diamond pro- 
file had a 5% thickuess ratio, a 2.05-in. chord, and a 
6.82-in. span. The aspect ratio was increased to ac- 
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commodate the larger nose shock cone issuing from 
the body at this lower Mach number. The 5-H body, 
which had a 12% thickness ratio and a 7-in. length, 
was truncated at one end to provide room for a */j¢-in. 
diameter sting. The wing was located centrally above 
the body; the gap between the centerline of the body 
and the wing was set at 17/gin. The design gap, based 
on favorable interference considerations, was 1.95 in.; 
thus an error of about ().1 in. was introduced during 
installation. All tests were performed at a Mach 
number of 1.6 and a Reynolds number of 131,000 per in. 

In regard to the boundary-layer flow, no significant 
separation on the wing was noted for any angle of 
attack up to 6°, the maximum angle used. Simi- 
larly, separation was never observed on the interior 
side of the body. However, on the exterior side of the 
body the flow usually separated after interacting with 
the impinging shock from the wing. 

A sequence of schlieren photographs showing the 
phenomena is presented in Fig. 18. (The pictures 
show some extraneous shocks due to flow nonuni- 
formity, and, in some cases, spurious streamlines due 
to compressor oil on the window.) The two numbers 
on each photograph denote respectively the wing angle 
of attack and the body angle of attack. The first 
photograph (00) exhibits the separation discussed. 
The manner in which the separation occurred suggests 
that the flow was still laminar up to the zone of inter- 
action. Therefore, a boundary-layer trip which was 
made of a loop of 0.005-in. diameter wire, was placed 
2!/, in. downstream from the nose. With this trip, as 
shown in photograph (OUT) in Fig. 18, the flow seems 
to have remained attached after interacting with the 
shock. (Note that the Prandtl-Meyer fan crosses the 
region of interest.) This trip was used for the re- 
mainder of the tests. When the wing angle of attack 
was increased to 1°, the fowseparated again. As shown 
by the photographs in the remainder of the left column 
in Fig. 18, the separation became increasingly pro- 
nounced as the wing angle of attack was progressively 
increased. 

In the next sequence of photographs, the angle of 
attack of both wing and body are progressively in- 
creased from zero. This arrangement eliminates 
cross flow due to wing angle of attack but results in 
cross flow due to the free stream. The flow on the 
exterior side of the body first separates as the angles 
are increased from 0° to 1°; then as the angles are 
increased to 4° the flow seems to stay practically the 
same. 

The Reynolds numbers for these tests were much 
lower than for the flight range. At the suggestion of 
the authors, the Transport Division of the Boeing 
Airplane Company kindly ran a similar test at a higher 
Reynolds number (in the range of 9 X 10° based on 
body length) and also at the higher Mach number of 
2.5. Their results were essentially the same as those 
reported here. It appears, however, that separation 
was eliminated entirely when the wing and body were 
pitched together. The separation which eccurs when 
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the wing alone is pitched resembles the well-defined 
vortex pair which appears on an isolated body at 
sufficiently high incidence to the stream. 

It is clear from the evidence that in this case the 
boundary layer managed to penetrate the incident 
and reflected shocks on the interior side of the body 
without separating, and yet it separated on the ex- 
terior side where the incident pressure rise at the shock 
is virtually canceled locally by the reflected cross flow. 
From this evidence and from what is known about 
boundary-layer behavior on isolated bodies, one is led 
to examine the pressure field established on the body 
by the reflection of the shock-induced cross flow. Con- 
sider the equivalent problem of flow past an open- 
nosed body (with independent internal and external 
flows) set at a negative angle of attack, —a, to the 
stream; this problem has been treated by Lighthill.' 
We seek the pressure gradient encountered by each 
fluid particle upon crossing the shock, in terms of the 
substantial derivative of the pressure coefficient C>,. 
We obtain at the shock 


DC,/Di = (Ua/B?R)|—cos 6 — 2Ba sin? 6 + 

O(t;? In t)) ] 
in which 6 is measured from the top-most point of the 
cross section. The terms of order ¢,? In 4, are related 
to the pressure coefficient of a slender body at angle of 
attack. It is seen that the particles entering the shock 
on the exterior portion of the cross section encounter a 
positive gradient, while particles on the interior por- 
tion encounter a favorable gradient. The adverse 
pressure gradient plus the accumulation of material 
due to the cross flow cause the flow to separate there. 

We now see that the separation may be avoided in 
practice by kinking the body to accommodate cross 
flow due to wing incidence and thickness, and/or by 
pitching the body. We may kink the body for wing 
thickness and pitch the body with the wing; the latter 
is automatic in an actual aircraft configuration, since 
the body and the wing would be connected by struts. 
However, if the body is pitched with the wing, the 
separation evidenced on an isolated body at incidence 
to the stream may occur. Therefore it seems more 
advisable to kink the body for wing thickness and 
design angle of attack without pitch; for angles of 
attack other than the design angle, the wing cross 
flow arising from this deviation is eliminated by pitch- 
ing the body with the wing. 

Lift measurements on the wing, both in isolated 
flow and in the arrangement with the body, are pre- 
sented in Fig. 19 in coefficient form as a function of 
geometric angle of attack ageom. The theoretical lift 
coefficient for the isolated wing C72 is calculated accord- 
ing to linear theory. The additional or interference 
lift on the wing due to the body flow field, Cz’ is 
calculated according to Eq. (22) as 0.05. (Estimates 
of the nonlinear wave propagation, which effectively 
shifts the range of integration in the equation, suggest 
that the correction is on the order of 0.005. This is 
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seen to be the order of discrepancy in Fig. 19.) These 
straight lines were shifted by about 0.1° to the true 
zero angle of attack. The measurements are seen 
to be in excellent agreement with the predicted values. 
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(Continued from page 546) 





This last equation gives the ratio of heats transferred 
with and without a magnetic field at the stagnation 
point, and is quoted from references | and 2 where it 
was originally established in order to show the full sig- 
nificance of the parameter Z in both the viscid and 
inviscid regions. 
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A Summary of the Supersonic Pressure Drag 
of Bodies of Revolution 


DEANE N. MORRIS* 
Douglas Aircraft Company, Ince. 


Summary 


A number of approximate theories for supersonic and hyper- 
sonic flow over bodies of revolution at zero angle of attack are 
appraised by a critical comparison with characteristics and 
second-order results, with the use of hypersonic similarity as a 
basis for the comparison. Most of the approximate theories 
are inadequate except over very limited ranges of fineness ratio 
and Mach number. The combination of second-order super- 
sonic theory and second-order shock-expansion theory provides 
consistently good results throughout the supersonic speed range. 

On the basis of exact (or nearly exact) supersonic solutions and 
a limited amount of test data and theory in the transonic region, 
summary design curves are developed that give the pressure 
drag of conical and ogive noses and conical and ogive boattails 
over the complete range of transonic, supersonic, and hypersonic 
Mach numbers. Other shapes can be analyzed in the same man 
ner, provided that an equivalent amount of data is available. 

The analysis is made with the assumption of inviscid flow, 
so that the effects of boundary-layer growth, shock-boundary- 
layer interaction, and flow separation are not included. The pres 
ent correlations provide a sound basis of inviscid-flow results 


from which these additional viscous effects can be evaluated 


: 
Symbols 
Cp, = pressure drag coefficient based on maximum cross- 
sectional area 
Cy, = pressure coefficient = (p — pao)/ge 
d = maximum diameter of body 


FR = fineness ratio = 1/d 
Ky = B/FRy (see Eq. 5) 
Ks = see Eq. (6) 


l = length of body 

M = Mach number 

p = static pressure 

q = dynamic pressure = (y/2)pM? 
r = local radius of body 

R = maximum radius of body 

U = velocity 


@ =~WM.?-1 
= ratio of specific heats 

#y = half-angle of body nose (see Fig. 1) 
Subscripts 


© = free-stream value 
BT = boattail 

C = center section 

N = nose 


Introduction 


A FEW YEARS AGO, the problem of calculating the 
pressure distribution and pressure drag of a body 
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of revolution at zero angle of attack in supersonic flow 
was a subject of considerable discussion among aero- 
dynamicists. Many approximate methods of solution 
were proposed, but the absolute accuracy of the approxt- 
mations remained in doubt. Reliable test data were 
lacking and, although the method of characteristics 
was available, the complexity of a characteristics 
solution was a strong deterrent to its widespread use. 

Today, the fundamental problems of inviscid axially 
symmetric flow can be analyzed on a rational and 
straightforward basis. This situation is due prin- 
cipally to the developments of the second-order and 
shock-expansion theories, the use of high-speed com- 
puting machinery, and the application of the hyper- 
sonic similarity laws. It is the purpose of this paper 
to reappraise briefly the accuracy of the approximate 
supersonic theories (as well as several approximations 
for hypersonic flow) and to develop from the analysis 
a set of semiempirical curves for the pressure drag of 
bodies of revolution. 


Hypersonic Similarity 


Summaries of the hypersonic similarity law may be 
found in numerous articles in the current literature'~*; 
the fundamental derivation of the rules was originally 
presented by Tsien*® and Hayes.*’ As applied to the 
pressure drag at zero angle of attack of slender bodies 
of revolution of the same family shape, the law states 
that the drag parameter M=*Cp, is a function only 
of the hypersonic similarity parameter Me/FR, pro- 
vided that the free-stream Mach number is large com- 
pared to one. The limits of application of the hyper- 
sonic similarity law were investigated,” * and it was 
established that similarity existed for much smaller 
Mach numbers and fineness ratios than the theoretical 
derivation would suggest. Later, Van Dyke* showed 
that the hypersonic and supersonic similarity rules can 
be combined artificially, thus extending the limits of 
validity of hypersonic similarity nearly to transonic 
speeds. This unification is accomplished simply by 
replacing the Mach number Me in the similarity 
parameters by 6 = V Ma? — 1. 
parameter is multiplied first by the square of the hyper- 
sonic similarity parameter, the net result is: 


If the pressure-drag 


Cp,(FR)? = f(8/FR) (1) 
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Fic. 1. Geometry of bodies of revolution. 


Strictly speaking, Eq. (1) is in supersonic similarity 
form, but the term hypersonic similarity will be retained 
in order to emphasize the extension to hypersonic flow. 
The term 8/FR will be referred to hereafter as the hyper- 
sonic similarity parameter. 

On the basis of a comparison with the results of New- 
tonian impact theory, Van Dyke has suggested a further 
improvement in the hypersonic similarity rules,’ par- 
ticularly for thicker bodies. Hypersonic similarity 
theory indicates that the pressures (and the pressure 
drag) are proportional to the square of the tangent 
of the body-surface slope. However, impact theory 
assumes that the pressures are proportional to the 
square of the sine of the surface slope of the body. Van 
Dyke shows that much better correlation is obtained 
for cones if the hypersonic similarity rule is modified 
to a sine-squared variation. For cones of semivertex 
angle 6y, the relations between the nose angle and the 
nose fineness ratio are (see Fig. 1) 


tan? dy = 1/4(F Rw)? 
sin? @y = 1/4[(FRy)? + (1/4)] 


For cones, then, the tangent implicit in the left-hand 
side of Eq. (1) is replaced by the sine simply by re- 
writing it in the form: 


Co, |(FRy)? + (1/4)] = f(8/FRy) (2) 


In the case of ogive nose shapes, the geometry is not 
quite so clear-cut, but Eq. (2) would represent a sine- 
squared variation with regard to the nose angle of the 
inscribed cone. Further modifications will be dis- 
cussed later. 

In correlating the data, it is desirable to cover as 
wide a range of variables as possible in a single plot. 
To accomplish this without loss of accuracy, the inde- 
pendent variable in all plots is the hypersonic similarity 
parameter 8/FR for a range from zero to one; then this 
parameter is inverted, to FR/8, for a range of values 


from one to zero. In this way the graph covers the 
Mach number range from one to infinity. 


Conical Noses 


The pressure drag of the axially symmetric cone in 
supersonic flow will be considered first, since exact 
results, in the form of the Taylor-Maccoll solutions, are 
readily available. A number of approximate solutions 
for cones are shown in Fig. 2, plotted in hypersonic 
similarity form. It should be noted that changes 
in the parameters such as the replacement of tan @y 
by sin #y are indistinguishable in all of the approximate 
theories, since the body is invariably assumed to be 
“slender.” 

The limits of validity of each of the approximate 
theories have usually not been specified by the authors, 
and in general cannot be exactly defined. Therefore, 
most of the theoretical curves have been terminated 
arbitrarily at 8/FRy = 1, although this is not meant 
to imply that the theories are valid to this value of the 
similarity parameter. The later comparisons with 
exact theory will allow an estimate of the regions of 
validity of each theoretical result. 

It will be noticed also that both “‘linear’’ and ‘“‘first- 
order” results are shown for a number of the theories. 
While it is now generally agreed that the linear solu- 
tions are basically incorrect, they are quite well known 
and it seems advisable to retain them in this com- 
parison. 

The linear slender-body theory is that presented 
by von Karman.’ Written in hypersonic similarity 
form, the expression for pressure drag is 

Co,y(FRwn)* = (1/2) In (4FRy/8) (3) 
The first-order result is obtained directly from linear 
theory by retaining an additional term in the expression 
for the pressure coefficient. The relation for pressure 
drag becomes 


Copy(FRy)*? = (1/2) In (4FRy/B8) — (1/4) (4) 
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Fic. 2. Supersonic pressure-drag parameter for conical noses, 
according to various approximate theories. 
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Fic. 3. Supersonic pressure-drag parameter for conical noses. 
Comparison of exact and approximate theories. 


The second-order slender-body theory of Broderick,*! 
when reduced to hypersonic similarity form, results in 


Cp, y(FRn)? = (1/2) In (4FRy/8) — (1/4) + 
(1/16)(8/FRy)?[3 In2(4FRy/8) — 
5 1n (4FRy/8) + (113/20)] (5) 
The supersonic linear and first-order theories are given 
by the respective formulas 
[cosh—!(2FRy/8) | 
Cop FR)" = \/4 — (8/FRu)? 
[cosh—! (2FRy/8) | 
V4 — (8/FRw)? 


(6) 


Co,y(FRw)?* — (1/4) (7) 


The supersonic second-order solution is derived from 
Van Dyke’s Eqs. (13) and (14) of reference 8. Neglect- 
ing higher-order terms and assuming that \/* = £6? 
in certain places, the pressure drag of a cone can be 
written (with Ay = 6/FRy) 

10 \ Ki" 2K 77? 


< l comment 


Co. Fany = > 9 
sas ' 7Ku* | 3 4— Kr? 


os 3Ku° 
10 
sech—! Ku 2 


xX 
V4 - Ay* 


Ky?(40 + Ku?) 
16(4 — Ky?) 


sech—! Ky/2 
V 4 sad Ki 
59K 77" (7, . 10 
40(4 — Ky?) \ 7K}? 
(8) 
The hypersonic small-disturbance theory is that pre- 
sented by Van Dyke’ in the form of a numerical solu- 
tion. 
The hypersonic approximate theory was proposed 
by Lees? on the basis of an expansion of the conical- 
flow equations when the shock lies close to the body. 


The result is: 
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Cy..(FRy)? Ks* — 4 4 
Ipn\ LAN) = > A 
ee , (y + LAr? 

1 (= i) | Y¥+1 | 

YA Kn = l ae (S Ks?) 

(9) 
. Yt 1, 
Ks = —— Kn + 
v¥+3 " | 


7+ 1 


Ve43) * v+3 | 


where the notation is changed slightly from that of 
Lees. 

In the Newtonian approximation it is assumed that 
the fluid momentum normal to the surface is entirely 
lost upon impact, so that the pressures are proportional 
to the square of the sine of the body-surface slope. 
For slender cones, this results in: 


Cr, y(FRwy)* = 1/2 (10) 


The second-order Newtonian theory has been pro- 
posed by Cole’® to account for the effect of a value of 
y different from one. The formula for the pressure- 


drag coefficient of cones is: 


aL 4 : coe. +(")] 11 
2 Hy + 1) F 


This is plotted in Fig. 2 for y = 1.4. 

It should be noted that all of the approximate theories 
presented above require that the flow be completely 
supersonic—i.e., that the shock be attached to the cone 
and that the Mach number behind the shock be greater 


Cp,y(FRwy)? = 


than one. 

It is apparent in Fig. 2 that, plotted in this form, the 
various approximate theories give considerably differ- 
ent results. The approximate theories are shown again 
in Fig. 3, and the exact calculations for cones made 
by Kopal'' have been plotted on the figure. For moder- 
ate supersonic Mach numbers the data for the more 
slender cones correlate well with the second-order 





Fic. 4. Supersonic pressure-drag parameter (modified) for 
conical noses. Comparison of exact and approximate theories. 
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Fic. 5. Pressure drag of conical noses at WM = 1.0. 


theory. At the higher Mach numbers, however, and 
particularly for the blunter cones, the data diverge 
into a family of curves. In Fig. 4, the ordinate has 
been modified to the sine-squared variation represented 
by Eq. (2), and the calculations of Kopal replotted 
with the new parameters. The improvement is ap- 
parent; the correlation is almost perfect, particularly 
at the higher Mach numbers. The exact results follow 
closely the theoretical curves obtained from the second- 
order and the hypersonic small-disturbance theories, 
and coincide with Newtonian theory at infinite Mach 
number. 

It would appear that the remaining divergence 
shown in the figure, which occurs as the Mach number 
approaches that for shock detachment for a given cone 
angle, cannot be removed by further changes in the 
parameters. In this regard it should be noted that the 
lowest Mach number shown for each cone angle is that 
at which the surface Mach number is equal to one. 
In order to obtain the pressure drag for Mach numbers 
less than this Value, it is necessary to resort to experi- 
mental data or to transonic theory. 

The available data for the pressure-drag coefficients 
of cone-cylinders at a Mach number of one are plotted 
in Fig. 5 versus the nose semivertex angle 6y. The 
points from references 12-14 have been obtained by 
extrapolation from results at higher or lower Mach 
number by use of the principle of stationarity of local 
Mach number in the transonic region (see below). 
The result from reference 15 is a value calculated by 
the relaxation method for a Mach number of one, cor- 
rected for an error as discussed in reference 16. The 
dashed curves in Fig. 5 are the transonic slender-body 
theories of Miles” and Spreiter and Alksne.** The solid 
curve follows the slender-body result of Spreiter and 
Alksne as far as it appears to be applicable, and then 
has been faired through the data and terminated at a 
pressure-drag coefficient of 1.0 for 6y = 90°, which is 
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a reasonable estimate for the force coefficient at sonic 
speed on the end of a circular cylinder in axial flow. 
These results for a Mach number of one can be ex- 
tended into the supersonic range by means of the prin- 
ciple of stationarity of flow near a Mach number of 
one.'* This principle states that the local flow condi- 
tions (Mach number, static pressure) remain constant 
as the free-stream Mach number varies through a value 
of one. Neglecting changes in total pressure, the local 
pressure coefficient at any point on the body is given by 


A 2 fs 
C, = | }- 


~ Mo? 
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It is convenient to convert this expression from a rela- 
tion based on local velocity to one based on local Mach 
number. By use of the compressible-flow relations,“ 
either for isentropic flow or for flow behind a normal 


shock, the local velocity ratio is given in terms of Mach 


(ZY -(4/)4# 
Uo} \Mo/ M?4+5 


Then Eq. (12) becomes (y = 1.4) 


0 [/Me? + 5\"? 
C, = ue ( 4 — | (13) 
7Mn*|\ M? + 5 


The pressure-drag coefficient of an axially symmetric 


number by 


nose of arbitrary shape (based on maximum frontal 
area) is given by 


1 r\? 
Copy = f Coa (7) (14) 


Substituting Eq. (13) into Eq. (14) gives 


10 1/ Mo? + 5\" (5) | ‘ 
— a asi i ~1| (15) 
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For a free-stream Mach number of one, this becomes 


(CDpw)Mo =1 = 


MTL ( 6 , (<) | 
= ——— dj—) —1 (16) 
( 0 \ A? + 0 R JIMo=] 


Now, Eq. (15) can be rewritten and expanded in 


series: 


a CT 6 fe 
Co, » = 7Ma«2\ a i? + 5 + 
(= - ‘\( 6 yo ‘ | ( a) A 
2\M2+5/\M?24+5 aa i f 
ce E © (Me? — 1 | 
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f( ra ) (: ) ) . 
, 7 d — 1? (Id) 
0\M?+ 5 R f 


where terms of order (Je? — 1) have been neglected. 
If it is accepted that the local Mach number is invariant 
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pressure drag of conical noses. 


for free-stream Mach numbers near one, then the in- 

tegral of Eq. (16) may be substituted into Eq. (17), 

giving finally 

7Mo? + 5 C 5 
1 Mew? ( Dpn) Mao = 


Copy sand 


bo 


(Max? — 1) + O[(Mw? — 1)?] (18) 


This result agrees with that of Bryson,’ in that the 
slope of the drag curve at a Mach number of one is 
given by 
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dCp,n 4 2 
Ta = = ~ (Cpgy) 
dMx J/yo=1 Yt! y¥+1 Mo =1 


0 ~o 
=a —— — (Co,,) (19) 
3 6 en eee 
The expansion of Eq. (17) can be extended to the 
second order in terms of (.Z@ — 1) for comparison with 


This gives 


8 


the result of Guderley.* 


6 72 


by 5 125 
Com = [1 (Me — 1) + — (ie — | x 


2 5 65 
Ca a. Chee - Se 


36 
(Meo — 1)2 | +O [(Me — 1)*| (20) 
J 


This expression differs from that of Guderley in the 
coefficient of the third term of the first bracket and in 
the sign of the second term of the second bracket. 

The Mach number region of validity of Eq. (18) is 
uncertain. On the basis of the limited information 
available from transonic theory and experiment, it 
appears that the flow conditions should be approxi- 
mately invariant for values of the transonic similarity 
parameter (8FRy) less than about 1/2. By means 
of Eq. (18), values of the pressure-drag coefficients for 
cones from Fig. 5 have been extrapolated into the super- 
sonic range up to the limiting Mach number defined 
by BFRy = 0.50. The results are plotted in Fig. 6 
in the modified hypersonic similarity form discussed 
above, these extrapolations extending from the left 
axis up to the small tick marks on the curves. Curves 
for the exact calculations of Kopal in Fig. 4 are shown 
as solid lines, and the intermediate portions of the 
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curves have been faired. It can be expected that 
Fig. 6 will give quite accurate values of the pressure- 
drag coefficient for cones of semivertex angles from 5° 
to 45°, and for Mach numbers from one to the order of 
infinity. Of course, at the very high Mach numbers 
it is necessary to consider effects such as the interaction 
of the shock wave and the boundary layer, but these 
considerations are outside the scope of this analysis. 
The calculations of Kopal have been compared many 
times with experimental data for the pressure drag of 
cones with attached shocks. The general validity and 
accuracy of these theoretical results are well established, 
and further comparisons will not be made here. Unfor- 
tunately, there are very few reliable results for cones 
tested in the Mach number range for a detached shock 
wave. A family of curves for the predicted pressure- 
drag coefficient of cones is shown in Fig. 7, along with 
the only test data available (obtained from measured 
pressure distributions) for the detached shock ré- 
gime.'* \4, 7% *? As a matter of interest, the calculated 
curves have been extended into the subsonic region as 
far as the arbitrary limit BFRy = 0.59), with the use of 
the principle of stationarity of flow as given by Eq. (18) 


Ogive Noses 


The general procedures developed for cones will now 
be applied to the determination of the pressure drag of 
ogive noses. Curves for a number of approximate 
theories for pressure drag of ogives are plotted in Fig. 8 
in the standard hypersonic similarity form. In general, 
closed-form analytic solutions are difficult to obtain 
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Fic. 9. Supersonic pressure-drag parameter for ogive noses 
Comparison of exact and approximate theories. 


from any of the approximate theories, and most of the 
results are determined from numerical solutions. All 
of the approximate treatments (except the Newtonian 
flow theory) are based on the assumption of a reason- 
ably slender shape, so the differences between a true 
ogive (circular arc) and an approximate ogive (parabolic 
arc) are negligible. 
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The supersonic linear theory, based on linearized 
source-sink methods, is given by Jones and Margolis? 
for a parabolic-are body of revolution. The super- 
sonic first-order theory is obtained directly from the 
linear results by the usual modification of the approxi- 
mate expression for pressure coefficient, with the result 
that the pressure-drag parameter from first-order theory 
is less than that from linear theory by a constant term 
(1/3). 

The supersonic and hypersonic tangent-cone theories 
have been proposed by Probstein and Bray.”? The 
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Fic. 10. Pressure drag of ogive noses at MW = 1.0. 


tangent-cone approximation is applied to slender para- 
bolic-are noses using the results for cones from the 
supersonic linear theory of von Karman, Eq. (3), and 
the hypersonic approximate theory of Lees, Eq. (9). 
The supersonic first-order result follows in the usual 
way. 

The first-order shock-expansion method was de- 
veloped by Eggers and Savin.?* The shock-expansion 
theory is not limited to slender bodies, but it is specified 
that the hypersonic similarity parameter (8/FRy) 
should be somewhat greater than one. In its complete 
form the theory cannot be expressed in terms of the 
hypersonic similarity parameters above, but a slender- 
body approximation does reduce to hypersonic simi- 
larity form. Both results have been plotted in Fig. 8, 
with the complete theory shown for an average of the 
computations for nose fineness ratios of 3 and 6. 

The second-order shock-expansion theory is given by 
Syvertson and Dennis." By use of the numerical 
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methods presented in the reference, the pressure drag 
of ogives of fineness ratios 3 and 6 has been computed 
for Mach numbers from 3 to 10. Excellent correlation 
was obtained between the two sets of computations 
using the hypersonic similarity parameters, and the re- 
sultant curve is that shown in Fig. 8. 

The hypersonic small-disturbance theory of Van 
Dyke’: © has been given in series form for ogive noses, 
but only the first three coefficients of the series solution 
have been calculated. Van Dyke has concluded that 
this is an insufficient number of terms to give satisfac- 
tory results, and the relative divergence of the hyper- 
sonic small-disturbance curve in Fig. 8 bears out this 
conclusion. 

The slender-body approximation of von Karman’ is 
shown at a Mach number of one, as given by the formula 


Cp y(FRw)? = 7/6 (21) 


Newtonian flow theory has been applied to ogive 
noses at very high Mach Numbers by Young and 
Huntzicker,* Ivey, Klunker, and Bowen,” and Eggers, 
Resnikoff, and Dennis.*! A direct application of New- 
tonian (or impact) theory to ogives results in a formula 
for pressure drag as follows: 

2 (FRy)? + (1/8) 
a ee (22) 
vo 


Co. = 
es [((FRw)? + (1/4) }? 


It has been pointed out by the above authors that the 
effects of centrifugal force around the curved meridian 
line have been neglected in this result. If it is assumed 
that the shock wave lies parallel to the body surface, 
centrifugal force effects are easily included, and it is 
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found that the predicted drag is just one half of that 
given by Eq. (22). However, this application results 
in pressures over parts of the nose that are less than 
zero. If it is then assumed that a vacuum exists in 
this region—i.e., that the shock “‘separates’’ from the 
body—Egq. (22) will be modified to give, approximately, 


1 (FRy)? + (1/8) 


. —— 23) 
2 [(FRy)? + (1/4))? (23) 


Cop _ 
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Eggers, Resnikoff, and Dennis*! have presented a modi- 
fied approach to the inclusion of centrifugal forces in 
which it is not assumed that the shock is parallel to the 
body surface, and they thereby avoid the problem of a 
vacuum and “‘shock separation.’’ For ogive noses, their 
result gives 


2 (FRx)? + (1/8) 
3 [(FRy)? + (1/4) 2? 


[ 1 (FRx)? + (1 > | . 
— (24) 
S (FRy)? + (1/8) 


This is very similar to the results obtained from Young 
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Fic. 12. Comparison of calculated and experimental values of 
the pressure-drag coefficient of ogive noses. 


and Huntzicker’s ‘‘Case 5’’.% For slender bodies this 
reduces to 


7 (FRy)’ + (1/8) 
12 [(FRw)? + (1/4)? 


Copy = (25) 
Eqs. (22), (23), and (25) have been shown in Fig. 8. 
It is not yet clear which of these expressions should be 
used. 

The form of Eqs. (22)-(25) suggests that the hyper- 
sonic correlation factor for ogives should be modified 
to [(FRy)? + (1/4) |?/|(F Ry)? + (1/8)]. This factor 
provides the same property for ogives as the factor 
1/sin? 6y = (FRy)* + 1/4 does for cones, namely, that 
perfect correlation for all fineness ratios is obtained as 
the Mach number becomes large. 

Just as in the case of the cone solutions, the various 
approximate theories in Fig. 8 show an appreciable 
variance in the prediction of pressure drag. In Fig. 9, 
the approximate theories are compared with exact 
calculations by the method of characteristics (from 
published literature) and with calculations made for 
the present report using Van Dyke’s second-order 
theory. These exact (or nearly exact) calculations 
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have been correlated in the figure by using the modified 
hypersonic similarity parameter for ogives which was 
discussed above. A single well-defined correlation curve 
is obtained which does not agree well with any of the 
approximate supersonic theories, but is in good agree- 
ment with the hypersonic second-order shock-expansion 
theory. 

It is somewhat surprising to note that the exact cor- 
relation curve appears to approach the result of slender- 
body theory at the left axis. The exact calculations 
are applicable only if the local Mach number behind 
the attached shock at the nose is supersonic; conse- 
quently, the correlation curve at the far left applies 
only to very slender noses (the last points shown for 
second-order theory were calculated for ogive noses 
of fineness ratio 15). Just as for the conical noses, the 
pressure drag for Mach numbers in the detached-shock 
region must be obtained from experimental data or 
transonic theory. 

Measured _ pressure-drag 
noses”®, *7, 4—87 at a Mach number of one are shown in 


coefficients for ogive 
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Fig. 10. 
experimental points, using the transonic slender-body 
theory of Spreiter®** as a guide for the more slender 
noses. 

The pressure-drag coefficients at a Mach number of 
one, obtained from the faired curve of Fig. 10, have 
been extended into the supersonic region by means of 
the sonic-stationarity principle represented by Eq. (18). 
These results are shown in Fig. 11, where the transonic 
extrapolation extends up to the tick marks on the 
dashed curves. The solid line is the single correlation 
curve through the exact calculations shown in Fig. 9, 
and the dashed curves have been faired into the solid 
curve at values corresponding to the Mach number for 
local sonic velocity behind the attached shock wave. 
Except for the uncertainty regarding the detailed 
values at a Mach number of one, Fig. 11 will give ac- 
curate pressure-drag coefficients for ogive noses of 
fineness ratios from 0.5 (hemisphere) to 10 and for 
Mach numbers from one to the order of infinity. 

Representative curves of pressure-drag coefficient 
versus Mach number are shown in Fig. 12, as obtained 
from the design curves of Fig. 11. 


Other Nose Shapes 


For application to nose shapes other than the cone 
and ogive, the correlations of exact and approximate 
theory have indicated that the second-order super- 
sonic method is reliable for the lower supersonic Mach 
numbers with attached shock and that the second- 
order shock-expansion theory provides valid results at 
the higher Mach numbers. With a few reliable meas- 
urements of pressure distribution in the detached-shock 
region, similar sets of design curves can be prepared 
for other families of nose shapes. 


Boattail 


The pressure drag of isolated boattails has been 
estimated by Fraenkl,** who used slender-body theory. 
In view of the deficiencies of slender-body theory for 
the prediction of the pressure drag of nose shapes at 
Mach numbers much above one, it was considered to 
be desirable to compute the pressure-drag coefficients 
of isolated boattails using the second-order theory. 
(In the present context, “isolated boattail’’ refers to a 
boattail attached to the rear of a cylinder that extends 
forward to infinity.) 

The results of the second-order calculations for coni- 
cal boattails are presented in Fig. 13 for various ratios 
of base diameter to maximum body diameter (some of 
these data were taken from reference 29). The usual 
hypersonic similarity parameters have been used, since 
there appeared to be no theoretical or numerical justi- 
fication for the sine-squared modification. Also shown 
in the figure is the curve for a conical nose (supersonic 
second-order theory) from Fig. 2, the curve representing 
a vacuum over the entire boattail, and Fraenkl’s slen- 
der-body results for several diameter ratios. The slen- 
der-body theory is surprisingly accurate over a large 
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Mach number range. The curves for various diameter 
ratios in Fig. 13 have been faired tentatively into a 
Mach number of one (6/FRgr = 0) but the extrapola- 
tions are only speculative at this time. 

Similar calculations have been made for ogive boat- 
tails and are presented in Fig. 14. It will be noticed 
that the Fraenkl slender-body theory for boattails is 
seriously in error except for very slender bodies near 
a Mach number of one. As in the case of the conical 
boattails, there are no theoretical or experimental re- 
sults for boattails with finite diameter ratios at sonic 
speed, and it is necessary to make tentative extrapola- 
tions of the curves to the left axis. 

The faired curves of Figs. 13 and 14 are presented 
without the calculated points in Figs. 15 and 16, which 
may be considered as design curves for boattail pres- 
sure drag. With regard to the parameters used in the 
figures, it should be noted that the boattail fineness 
ratio is based on the length from the start of the boat- 
tail to the blunt base, and the drag coefficient is based 
on the maximum cross-sectional area of the body. 
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Fic. 15. Supersonic pressure-drag parameter of conical boattails. 


It should be borne in mind that viscous effects will often 
produce large effects on boattail pressures and there- 
fore on pressure drag, and that even the results of 
second-order theory do not always agree well with the 
small amount of experimental data that is available 
on boattails. Whether this disagreement is a result 
of viscous effects, deficiencies in the theory, or errors in 
the test data is not clear. 


Nose-Boattail Interference 


Design curves have been presented for the pressure 
drag of isolated noses and isolated boattails. It re- 
mains to determine the interaction effects for the finite 
body of revolution—i.e., the nese-boattail interference 
drag. This subject has been investigated by Fraenkl” 
by use of slender-body theory. Since the previous re- 
sults have shown that slender-body theory is accurate 
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only near a Mach number of one, further investigation 
by a more nearly exact theory is required. This work 
has not been completed, but a sample of the type of 
curve which is obtained is shown in Fig. 17. The drag 
coefficient ACp,,, represents the additional drag acting 
upon the boattail due to the presence of the nose. The 
data were calculated by second-order theory for an 
ogive-boattail combination with no cylindrical center- 
section, and are presented for several values of the ratio 
of boattail length to nose length (this is an extreme 
case, and the interference drag is larger than might 
normally be encountered). The curves have been 
faired into the values obtained from Fraenkl’s results 
at the left side of the graph. Because of the great 
number of variables involved, many similar figures 
would be required in order to give a complete presenta- 
tion of the results. 


Conclusions 


The summary curves presented in this paper give 
the pressure drag of conical and ogive noses and conical 
and ogive boattails over the complete range of transonic, 
supersonic, and hypersonic Mach numbers. Other 
shapes can be analyzed in the same manner provided 
that an equivalent amount of transonic experimental 
data is available. The analysis has been made with 
the assumption of inviscid flow, so that the effects of 
boundary-layer growth, shock-boundary-layer inter- 
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action and flow separation are not included in the re- 
sults. Regions of separated flow will be present to some 
extent over the rear of most boattail configurations and 
will affect the pressure drag. At the higher supersonic 
Mach numbers, interaction between the nose shock 
wave and the boundary layer will exert an influence 
on the pressure drag of the nose. However, the present 
results provide a reliable basis for evaluating the effects 
of these viscous flow phenomena. 
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Buckling of a Fmite-Length Cylindrical Shell 
Under a Circumferential Band of Pressure’ 


B. O. ALMROTH,* ano D. O. BRUSH** 
Lockheed Missiles and Space Division 


Summary 


This paper is concerned with buckling of a circular cylinder of 
finite length subjected to a symmetrical band of external pressure. 
Both experimental and theoretical results are presented. The 
experimental data were obtained from tests of three thin-walled 
steel cylinders subjected to external pressure by a pneumatic 
tube encircling the test cylinder at mid-length. The theory is 
based on the principle of minimum potential energy, and the 
Rayleigh-Ritz procedure is used to expand the displacement 
components in trigonometric series. 

Theoretical results are given in the form of graphs which show 
buckling pressure as a function of the following ratios: 


Cylinder radius/thickness 
Cylinder length /radius 
Pressure bandwidth/cylinder length 


Theoretical results are in close agreement with existing solu- 
tions to special cases in which (1) the pressure is applied over the 
entire lateral surface, and (2) the pressure is concentrated along a 
circumferential line. The theoretical results are also in agree- 
ment with the test results. 


Symbols 
a = mean radius of undeformed shell 
h = thickness of shell 
8 = length of shell 
b = width of pressure band 
a = is ‘a 
7 = a/h 
6 = b/2a 
E = Young’s modulus 
v = Poisson’s ratio 
p = external pressure 
Per = critical external pressure 
On = membrane strain energy 
Us = bending strain energy 
Q = potential energy of external forces 
J = total potential energy (Um + U, + Q) 
x,2,@ = dimensionless cylindrical coordinates (see 


Fig. 1) 
= strains at a point in the shell’s middle surface 


€x, €o, Yro 


Xz, X¢> Xr@ = Changes of curvature at a point in the shell’s 
middle surface 
u,v, W = dimensionless axial, tangential, and radial 


displacements of a point in the shell’s middle 
surface. The corresponding distances are 
au, av, aw. w is positive inward. 
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t,n, ¢ = virtual increments in u,v,w, respectively 

1p = dimensionless displacement parameters |Eq 
(7)] 

Nn, 2; = dimensionless displacement parameters [Eq 
(7)] 

n, q, p, k = see Eqs. (6), (7), and (8 

K » = see Eq. (12) 

Ji(p, k,n ) 

Jop,k) 7 = see Eq. (10) 

5(p, k) 

r = Pell — v?)/r*E 

Ae = prebuckling deflection parameter 


Subscripts x or ¢ on displacement components denote partial 


differentiation. 


Introduction 


N MISSILE APPLICATIONS, large-diameter cylindrical 
I shells must be designed to resist loads introduced 
by handling-and-storage support frames. It is fre- 
quently undesirable and sometimes impossible to add 
internal stiffeners in the shell at the point of load appli- 
cation. In the absence of stiffeners, compressive hoop 
stresses resulting from the support frame loads can 
cause the structure to fail by buckling. 

This paper presents the results of a theoretical and 
experimental investigation of the buckling of a circular 
cylindrical shell subjected to an axially symmetrical 
band of external pressure. The problem was treated 
theoretically in reference 1. 
ever, the shell was assumed to be infinitely long, and 
as a consequence the stabilizing effect of end rings was 
Shells of finite length are considered 


In that analysis, how- 


not included. 
here. 

The present theoretical analysis is related also to the 
analysis in reference 2. In that analysis, a more 
general loading was considered in that the external 
pressure was not restricted to an axially symmetrical 
band. The present paper explores more fully the re- 
sults of the theory for the symmetrical loading than did 
reference 2, and also includes the results of experimental 
investigation. 

Only a small-displacement buckling theory is con- 
sidered. The shell is assumed to be quite thin (say, 
a‘h = 100) and stresses are assumed to be elastic. 


Experiments 


In the experimental investigation, three thin-walled 
steel cylinders were subjected to pressure applied by 
means of a pneumatic tube encircling the test cylinder 
at mid-length. Prebuckling stresses as well as buckling 
mode and critical pressure were observed. 
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Fic. 1. Coordinate system. 


Test Equipment 

The test cylinders were fabricated from 4150 steel 
sheet formed into cylinders 22.00 in. in diameter and 
36.0 in. long, and welded with one longitudinal seam. 
Bulkheads were tack welded in place 1 in. from either 
end. 

The measured wall-thickness of Cylinder Nos. 1 and 
2 was approximately 0.035 in. The seam weld on these 
two cylinders were ground flush. Cylinder No. 3 was 
fabricated subsequent to completion of tests of the 
first two cylinders. The measured wall-thickness of this 
cylinder was approximately 0.029, and the seam weld 
was not ground. The test setup is illustrated in Fig. 2. 

A pneumatic tube made of natural, pure gum rubber 
was used to apply an axially symmetric band of fluid 
pressure to the outside of the shell at mid-length. The 
width of this pressure band was 5 in. The tube was 
backed up by a restraining frame. A cross section 
through the tube and backup frame is illustrated in 
Fig. 3. Air was used to apply pressure after the tube 





Fic. 2. Test setup. 


TwVULY 1661 


had been filled with water. Forty '/4-in. SR4, con- 
stantan foil, epoxy-backed strain gages were applied 
to Cylinder No. 2 at the locations shown by Fig. 4. 


Test Procedure and Results 


The first test cylinder was not instrumented with 
strain gages. Pressure applied through the pneu- 
matic tube was increased slowly until buckling occurred. 
The first buckle formed in an area adjacent to the weld 
at a pressure of 14.7 psi. Three additional buckles 
formed on the opposite side of the cylinder at 21.5 psi. 
By 30.0 psi, a complete series of ten equally spaced 
buckles had formed. A photograph of the buckled 
shell under a pressure of 30 psi is shown in Fig. 5. 


PLYWOOD 







ALUMINUM 


- BOLT 


PNEUMATIC TUBE 


RUBBER WEDGE 

ry 1/8 IN 

= 5 1N — 
WALL OF TEST CYLINDER 


Fic. 3. Typical cross section through pneumatic tube and 
backup frame 


The second cylinder was instrumented with strain 
gages (Fig. 4). Pressure was increased slowly, and 
strain gage readings were recorded at pressure incre- 
ments of about 5 psi. Recorded circumferential 
membrane strains at test pressures of 5, 10, and 15 
psi are shown in Fig. 8. The first buckle formed near 
the weld seam at 15.5 psi. A second buckle formed 
at 19.8 psi, a third at 21.5 psi and a fourth at 22.4 
psi. By 30 psi, five additional buckles had formed. 
The portion of the cylinder in the region of gages 33-34 
did not buckle. Pressure/strain plots for two pairs 
of back-to-back gages which were located at the centers 
of buckles are shown in Fig. 6. 

The third cylinder, unlike the others, was loaded to 
buckling three times. Between loadings, the applied 
pressure was reduced to zero and the cylinder was 
rotated relative to the test equipment. In the first 
loading, three buckles appeared at 9.8 psi, and a fourth 
at 10.8 psi. One of the first three was in the region of 
the seam weld. In the second loading, three buckles 
again appeared at 9.8 psi. These three, however, were 
in a different portion of the cylinder circumference from 
the first three, and the region of the seam weld was not 
included. In the third load application, the buckles of 
the preceding loading reappeared, together with a 
fourth one, at between 9.0 and 9.1 psi. 

In all of the tests the buckles extended only a small 
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CYLINDRICAL SHELL BUCKLING 


distance outside the pressure band and were of uniform 
size. Each buckle popped in abruptly. The measured 
pressure was the same immediately before and after 
The cylinders continued to accept 


When 


each buckle formed. 
increased pressure in the post-buckling range. 
pressure was removed, the buckles disappeared. 


Theory 


A theoretical solution to the problem is derived 
from the same general theory as was used in references 
1 and 2. 


ciple of minimum potential energy. 


The stability criterion is based on the prin- 
Nonlinear terms 
necessary for a buckling analysis are introduced in the 
basic geometric relationships, and the Rayleigh-Ritz 
procedure is used to reduce the second variation of the 
potential energy to a quadratic form in the virtual 


displacement components. 


The Potential-Energy Expression 

With the notations u,v,w for prebuckling middle- 
surface displacement components, and £,y,¢ for incre- 
mental virtual displacements, the strains and changes 
of curvature in terms of displacements are 


€ = Up t & H+ lot? = pwt & + /ot? 


y% = -9 +s, -f + Ate 
Yzo = 4 i Nr = x(n = &%) 
‘ (1) 
Xz = Wrr — Srr 
~~ = = a C4 = ¢ 
Xzo = eo 


These expressions are based on Eqs. (5) and (6) in 
In deriving them it is assumed that pre- 
accurately 


reference 1. 
buckling displacements are determined 
enough by linear theory; thus, second-order terms con- 
taining “,v,w and their derivatives are omitted. Fur- 
thermore, advantage is taken of the fact that pre- 
buckling displacements are axially symmetrical and 
that axial stresses before buckling are zero. 

The customary expression is used for the total po- 


tential energy of the system: 
V=02+0U,+4+ U, (2) 


where (2 is the potential energy of the external forces, 
and U’,, and LU, are the membrane and bending strain 
energies, respectively. The latter are given by the 


equations 


Ea*h x } 


2(1 — v?) 


i se 
Mle + €,? + 2vere, + _ 


Eath l 
2(1 — vy?) 127° 
f fixe é-)4+ tien +e ~ deen 
(3) 


: ve | dx do, 


For the present analysis, 2 is based on Eq. (10) of 
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Fic. 4. Strain gage locations 


reference 1. Integration by parts with respect to @ 


in that equation yields 
Q = pat J Sf [—w t+ '/2(w? + 2vw, + v*)] dx dp (4) 


where v and w represent total displacement com- 


ponents. Introduction of », w + ¢, and ¢,, respec- 
tively, for vu, w, and w, yields 
2= pe? ffi-wtotialwt ot 


2né, + n*|{dx do (5) 


Insertion of Eqs. (1), (3), and (5) into Eq. (2) yields 
an expression for the total potential energy from which 
the second variation is readily obtained in terms of the 
prebuckling displacement w and the virtual increments 
Ene. 

Trigonometric Expansion of Displacement Components 

The prebuckling displacement w may be determined 
from Eqs. (1) through (4) by use of the principle of 
stationary potential energy. 
metric series, and observing that the vector (£,n,¢)=0 


Expanding w in trigono- 
£ £ 


before buckling, we obtain 


Pu —s : nw 
w= : X, COs 3 (G6) 
whe 2 =1,3,5 a 





Buckled shell, 


Fic. 5. 
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where X — An B, 62V = y r if K (p)A pA s) 
and A, = 4(1 — v?)y(1/n) sin (n76/a) pk = 1,3,5 


B, = 1 — 22 + (1/12y2)(n/a)! In the following, the two differential equations ex- 


pressing equilibrium in the plane of the shell surface 


For the virtual displacement components, the following are used to determine u, and v,. This procedure re- 
forms are assumed : duces the size of the final matrix considerably. The 
source of these equations is discussed in reference 2. 
> : pr } nh 
€ = cos (qg¢) UyA » Sin x They are: 
p=1,3,5 a ” 
pr Vig _ Vere Sra 
7 = sin (9) >» VpAy cos x (7) Vin = (2 = V)Srr¢ + Cog (9) 
p=1,3,5 a > e\ 1. (re 2 of 
| where Vv! [(0?/ Ox") + (07/0?) |? 
> pr \ | - — 
¢ = cos (q¢) ~ A, cos x Phen, with Eq. (9) and the relations 
p= 1,3,5 a 


2 


Ji(p, k, n) = (—1)etktn-D?2 





These expressions are periodic in ¢, and they satisfy 
the forced boundary conditions n = ¢ = UV atx = 4pkn 
+a/2. They also satisfy the natural boundary condi- (pt + kt + n*) — 2[(pk)? + (kn)? + (np)?! 
tions {rr + rfgg = & + v(ng — $) = Vatx = +a//2. , — ore ee 
The parameter q is an integer greater than one; q = 0 (m6/a) + isin (2p76/a)/(2p) if p = k 
and g = | represent respectively axially symmetric sin{(p + k)(p6/a) | 
buckling and buckling of the cylinder as a column with Jx(p, k) = pt+ek 
undistorted cross section. These two cases are not sin|(p — k)(p6/a)].. 
included in the analysis. p—k itp ~k 
Insertion of Eq. (7) into the expression for the total 
potential energy yields for the second variation the 5(p, k) = jl if p=k 
quadratic form \0 ifp#k 
(10) 
and A = p(1 — v?)/wE (11) 





1 2 
e») = 


1 4 2 ; 
cron |(=) + (q* — 1)? + 2 (2) (q? — »|\ — i Jog — vy)(q — %) Dy XrJi(p,k,n) — 


2 b l—v ) 
we have K,, = 6(p, k) {(u, et) + (qvp — 1)? + 2v S wn - Hs (au, + : 
a 


a 


2ry(1 — q\vp = i Vx) = Vptx |J2(p, k) ; (12) 


The solution to the buckling problem is thus reduced 
to finding the eigenvalue, 4, of a matrix of infinite 
order. 


Numerical Analysis Table I 





For numerical analysis the matrices must be reduced COMPARISON OF RESULTS FROM DIFFERENT ANALYSES FOR b/L 1 


to a reasonable size. Accordingly, the sequence of 
































values p and k is truncated as follows: 10° -¥ Pop/E 
a | Fligge | Present Analysis | staore 
Dee 1S, Oe a y | (Ref. 4) | Column A” | Column B”'| (Ref. 5) 
| 
where P is chosen to give sufficient accuracy in the final i - | ee Wa Miata _— 
result. The numerical analysis shows that good accu- Be = tia | = “— saat 
racy is generally obtained with a matrix of rather wa _ | — .. = 1a 1am 
moderate size. 2r 100 157.7 | 158.0 168.6 165.9 
For each g we have one lowest eigenvalue, \,. The wal wai 18.11 | = 18.12 18. 88 18.52 
smallest of all the A, represents the solution of the an 1000 4.610 4.611 4.743 4.664 
problem. r/4 100 1395 1377 1397 1400 
In the general expression for the coefficients of the r/4 400 159.9 159.6 160.4 160. 2 
eigenvalue problem [Eq. (12)], an infinite summation r/4 1000 38.94 38.93 39.06 39.00 
of prebuckling deflection parameters appears. It was 
shown through several numerical examples that, within pe re Ee ER. LER ee ee 
** Column B. Second-order terms of Eq. (4) not included 





the range of parameters investigated, extension of this 
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Fic. 6. Pressure/strain plots 


summation beyond the twelfth term has practically no 
effect on the final result of the computations. Conse- 
sequently, in the numerical analysis the sequence 
n= 1,3,5... was limited ton = 23. 

The most accurate solution to the buckling problem 
is represented by a matrix of infinite order. The com- 
puted buckling pressure decreases with increasing size 
of matrix. By successively increasing P? to a point 
where the buckling pressure is not significantly affected, 
we can determine an acceptable minimum size of the 
matrix. For the same accuracy in the final results, the 
order of the matrix increases with increasing cylinder 
length and decreasing pressure bandwidth, In the 
limiting case b L = 1, the buckling pattern is purely 
sinusoidal in the axial direction (one half wave over the 
entire length). 

A parameter study based upon this formulation has 
been performed. A _ digital computer (Remington 
Rand 1103A) was used for the numerical work. Re- 
sults are given in Figs. 7a—7c for the parametric range 

r/4SgaZ8r 100 < y < 1,000 
0.01 S$ 6/L € 1.0 


For smaller values of b/L, the buckling pressure is 
inversely proportional to the pressure bandwidth. 


Discussion 


Present theoretical results for the prebuckling dis- 
placement w [Eg. (6)] were compared with corre- 
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sponding results based on an analysis in reference 3, 
page 397, and were found to be in close agreement. 
The theories indicate that the compressive hoop stress 
is essentially uniform in the loaded region and that 
it drops to zero at an extremely short distance outside 


the loaded region. 
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The circumferential distribution of the prebuckling °K | a a 
circumferential membrane strain ¢, is shown in Fig. 8 J | L ) 
: ; . ‘agai aa ae ~ 
for the cylinder cross section under the middle of the REF 4L/b=! | | fin — 
pressure band. Theoretical values [based on Eqs. (1) 6 pf fg, 4 
and (6)] and test data are shown. Except for two - at 

. . ° - ° Ww 5 $$$ ttt 
points along the circumference of the cylinder, test > ake 
data and theory are in close agreement. The disagree- & le | —+ 

‘ 4 ++ —__ 


ment at 25-26 may be attributed to the nearness of the 
weld. The second discrepancy (at gages 33-34) was 
at the one region in the cylinder circumference which 
did not buckle (see above). 

In the test of Cylinders No. | and 2, a buckle formed 
in the region of the seam weld at about 15 psi, and three 
additional buckles formed in a region remote from the 
weld at about 21 psi. The theoretical buckling pressure 
for these cylinders is found from Fig. 7 to be 15.9 
psi. The sequence of formation of buckles in these 
two tests suggests that geometric irregularities in the 
welded region may have precipitated buckling in that 
region at a reduced pressure. Furthermore, other 
tests have shown that objects which limit outward 
bulging of the shell during buckling tend to exert a 
stabilizing influence on the shell, and the two rubber 
wedges shown in Fig. 3 partially restrict such bulging 
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Fic. 8. Circumferential membrane strain distribution. 


at the edges of the pressurized zone. Consequently, 
Cylinder No. 3 was fabricated and tested with the 
rubber wedges omitted from the loading device. 

As noted above, the third cylinder buckled twice 
at 9.8 psi, and a third time at about 9.0 psi. Further- 
more, three or more buckles appeared simultaneously 
each time, and the welded region did not participate 
in the initial buckling in the last two load applications. 
The theoretical buckling pressure for Cylinder No. 3 
is found from Fig. 7 to be 10.0 psi. Therefore, the 
theoretical and experimental values for the buckling 
pressure are concluded to be in reasonably close agree- 
ment. 

Results of the theoretical analysis were compared 
with available solutions of related problems. In the 
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Fic. 9. Effect of end-ring restraint for a typical example 


case where the pressure extends over the entire lateral 
surface of the cylinder (6/Z = 1), the solutions obtained 
are in close agreement with those by Fliigge.‘ | The 
solutions are compared in Table 1. Also shown are 
results based on Batdorf’s analysis? (the Donnell 
equations). 

The solution of a second limiting case is available 
in Hahne’s analysis of the buckling of a circular 
cylinder subjected to a radial, circumferential line 
load. An approximate line-load solution is obtained 
from the present analysis if, in a case where the pres- 
sure band is narrow, the buckling pressure is multiplied 
by the pressure bandwidth. For the example a t, 
y = 100, we have the buckling load 4.2 & 1074 Eh 
pounds/inch of circumference. Hahne found this 
value to be 4.6 K 10-7 Eh. 

Out of curiosity, the final equations in the present 
analysis were also developed without the inclusion of 
quadratic terms in (v,w) in the potential energy of the 
external forces [Eq. (4)]. Such terms are frequently 
omitted in cylindrical shell analyses but are known to 
be important in the analysis of a circular ring subjected 
to uniform external pressure (see reference 7). A few 
values for the case 6/L = 1 were computed by use 
of these simplified equations. Results of the calcula- 
tions are shown in Table 1, Column B. It may be seen 
that deletion of these second-order terms increases 
the buckling pressure of the longer shells. 

Fig. 9 shows, for comparison, certain results of the 
analysis in reference | for infinitely long cylinders and 
corresponding results from the present analysis. It is 
seen that the buckling pressure from reference | 1s 
about 20 percent higher than the value shown by the 
present analysis for long cylinders. This difference 
reflects the fact that in applying the Rayleigh-Ritz 
procedure, more degrees of freedom were used in the 
displacement functions of the present analysis than 
were used in the analysis in reference 1. 

The restraining effect of end-rings is also illustrated 
in Fig. 9. The ratio L/b = 1 represents the load case 


(Continued on page 592) 
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Integral-Method Compressible Boundary- 
Layer Calculations for Strong Favorable 
Streamwise Pressure Gradients 

Over an Insulated Surface* 


Martin H. Steiger 

Research Associate, Polytechnic Institute of Brooklyn, 

Freeport, N.Y 

November 17, 1960 

4 es CALCULATIONS of constant-density, laminar-boundary- 
layer properties under strong favorable pressure gradients 

by means of a one-parameter integral method has been treated 

In that paper a discussion was held 

to other 


by Steiger in reference | 


on the relation of the integral-method solution, as given, 


suggested methods of analysis, such as that of local similarity and 


the ‘‘cold-wall’”’ approximation. The breakdown of the usual 


polynomial-profile assumption was also discussed 


This paper presents an extension of the method of reference 1 


to the compressible case for insulated surfaces. It was con 


sidered of interest to compare some of the results with the exist 


ing ‘‘cold-wall”’ estimates given by the method of Lees,? and the 
local similarity.? Experimental evidence* ® 


techniques of 
has shown that the cold-wall approximation is reasonably satis- 


even under conditions for which the ratio of surface to 


is not much less than one 
ent analysis, cold-wall 
thus provides an interesting comparison with the cold-wall analy 
not in the cold-wall category 


factory 
local inviscid temperature The pres- 


which does not utilize the restriction 


sis in a case which is distinctly 
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Surface skin friction. Flow past a circular cylinder. 
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The pertinent comparison is made for the surface shear stress. 

The case of Prandtl number equal to unity is considered here 
to permit the use of the Crocco integral relation, which in the 
pressure gradient case specifies that the stagnation enthalpy is 
constant throughout the flow field. The nomenclature is that 
used by Schlichting.® 

By introducing a Dorodnitsyn transformation, the momentum- 
integral differential equation that governs a two-dimensional, 
compressible boundary layer in the transformed plane is given by 


dvd 4 din pitty? 4. 5° din 1 _" ( pu oufen ( (1) 
dx dx dx uy Piki J u oN w 5» 
where (p/pi)dy = 6, dN (la) 


The boundary conditions are 


P pu Ou/u 6,2, du 
y=0, N=0, u=0, | ] os creas 
= 


pm ON Pwr, dx 

2a) 

y>b, N=1, «=m, OU/ON = O'%Uu/ON? = 0 (2b) 

Let u/u. =1—(1— NY)? (3a) 

Condition (2a) yields 

oO In( py / 

ajia—-1)=Ar+ta pul pus) (3b) 

oN wu 

by? 11 

where A = == : ris (3c) 
Pwri( pu/ pits )w Ax 

while (2b) yields 2<4s5 @ (3d) 


However, sitice the surface is insulated 


I6r 
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Fic. 2. Boundary-layer thickness. Flow past a circular cylin- 
der. Uo = 20,000 ft/sec at 200,000 ft. Nose radius 5 ft. 


h, = h + (u?/2) = constant throughout 1) 

it follows A = a(a — 1) (5a) 
3/5, = a/[((2a + 1a + 1)] > 0 (5b) 

5*/5, = (h,/m){ (8a + 1)/[(2a + 1)(a + 1)]j — 9/6, (5c) 

5/5, = (hs/h) — a/[(2a + 1)(a + 1)](m?/hy) (5d) 


The present theory is shown to agree very closely with the 
exact solution in the stagnation region of a blunt body.’ 

The analysis is applied to the flow about a cylinder. The free- 
stream flight conditions are U’,, = 20,000 fps at 200,000 ft and 
a 5-ft nose radius. The solution is presented in the following 
figures. It should be noted that there were no difficulties with 
regard to singularities in the solution. 

For a further discussion of this analysis see reference 7. 
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On a Generalized Solution for the Torsion of 
a Prismatical Beam of Any Convex Polygonal 
Cross Section 


Howard E. Brandt 
Massachusetts Institute of Technology, Cambridge, Mass. 
September 27, 1960 


oe the author published a paper concerning the 
torsion of a prismatical bar whose cross section is a regular 
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Fic. 1. Geometry of the section. 


\-sided polygon.! This paper is a generalization of the theory 
there presented to any convex polygonal cross section—e.g., 
trapezoidal, rectangular, etc. This generalization is a result 
of the topological character of the stress function in the Prandtl 
sense and is associated with the properties of the Laplacian 


operator. 


MATHEMATICAL FORMULATION OF THE PROBLEM AND ITS 
SOLUTION 
The solution will be found by means of the Prandtl stress 
function, applying the Galerkin method. This method for 
determining stress and displacement fields was outlined by 
Stanisié and Shirely.? 
The Prandtl stress function is of the form 


Vix, y) = , i 
: Lo 
and the Galerkin form of the problem becomes 


ff V2[W(x, y) + 2] n(x, y)dS = 0, vue i2s.....°0) 2 
5 


where »-(x, y) are coordinate functions which must be found, 
such that the boundary condition, ¥(T) = 0, is satisfied on the 
boundary IT of the polygon; A,, are constants satisfying Eq. (2); 
ind S is the domain of integration. 

A function, V(x, y), must first be found which satisfies the 
boundary condition. Refer to Fig. 1. The convex polygonal 
cross section has N sides which divide it into N-triangular 
surfaces, S;, So, ..., Sy-1, Sv. The equations of the respective 
N sides, Lm, must be expressed with respect to a coordinate 
system, so chosen that the point of intersection of the first and 
Nth side of the polygon is a point on the positive x-axis; the 
equations are then expressed as functions of their perpendicular 
distances, pm, from the origin and the angles, 0,, which these 
perpendiculars make with the positive x-axis. It is readily 
shown that 


Lm =X COS Om + ysSin Om — pm = 0, (m = 1,2,3,...,N) (3) 


Thus ¥(x, y) will take the following form: 
N 
V(x, y) = A; IL [x cos 0, + y sin 0, — pn] (4) 
n=l 


where A, is a constant of first approximation to be determined 


later and II denotes the product of expressions with indicated 


index. Eq. (2) becomes 
ey 
SJ [A,V2m(x, vy) + 2]m(x, y)dS = 0 (5) 
S 


v cos 0, + y sin O, — p, (6) 


a 

ll 
—_ 
_ 


where nil 


Solving for A;, one finds 


> ff i %,(x, y)dS 


2=1 


9 m=1 Sm 3 
a ee N \ 
p SJ I] Pr X, VY) Vv? II PAX, V dS 


m=1 Sm n=! n=1 
where the integrals are evaluated over the respective triangles, 
Sm, and 


n(x, Y) = x cos O, + y sin On — pn (8) 
But it is readily shown that 


N N J 
V2 II on(x, y) = IL odnlx, y) \~ >» on (xX, ¥) + 


n=] n= 1 
On, | ; ino, |) 
cos SIN Gn 
— i +1 > (9) 
n=1 Pn(X, y) n=1 On(X, y) f 


The limits of integration for S,, in polar coordinates, r, 6, are 
easily determined. The lower and upper limits with respect 


to r for each S,, are clearly rmz_ = O and rm, = R sec(@ — On], 
respectively. 
The lower and upper limits with respect to @ are (note % = @n 
and po = pwn) 
Pm COS Om—1 — Pm—1 COS Om’ 
Omnz_ = tan='{ - - - 
Pm-1 SIN Om — Pm SIN Om—1 
Pm+i COS Am — Pm COS Omi 
and Omy = tan™ — —— 
Pm SiN O0mi1 — Pm4i SIN Om 


respectively. Eq. (8) in polar coordinates becomes 


¢,(0, r) = rcos (8 — On) — pn (10) 


Hence Eas. (4), (7)-(11) imply 
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n=1 





, mL 
ee 
V(x, y) = Omu R sec (0—0m) N 
yi i) f r| IL ¢4(6. r) 
m=1 OmL 0 n=1 


Evidently the stress and displacement field are uniquely de- 
termined (see Eqs. 14—17).! 
The torsional rigidity is given by 


N Emu R sec (@—@») 
D=26 >> { rW(r,0)drd@ (12) 
0 


m=1 OmL 


where V(r, 0) is obtained from Eq. (11) by substituting x = 
r cos 0, y = r sin6@; and G is the shear modulus. 
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Effect of Thermal Radiation From a Hot Gas 
Layer on Heat of Ablation 


Robert T. Swann 

Aeronautical Research Engineer, NASA Langley Research Center, 
Langley Field, Va. 

October 6, 1960 


SYMBOLS 


C = stagnation-point velocity gradient 
H = enthalpy 

Ha = heat of ablation 

Sa = heat required to inject unit weight into boundary layer 
m = rate of mass injection 

Nu = Nusselt number 

ad = aerodynamic heating rate 

dr = net radiant heating rate to surface 
Re = Reynolds number 

S = Sutherland constant 

T = temperature 

p = density 

bu = coefficient of viscosity 


Subscripts 


e = external flow condition 
O = no injection 
w = wall value 


_— RADIATION from the hot gas layer may constitute an 
appreciable fraction of the total heat input to the stagnation 
area of a vehicle entering a planetary atmosphere from super- 
circular velocity. Unless the gaseous products of ablation are 
opaque to thermal radiation, the ablative mechanism provides 
no shielding against radiant heating, and the heat of ablation will 
be less for combined aerodynamic and radiant heating than for 
aerodynamic heating only.! The purpose of this note is to indi- 
cate the effect of radiant heating on heat of ablation for the special 
case in which the gaseous products of ablation are similar to air. 
Results for other transparent gaseous products may be inferred 
from this case. 

The heat of ablation will be defined as the rate of heat transfer 
to the surface with no injection, divided by the rate of injection, 
1.e., 


Ha = (qo + Qr)/m (1) 


rN Omu R sec (@—@m) N 1 N 
y f f r Il ¢,(8, Lies | II ¢n(x, y) 
ml 6 0 J n=l 
; j N 2 
5 = ; cos #, \~ 
dD on-(4,7) — ( OO -) — 


n=1 


- 2 \ 11) 
so dr do 


n=1 onr(7, 0) n=1 On(7, 8) 





The rate of injection is determined from an energy balance at 
the surface to be as follows: 


m = (q + qr)/Ah (2 


The aerodynamic heating rate is related to the rate of injection 
through the boundary-layer equations.? Exact solutions of the 
boundary-layer equations for a highly cooled axisymmetric 
stagnation point with air injection are shown in Table 1 in terms 
of dimensionless parameters. [The wall temperature has little 
effect on the heat of ablation defined by Eq. (1). | 

The heat of ablation was calculated from Eqs. (1) and (2) and 
Table 1 for a range of values of g,/qo, and the results are plotted 





TABLE I 





Heat Transfer as a Function of Mass Injection Rate 








T/T. = 0.05 S/T, = 0.02 

m/V 2 pwltwC Nu/V Re 
0 0.4096 
0.25 0.2526 
0.50 0.1271 
0.60 0.0882 
0.70 0.0568 
0.80 0.0332 
0.90 0.0175 
1.00 0.0091 





on Fig. 1. With no net radiative heat transfer, the heat of 
ablation is somewhat less than the value predicted by the usual 
linear theory. When radiation to the surface is present, the 
heat of ablation has the limiting value, Ah[{1 + (qo/q,r)], ap- 
proached for large values of (H, — H,)/Ah, corresponding to 
complete blocking of the aerodynamic heating. Note that 
the limiting value of the heat of ablation is independent of the 
properties of the injected gas. Since radiation from a hot gas 
layer will normally be significant only when the enthalpy dif- 
ference across the boundary layer is large, the properties of the 
injected gas are diminished in importance when radiation from a 
hot gas layer is significant. Therefore, the choice of an ablating 

















Fic. 1. Heat of ablation, T7,/7T. = 0.05; S/T. = 0.02. 
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material for such an application may be based on other con- 
siderations. 

Radiation from a hot gas layer to an ablating surface not only 
increases the total heat input to the surface, but also reduces the 
efficiency of ablation in blocking and absorbing the heat input. 
Radiation to an ablating surface is thus a much more serious 
problem than is indicated by the percentage increase in total 
heat input. For example, consider a case with air injection for 


which 
(H, — H,,)/Ah = 10 


If a-/g 0.2, the heating rate is increased by 20 percent, but 
the mass injection rate is increased by 44 percent. If q,/g¢ 
1.0, the heating rate is increased by 100 percent, but the injection 
rate is increased by over 500 percent. 

These increased rates of mass injection will cause large in- 
creases in the required weight of ablating material. However, 
there are several possibilities for reducing the radiant flux to the 
body, such as changes in vehicle configuration and selection of 
ablating materials having favorable values of absorptivity and 
emissivity in both the solid and gaseous states. 
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Transformation of Coordinate Axis Systems 
of Aircraft 
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SYMBOLS 


6 = attitude angle, angle between the x axis of the air- 
plane and its projection on the horizontal 

y = heading angle, angle between the initial heading 
(x axis of the airplane) and the projection of the 
x axis on the horizontal 

rr) = bank angle, angle between the y axis of the airplane 
and a horizontal line in the y-s plane. 

a = angle of attack, angle between the x axis and the 
projection of the wind vector on the plane of 
symmetry 

B = angle of sideslip, angle between the wind vector and 
its projection on the plane of symmetry 

n = angle of attack of longitudinal principal axis 

] = matrix 
}-! = inverse of matrix 

bel = matrix rotation from wind axis system to earth 
axis system 

n’| = matrix rotation from principal axis system to body 
axis system 

a’| = matrix rotation from body axis system to stability 
axis system 

B’) = matrix rotation from stability axis system to wind 


axis system 


X, Y,Z= components along the X, Y, Z axes, respectively 


Subscripts 
e = refers to earth axis system 
pb 


refers to principal axis system 


Il 


= refers to stability axis system 





Aj = cos @ cos 


FORUM 583 
u refers to wind axis system 
B = refers to body axis system 


§ bene are primarily five coordinate axis systems in which an 
airplane’s dynamic motions can be described: earth axis, 
body axis, stability axis, principal axis, and wind axis (see refer 
ence 1 and 2 for a description of these axis systems They are 
related to each other by orthogonal transformations Since 
analysis of the dynamic motions of an airplane in various degrees 
of freedom usually require changes of coordinate systems, this 
note presents a general matrix equation which relates in a single 
equation the five coordinate axis systems of an airplane. The 


general transformation equation can be written in the following 


form 
XxX, xX; 
_? Lw||n'\[a’| |B } l 
Pa bs 


which expresses the transformation from the principal axis system 


to the earth axis system. The first four matrices on the right 


hand side of Eq. (1) are: 


cos #@cos yy singsin#sin Y cos yy cos dsin # 


sin ¥ cos @ sin y sin @ 





L cos @6sin y sin ysin@sin@ sin ¥ cos @ sin 6 2 
cos y cos @ sin y sin @ 
—sin 6 sin @ cos @ cos @ Cos @ 
COS 7 () sin 7 
9’| = 0 ] 0 2a 
sin n 0 COS 79 
COS a 0 sin a 
a’ 0 l 0 2b 
|_ —sin a 0) cos a 
cos 8 sin 6 0 
[8’] = sin 6 cos 8 0 ae 
| 3260 0 1 
and 
1 A, 1 
[L..][n’][a@’][B lo, Ay» 1 } 
1 A 1 
where: 
Ay cos 0, COS Ww COS B COS @ COS 7H 


—COS Oy COS Pe COS B Sin a sin 7 


sin ¢, sin 9, COS y sin 8B COS @ COS 7 
+sin o, sin 6, cos Wy, sin 8 sin @ sin 7 


+sin Yw COS d» Sin B COS a COS 7 


—sin ~» COS dy» sin 8 sin @ sin 7 
— cos Py» COS gy» Sin 0 sin @ COS 7 
—Ccos Py COS d sill 0, COS @ Sin 7 
—sin Py Sin dy» SIN @ COS n 

—sin WY» SiN dy» COS @ sin 7 


sin 8 + sin dy sin 0, cos Py, cos B — 
sin y COS dy» COS B 


Aj3 = COS Oy COS Py COS B COS @ Sin 7 


+ cos 0, COS Py» COS B Sin @ COS n 
—sin dy» Sin 6, COS Px Sin B COS @ Sin 7 
—sin dy» sin 0, COS Y» sin 8 sin @ COs 7 
+ sin Y, COS dy» sin 8 cos @ sin 7 

+ sin Y» COS dy» sin 8 sin @ cos 7 

— COS Pw COS dy Sin 9, Sin @ sin 7 
+ cos Wu cos ¢ sin 0, COS @ COS n 
— sin ¥» sin d» sin @ sin 9 


+ sin Py Sin dy» COS @ COS 7 


Ax = cos 0, sin Wu cos 8 COS a@ COS 7 
—cos 6, sin Ye cos B sin @ sin 7 
—sin YP» Sin O. Sin @» sin B COS @ COS 7 
+sin Wy» sin 0, sin dy sin 8B sin @ sin 7 


—CcOos Pw COS dy» Sin B COS @ COS 7H 
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+ cos Pw COS dw sin B sin @ sin n 
— sin Ww COS dw SiN Oy Sin a COS 7 
— sin Pw COS dw Sin 0, COS @ sin 7 
+ cos Pw Sin dw» sin a Cos 7 
+ cos Pw sin d» COS @ sin 7 


Ao. = COS Oy sin Pw sin B 
+ sin Pw sin Ov Sin dw cos B 
+ cos Pw COS dw cos B 


Av, = COS Oy Sin Pw COS B COS @ sin 7 

+ cos 0» sin Y» Cos B sin @ Cos 7 

— sin Ww» sin 0. sin d» sin B cos @ sin n 
—sin Ww sin 0, sin dy sin B sin a cos n 
—COS Ww COS dy» sin B Cos @ sin n 

—COS Pw COS dw Sin B sin @ cos n 

—sin Pw COS dw Sin 0» sin @ sin 7 

+sin Pw COS dw SiN Oy COS a COS 7 

+ cos Ww sin dw sin @ sin y 

—COS Pw SiN Pw COS a COS 7 


Ax, = —sin @ cos B cos a@ cos 7 
+ sin 0 cos 6 sin a sin 7 
—sin dy COS 0, sin B COS a COS 7 
+ sin dy» COS 0» sin B sin @ sin n 
—COS dy COS I» SiN a@ COS 7 


—COS dw COS Ow COS a Sin 7 


Az. = —sin 0» sin B + sin dw Cos 0, cos B 


Azz = —sin @ cos B cos @ sin n 
—sin 0, cos B sin @ cos 9 
—sin dy» COS Oy sin B cos @ sin 7 
—sin dy COS 6» sin B sin a cos n 
—COS dy» COS Oy Sin a sin n 
4+-COS dw COS Ay COS @ COS n 


The inverse of Eq. (1) is 


Xy Xe 
Yp | = [Lu] [n’] [a] [8] "| Ve 
Z» Ze 


In Eq. (1), 


Xp Xp 
[n’] rotates Y> to Vp 
Lp Zp 
Xp ae 
[a’] rotates Yg to y. 
Zp he 
Xs Xw 
[B’| rotates ¥, to 7. 
Bi Zw 
py [Xe 
and [Ly] rotates Ve to i & 
ws es 


Similarly in Eq. (4), 


as | F Xu 
[Lw]-! rotates E to | 
; Lhe 








E a ie 7 
[B’|-! rotates eo to Ys 
Zz. = 





6 Fa Xp 
[a’]~! rotates E to Vz 
: ae 


Xs ge 
an< [n’]-! rotates Yeg to Ke 


Ze J Lilt. 








1961 


Therefore from Eq. (1) and its inverse, Eq. (4), the transforma- 
tion between any two of the five axis systems considered can be 
obtained by setting the appropriate angles equal to zero. For 
example: 

(1) If we want the coordinates Xx, Ys, Zz in terms of X,, Y,, 
Z., then set 0 = @ = yy = B = n = Oin Eq. (4). 

(2) If we want the coordinates X,, Yu, Z. in terms of Xz, 
Ys, Zz, then set 7 = 0 = ¢@ = y = Oin Eq. (1 
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The Propagation of a Plane Shock Wave Into a 
Moving Fluid* 


Roy M. Gundersen 
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B* THE USE OF a linearization based on small area variations 
the perturbation of a plane shock wave propagating through 
a tube of varying cross-sectional area into an ideal compressible 
fluid at rest was determined by Chester,! and also by the pres- 
ent author?; the former started with the full three-dimensional 
theory, and then used an averaging process, while the latter used 
a one-dimensional analysis. In the present note, these results are 
generalized to the case that the flow in front of the shock is a 
uniform state and thus applicable to an initially stationary shock. 

When the equations, which govern the one-dimensional, un- 
steady flow of an inviscid ideal compressible gas in a tube of vari- 
able cross-sectional area, are linearized in the neighborhood of a 
known (isentropic) solution, the general solution for an initially 
uniform state is: 


A =4%,/2+ a/(7 — 1) = Flx — (uo + cot] + 
co R(x — uot) — uccoE/2E (to. + co) 
B = —wm/2 + a/(y — 1) = Glx — (uo — co)t] + 
co! R(x — uot) — uocoH,/2Eo(uo — co) 
Ho = co? Sizr/Coy(y — 1) = 2 R(x — ult) (1) 


where u, c, s, Eo, E;, and y are the particle velocity, local speed 
of sound, specific entropy, original uniform cross-sectional area, 
area perturbation, and ratio of the specific heats at constant 
pressure c, and at constant volume c,. The subscript zero 
denotes the basic (known) flow, the subscript 1 the terms of 
first-order and F, G, and R are arbitrary functions of one variable 

Let a shock of arbitrary strength be propagating along the 
tube of varying cross section, and let the flow in front of the shock 
be denoted by subscript 1 and the flow behind by subscript 2 
Denote perturbations of a base quantity by a bar. Let P 
denote the pressure, w the shock velocity, v = w — u, m = 
u/c, M = v/c, so that the following relations hold: 


ue = Uy a 2Mia (1 —_ M,~*) ( Y + 1) (2a) 
O1/ce = (y + 1)Mo/M (y — 1 + 2,7?) (2b) 
P./Py = 2yM2/(y + 1) — (¥ — 1)/7 + 1) (2c) 


M2? = (y¥ — 1+ 2 M~?*)/[27 — (7 — 1)M-?3] (2d) 


* This work was supported by a grant from the National Science Founda- 


tion. 





com 
by F 
shor 


in tl 


Sine 


If 
tube 
and 

Fx 


whic 
is al 
as tl 
weak 
T 
Eq 
notil 


Eval 
equi 
G(A) 


(i + 


sforma- 
can be 
) For 


nstonal 
l Load 


Trans- 
lircraft, 


ito a 


tions 
rough 
ssible 
pres- 
‘ional 
used 
[s are 
risa 
10ck, 
, un- 
vari- 
of a 


ially 


da- 





READER 


me = (y + 1) Mom /Mi(y — 1 + 247?) + 


(y¥ + 1)ve = Mial(y — 1 + 244?) (2f) 


From the formulation of the problem—i.e., with the shock 
coming from minus infinity, so to speak, the contribution given 
by F is not included for there is no mechanism downstream of the 


shock which could give rise tosuchaterm. Thus, setting F 0 
in the second of Eqs. (1) gives: 
fi2/2c2 + C2/Cy — 1) — $2/2¢ery(y — 1) 4+ 

mek, /2Eo (1 + m) = 0 (3 


) 


Replacing the second two terms of Eq. (3) by their equivalent 
co from (2b), V/, from 


P,/2yP2, determining f2/c2 from Eq. (2a), « 
» 


2c) gives: 
P, = —(y + 1)(P2 — Pi)ti/Mic(1 — M7?) X 
(1 + 2M. + An?) + (vy + 1001 + MN?) X 
(P2 — Pi)P\/2yP\(1 — M2?) (1 + 2M2 + M,-?) — 
2(P2 — Pi)&/a(l + 2Me + M7?) — (y + 1)(P2— Pi) X 
my MoE\(x)/E (1 mo) (1 — Ms?)(1 + 2Mz M,~? (4) 
Since it follows from (2e) that 
+ 1)meMe/(1 — Meo?) = 24+ (y + 1)m/Mi(1 — M 
Eq. (4) may be written as: 
P. = —(y + 1)(P2 — Pi)t,/Mial(1l — M72) X 
(1 + 2M. + M,~*) + (4 1)(P2 — Pi) X 


(1 + M~*)Pi/2yPi(1 — Me?)(1 + 2Me + M,-*) - 
2(P2 — Py)&/o(1 + 2Meo + M7?) — (y + 1l)m X 
(P2 — P,)E\(x)/Eo(1 + mz2)M(1 — 147?) X 
+ M,-?) — 2(P2 — P))E\(x)/Eo(1 + me) X 
(1 + 2M. + M,~*) (5) 


1 + 2M, 


If the area variations are restricted to a finite length of the 
tube, the ultimate effect is to produce a change in shock strength 
and the shock again becomes uniform 

For the case that the state one is an unperturbed rest state, 
m = 0, M, = w/a, ete., all terms on the right-hand side of 


Eq. (5) vanish except the last so that 


Py = —2(P2 — P)Ei(x)/Eo(1 + me)(1 + 2M2 + M,.-?) 
—K(P2 — P,)E\(x)/Eo 


which is the result first obtained by Chester. The parameter K 
is a monotonic decreasing function of the shock strength defined 
as the pressure ratio P2/P,, and varies between 0.5 for a very 
weak shock to an asymptotic limit of 0.3941 (for 7 = 7/5) 

To complete the present solution, the arbitrary function G of 
Eq. (1) must be determined. This is done most readily by 
noting that the system (1) may be written as: 


~ 


te/c2 + yP2 + meE\(x)/Eo(1 + me) = 0 


—te/co + Pe/yP2 + meE\(x)/Eo(me — 1) = 


2G|[x — (te — c)t] c. 


By addition, it follows that 
P, = —me*(P2 — Py)[2y — (y — 1)Mi~?| Ex(x)/2(m22 — 1) X 
(1 — M,~*)Eo + [2y — (vy — 1)Mi7?7)(P2 — P)G[x - 


> 


(t2 — c2)t]/2ce(1 — M,-?) (6) 
Evaluating this on the shock, x = wt, and replacing P: by its 
equivalent from Eq. (5) gives: 
G(\)/c2 = me[2Me + m1 + M,~*)] Ex(rd)/(me2? — 1) X 
(1 + M,-? + 2M2)Eo — 2(4 + 1)thi[rA, A/ex(1 + Me)] + 
Mial2y — (y — 1)M-7)(1 + Mi? + 2M2) - 
411 — M,~*)&[rA, A/c 1 + Me)|/(2y — 
(y —1)M,~?] (1 + M-? + 2Me2) + (1 + M-*) X 
P,[rd, /cX1 + Me)]/yPi(1l + Mi-? + 2M2) (7) 


FORUM 585 
where 


Thus substituting Eq. (7) into Eq. (6) and noting that 
[2y — (y — 1)M7)/201 — AN?) = (y + :1)/2(1 — M2?) = 


}moMe + (y + 1)m/MN|2y — (y — 1)dNn7) 5 


the complete solution for Pris: 


P./( Ps — P,) = (y + 1)me[2M. + m(1 + WM, 
E\(rd)/2(1 — Me?)(m2? — 1)(1 + M 2Mz)Ey — 
m + 1)E,(x)/2Eo(m2? — 1)(1 — M - 
+ 1) mHIrA, A/oo(1 + M Wail — M x 
(1 + 347? + 2Me) — 2% [rA, A/cex1 + Me)]/a(1 
M,-? + 2M2) + + 1)(1 V P| rd,X/co 1 


M2z)|/2yP(1 — Ms?)1 + M7? + 2M2) (8 


For the case that state one is an unperturbed rest state, all 
terms on the right-hand side of (8) vanish except the first two so 
that 


P./(P; — P,) = 


(1 + M7? + 2Me) — meE\(x)/EoMe(m.? — 1 
— NE\(rA)/Eo — moE\(x)/EopMo me? — 1 


which is Chester’s result He has calculated the values of the 
parameters 7 and N and found that the parameter 7 was a mono 
tonic decreasing function of the shock strength and took, respec 
tively, the values 2, 1, and 0.6076, when the shock was very 
weak, when the flow behind the shock was sonic, and when the 
shock strength became infinite. For subsonic flow behind the 
shock, he found that V increased rapidly from a value of 0.5 for 
a very weak shock and became singular when the flow behind 
the shock was sonic. For supersonic flow, \V became negative 
(singular at ms = 1) and increased monotonically to an asympto 
tic limit of —0.1550 

In the region directly behind the shock, the pressure change 
is given by Eq. (5 When the flow behind the shock is subsonic, 
a disturbance will be reflected downstream of the shock, and the 
pressure perturbation is given by Eq (8 If the area perturba 
tion is such that E;(x) = O for x < O (say), the second term on 
the right-hand side will ultimately become zero 

When the flow behind the shock is supersonic, this disturbance 
will propagate upstream of the shock, and the pressure perturba 


tion is given by the entire Eq. (8 
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Surface Streamlines in Three-Dimensional 
Hypersonic Flows? 
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"T NHE COMPUTATION of boundary layers generally is most con- 
veniently carried out along the surface streamlines of the 


inviscid flow.4s2, This, of course, requires an a priori knowledge of 


+ Based on work performed under the auspices of the U.S. Air Force, 
Ballistic Missile Division, Contract AF 04(647)-269 

* Engineer, Theoretical Aerodynamics. 
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the shape of these streamlines—a knowledge which at best is very 
difficult to obtain by exact mathematical solution even for most 
simple geometrical configurations. The object of this paper is to 
compare the results of three approximate methods of determining 
streamlines to the results from one of the few available exact 
solutions. 

The first of the approximate methods is based on a result ob- 
tained by the Newton-Busemann theory of hypersonic flows.* 4 
This result is that a particle in the shock layer is acted upon only 
by forces normal to the body surface and thus follows a geodesic 
trajectory. Therefore the inviscid flow in contact with the surface 
determines streamlines which are the family of gecdesics originat- 
ing from the stagnation point with assigned initial directions 

In the second method the direction of a streamline is in the 
plane containing the normal to the surface and the direction of 
the free-stream flow. Streamlines determined in this way denote 
the family of paths of steepest descent from the stagnation point 
with respect to the direction of the flow at infinity. This method 
will, therefore, be called the steepest-descent method. Note that 
in the Newton-Busemann theory the initial direction of a particle 
entering the shock layer is along a path of steepest descent. 

The third method utilizes some of the exact flow-field equations 
to determine streamline shapes from a known or assumed pres- 
sure distribution—e.g., the modified Newtonian sine-squared law. 
The relevant equations are the 
velocity to pressure along the isentropic body surface, 


energy equation relating fluid 


1/2u2 + h(p, So) = ho (1) 
the thermodynamic state relation giving density, 
p = p(p, So) (2) 
and the component Euler equation, 
pu*x = —Op/On (3) 


where dy is an elemental distance in the body surface normal to 


the streamlines and «x is the surface component of streamline 


curvature. This method is called the pressure method. 
APPLICATION TO A SWEPT CYLINDER-WEDGE 

Consider a wedge with cylindrical leading edge and swept at an 
angle Q with the free-stream flow. Dimensionless plots of the 
streamlines for two wedge angles are shown for an exact calcula- 
tion as well as for the three approximate methods in Fig. 1. The 
exact solution was determined by superimposing a constant span- 
wise velocity on the velocity distribution obtained from a two- 
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dimensional flow field.® 
Since the stagnation point is at an infinite distance 
about how the 


along the 


leading edge, there is some question geodesic 
method is applied. 
wedge surface the direction of the geodesic and the path of steep- 
est-descent coincide. Since the is developable, the 


geodesic is a straight line on the developed surface of the cylinder 


It can be reasoned that along the planar 


surface 


as well as on the wedge. 
The equation of the path of steepest descent on the cylinder is 


y/D tan 2 = 1/2 In tan (x/D) }) 


where y is the distance along the leading edge, x is the distance 
along the developed surface, and D is the diameter of the cylinder. 
On the wedge of half-angle 8 the path of steepest descent is the 
straight line given by y/tan Q = x/cos p. 

Using modified Newtonian pressures in the pressure method, 
we have 


Pp = po + (fo — fo) cos*(2x/D) 5) 
on the cylinder and an equivalent formula with cos? (2x/D) re- 
placed by sin? 6 on the wedge. The leading-edge pressure pp) and 
are determined by the normal-shock equations 
as functions of po, po,and Uw cos 2, so that p(x/D) is known. 


Since pis a function only of x/D, an integral of Eq. (3) is obtained 


surface entropy S 


as 
§ _ [d(y/D)/d(x/D)|*_ (7 const. 
P (6) 
VW + [d(y/D)/d(x D)j\24 u[p(x/D), 
The conditions that uw = U,, sin Q and d(y/D)/d(x/D) ~ @ as 
x — 0 give the value of the constant in Eq. (6); viz., const. l 


sin 2, and Eq. (6) may be integrated numerically to give y/(D tan 
Q) as a function of x/D. 


CONCLUSION 


method does not give the indicated physical 


Considering for the moment a cone 


The geodesic 
behavior at the leading edge 
attack, geodesics from the 
three-dimensional 


vertex are simply the 
effect, 


three-dimensional effects 


at angle of 
straight generators and there is no 
whereas the other methods do give 
The Newton-Busemann theory in general, leads to many anomal- 
ous results. 

Both the steepest-descent 
essential three-dimensional effects for the swept cylinder-wedge. 
The pressure method agrees more closely with the exact solution 
method 


and pressure methods display the 


than does the steepest-descent method. The pressure 
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Peo = 47 Ib/ft? 
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improves with increasing wedge angle and is very close indeed 
for the 30° wedge. The steepest-descent method increasingly 
diverges from the exact solution with increasing wedge angle. 
The pressure method is somewhat more complicated to apply to 
the blunt wedge than the steepest-descent method, and both are 
correspondingly more complicated for more complex bodies. 
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A Generator of Stable Smoke* 


Donald F. Clausen 

Research Associate, Dept. of Physiology, 
University of Minnesota, Minneapolis 
October 19, 1960 


i WAS DEEMED DESIRABLE to study at low Reynolds numbers 
some of the properties of aerosols made of biological mate- 
rials. Under these conditions, a well-known instability of gas 
flow may be observed such that separation of aerosol particles 
according to size becomes impossible. Since smoke streamers 
are useful as indicators of the state of the flow of a gaseous 
system, smoke from various sources such as burning tobacco or 
tar was injected into the aerosol stream according to suggestions 
offered by members of the Departments of Aeronautical and 
Mechanical Engineering of the University of Minnesota. 
Unfortunately, moisture and tars from these smokes soon 
plugged the metal tubes used for the injection and also covered 
the inside of the apparatus with deposits. Therefore, a smoke 
generator was designed and built to produce quantities of clearly 
visible smoke, free from tars and water, which is obtainable 


from one charge for at least two hours 


DESIGN OF SMOKE GENERATOR 

As shown in Fig. 1, the generator consists of a metal pot made 
of welded, heavy-gauge sheet metal, 10 in. in diameter and 14 
in. high, capped by a 6-in. pipe nipple, of which the cap can be 
tightened by means of a metal rod thrust through pairs of metal 
eyelets welded to it. A rubber gasket provides for a gas- 
tight fit. The upper four inches of the pot is cooled, to protect 
the gasket and to allow for handling, by means of a coil of 1/4-in. 
copper tubing soldered to the pot wall and through which cold 
water passes. The same water flow cools a 24-in. metal con- 
denser through which the smoke passes, and which removes 
water of combustion and some tars, both of which can finally be 
drained from the bottom of the condenser by means of a petcock. 
The partially cleaned smoke is forced through a water scrubber 
several inches deep (not shown) to remove remaining tars and 
then through a second (glass) condenser (also not shown) to re- 
move water introduced by the scrubber. Rubber tubing con- 
ducts smoke around the room to points of use. The smoke 
is injected into the aerosol stream by means of short lengths of 
metal tubing either 0.5 mm or 1.0 mm in. inside diameter. 





* This investigation was supported in part by Public Health Service 
Grant C-3583 from the National Cance: Institute of the Public Health 


Service. 
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Fic. 1. Smoke generator. Entrance of nitrogen and exit of 
smoke are shown by means of arrows 


The pot is heated by a 6-in. gas hot plate. Because it was 
felt that oxygen in the presence of heated combustibles might 
lead to an explosion in the pot, the incoming gas, originally air, 
was changed to tank nitrogen without noticeable effect on the 
vield of smoke. The stream of nitrogen is conducted to the 
bottom of the pot by means of a length of 3/8-in. copper tubing 


run through a packing gland 


OPERATION OF SMOKE GENERATOR 

The pot is packed with about 400 g. of either sanitary wipes 
(Kleenex-type tissues) or surgical sponges (gauze) weighted 
down by means of eight 4-in. spirals of plumber’s tape to prevent 
unburned material from remaining as an elevated dome over the 
ash of burned material. The pot cover is tightened in place, 
water is run through both condensers and the cooling coil around 
the top of the pot, and the gas hot plate is turned on full force 
Smoke is available in quantity within five minutes. By ad- 
justment of the level of the flame of the hot plate and of the ni- 
trogen flow, generous supplies of smoke are obtainable for at 
least two hours. By means of a T-tube in the rubber tubing used 
to conduct smoke to the site of use, excess smoke is shunted into 
a fume hood. No tars have clogged the system over months of 
use, and no deposits form on the inside of equipment. Yet the 
smoke is visible enough to visualize the gas flow easily, particu- 
larly if the room is darkened and if a flashlight is used. Infre- 
quently, very small drops of water will form on the end of the 
metal injection tubing, but these seem to evaporate spontane 


ously in the aerosol stream 


Effect of Stability Augmentation on Frequency 
Response of an Aircraft in Atmospheric 
Turbulence 


Herbert L. Rothman, and Max Chernoff, 

Senior Engineer and Principal Engineer, Respectively, 
Republic Aviation Corp., Farmingdale, N.Y. 

October 14 1961 


SYMBOLS 


w = angular frequency, rad/sec 
C(w) F + iG, Theodorsen function, represents the change 
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| : where 
a i ee A = M,(Zq + Zirkg) — Z.(Mq + Mirkg) 
| ? 
$x 104] B = (Mq + Mirky) + Z Eo= MV (7a 
‘ 
: Anz 4 Eq. (7) is plotted in Fig. 1 for a typical fighter-bomber aircraft, 
a lel keq*.2 \ 0.9M at 35,000 feet and kg = 0.2. Two additional curves are 
5 Saal kg: .4 \ plotted in Fig. 1, for kg = 0.4 and ky 0, in order to show the 
\ effects of doubling the pitch-damper gain and having no sta- 
5.10% bility-augment system. 
It is clear from the three curves in Fig. 1 that stability aug- 
aus? mentation has a significant effect in reducing the dynamic mo- 
tions of an airplane responding to atmospheric gusts 
~5 
a w (rad/sec) REFERENCES 


Fic. 1. Frequency response functions for incremental normal 
load factor for typical fighter-bomber airplane for varying values 
of the pitch-damper gain. 


in static lift forces acting on an airfoil when it is 
performing harmonic oscillations 

F = standard notation for real part of Theodorsen function 

G = standard notation for imaginary part of Theodorsen 
function 

¢(w) = Sears function, represents the change in static lift 
forces due to the sinusoidal gust 

w = vertical component of total airplane velocity, ft/sec 


Wg = vertical component of gust velocity, ft/sec 

0 = body-axis attitude angle, rad. 

q = pitch velocity, y component of body-axis angular 
velocity, 6, rad/sec 

ir = deflection of horizontal tail, rad 

V = total aircraft velocity, ft/sec 

Anz = incremental normal load factor from 1 ‘‘g’’ trim at 
the c.g. 

kg = pitch-damper gain, sec 


= d/dt, time-rate-of-change operator, sec ~! 
|“ = modulus squared 


A N ANALYSIS of the response of an airplane due to vertical at- 
mospheric gusts is presented to show the effects of stability 
augmentation in relieving the dynamic motions of the airplane. 
The longitudinal two-degree-of-freedom (pitch acceleration and 
vertical acceleration ) dynamic equations of motion can be written 
in the following form: 
(Ww — VO) — ZyC(w)w — ZqC(w)q — 
ZizC(w)ir Zuh(w)w, (1) 
6 — MwC(w)w — MqC(w)q — MirC(w)ir = Mud(w)w, (2) 


where the horizontal-tail deflection (77) is a function of the pitch 
velocity (q) and the pitch-damper gain (kg) of the stability- 
augment system—.e., 
ir = kg q (3) 
Therefore Eqs. (1) and (2) can be rewritten as 
(w — V0) — Z.C(w)w — (Zq + Zirkg)C(w)g = Zud(w)w, (4) 
6 — MuC(w)w — (Mq + Mirky)C(w)g = Mud(w)wyg (5) 


From Eqs. (4) and (5), and the auxiliary equation 
1 = 
An, = —- (w — V0) (6) 
g 


the modulus squared incremental normal load factor per unit 
gust velocity can be derived; it is given by the equation: 


Flodin, K., and Sundstrom, M., Jnvestigation of the Response to Rand 
Wind Gusts of a Typical Subsonic Fighter Aircraft in Level Flight with Mach 
Number 0.9 at Sea Level, KTH Aero TN47, Stockholm, 1958. 





An Approximate Solution of the Lifting 
Problem of Thin Airfoils at Sonic Speed* 


Iwao Hosokawa 
National Aeronautical Laboratory, Tokyo, Japan 


March 14, 1961 


SYMBOLS 


Uo, M velocity and Mach number, respectively, of the free stream 
parallel to the x-axis 

Y = ratio of specific heats 

& = perturbation velocity potential divided by U which is 
governed by the equation 

(1 — Mo) Orr + Syy (y + )Mo2br%z7 1) 

¢ linearized potential due to Oswatitsch and Maeder, which 

satisfies 
(l — Mo2)¢grz + gyy = Kez (K: const.) (2) 

a = K/2g2 or —K/2m? 

B? = —-m?>=1-— M 

Ko, Io = modified Bessel functions of the zeroth order 

F(x) = distribution of the cross-sectional area of the airfoil, whose 
chord length is assumed unity 

s(x) = slope of the camber-line of the airfoil 

A = symbol which means the difference of the subsequent 
quantity between the two surfaces of the airfoil 

c* = sonic point of the flow on the symmetric airfoil 

t = thickness ratio of the airfoil 

C7) = incidence angle of the airfoil 

Cp = reduced pressure coefficient 


D” to the recent development of the linearized transonic 
flow theory,':? the analysis of the transonic flow around 
airfoils may be possible. Spreiter’s recently published work, 
particularly, showed results in striking agreement with experi- 
ments in sonic flow case (17,, = 1). The present author, on 
the other hand, devised another method of approach for solving 
the nonlinear transonic-flow equation near the obstacle, making 
use of the results of the linearized transonic-flow theory as a 
starting point.¢ [Eq. (13) is a final formula given by the 
method.] This method has an advantage over Spreiter’s in 
covering the whole transonic range of J/,, and in predicting the 

* The same problem was treated in a similar manner by the author at the 
10th National Congress for Applied Mechanics, Japan 1960. However, 
the previous work should be replaced by the present one, as the former 
included an unnecessary approximate procedure and some erroneous calcu- 
lations were made with it. 


An; |? $(w) |? w*[(AF)? + (wZ» + AG)?] ’ 





| 
eT vs | 
| Wg | i 


: (7) 


[(wBF + G(2AF + E)}? + [w(—w + BG) +(G? — F)A — EF}*f 
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occurrence of a sort of shock wave. Here some results on necessary at first to know ¢ in the lifting case. As is well 
the lifting problem of two-dimensional airfoils at sonic speed known, the doublet method is the most convenient approach to 
will be reported. the linear theory of the thin airfoil with an arbitrary camber. 

In order to develop the author’s method for this problem, it is This can be applied also in the present case although Eq. (2) is 


not a Laplace type; application of the extended Green's theorem 


leads to 


l = = 0 ( Kolaj(x — &)? + B*y?j!? 
g(x,y) = -- ex\*~s) < F(t) + Ag(é) ' = dt,M.2S1 (3) 
2r Jo ' oys B 
1 *x— my . ; Pa) ( Io{'a i (x cals ¢)? — my? {} 2 
oe es J F(z) t Ag(§) P ; &, M, & 1 (4) 
a | oy m 
and particularly to 
1 oh > ) exp [—Ky?/4(x — £)] 
az) += = 7 <4 F’(é) + Ag(é) ( 7 dt, M = ] (5) 
2V 7K Jo { Oy Vx-—é 
rhe first term in the integrand is symmetric with respect to the x-axis; while the second is antisymmetric, and hence Ag should be 
determined so as to fulfil the boundary condition regarding the camber of the airfoil (including the incidence). But, generaliy speaking, 
this is not easy except in the purely subsonic case (K 0). Here the integral sign with a shilling fraction symbol signifies the finite 
part of the Hadamard integral 
However, the exact conversion formula which follows can be found in the case of WV. 1. Then, 
i Oo f ] - o exp|—Ky?/4(x — &)] o “ 
s(x) = lim <— Ag(é) daj;,0OSx*x51 (6) 
0 Oy | 2 «K J 0 oy Vx —€ . 


IIA 


] ; ‘ dé 
= K Ae(t)———_ ,0Z x (7) 
i ‘ (x — §&)3/2 
r J0 


Taking into account the continuity of ¢ at the leading edge: 
Ag(0) = & (8S) 
Eq. (7) is reduced to characteristics of the lifting airfoil. The pressure coefficient 
; sas ii in the lifting case is given by the sum of C, for the symmetric 
K | {— Av'(t)} te i(<s = ) oH airfoil and its deviation due to the camber (+1/2)AC,. (But 
0 


2Nr. 


Fortunately, Eq. (9) is the well known Abel’s integral equation, 


s(x) 
Vx —¢ in order that this treatment of approximation be valid, it is 


necessary that the airfoil have a much smaller camber ratio in 
comparison with the thickness ratio.) 


and so the converse transform can be made immediately as 2 Bocii tueopye ' 
For a symmetric circular-are airfoil at small incidence, we have 


follows: 
: ld © (2V 2/K)s(x) a 4 i oo 1 /x — (1/2)} 6 
— Ay'(£) “= dx, OStS1 (10) al, = C, + 4a = 16) 
w dé J, Vt—x V rl @ Cy t 
Therefore, in which 
2 $ s(x) = qi 8 x 
ee dx, OStS1 11) C, = +1-—, ( ae c* =.1/4) (17) 
VsK Jo Vt-=x 7 3 = 
Here it should be noted that this conversion be unique so and 
that another physical condition such as the Kutta-Joukowsky a? = K%2/{77+ 1)} 16/\/r (a similarity constant) (18 


condition in subsonic flow need not be imposed; but also no 
infinite pressure appears at the trailing edge. Particularly, 


in the case of s(x) = —8@, we have 


Ag(x) = (4/V rK)ovx, OSxE1 (12) 


® Now, in order to formulate the lift distribution AC, in terms 3 
of Ag, we can use the resulting formula which was established | co 
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-4 + 
in the previous work?: 
I ri 
- (y + 1)'3 A 
_ — _ 
o, = (900) * ental 
° 
, ; ° 
(y + 1)! | 2K HAZ 
= pia , {o(x) — ¢(c*)} 0 7 
py ae ay 3 oy 06 08. x/L 
for M, =1 (18) S | 
, fy 
where the double sign corresponds to x S c*. Hence 2 by — = = . 
° 
d : ——- Incidence 6= 0° 
‘s 1 ~ AK j 1K a? ' - 
AC, = ee SS — > = jAg(x) — Ag(c*)s (14) ° ee 6-2 
2 K Ves + 1)13f C, ‘ 
, —— ° ONERA'’s Exp. 
and further, based on the relation between K and c*,* \ fl nl | 
AK/K = }—@er2(c*)Age(c*) + prz(C* )Agr2(c*) } Prr*(c*) (15) ° . . . . 
a Seer ro - siliciaas Fic. 1. Pressure distribution on the two surfaces of a circular-are 


Thus, these results constitute the basic formulas giving the airfoil at incidence 2° 








590 JOURNAL OF THE AEROSP 


as easily derived from reference 4. Eq. (16) shows that the 
lift distribution should be essentially influenced by the thickness 
effect (through a and c) in the case of WM. = 1, as a result of 
the nonlinearity. As an example, a result from these formulas 
[Eqs. (16)-(18)] is shown in Fig. 1, which is in comparatively 
good agreement with the Michel-Sirieix experiment (at ONERA).® 
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Comment on ‘‘Torsion Analysis of 
Multicell Tubes’’ 


W. H. Wittrick* 
Lawrence Hargrave Professor of Aeronautics, 
University of Sydney, Australia 


October 21, 1960 


a has derived, theoretically, equations for the 
shear stresses and for the torsional rigidity of a multicell 
torque tube of rectangular section containing an arbitrary num- 
ber of identical and equally spaced internal webs. This same 
problem was previously solved by the writer? with essentially 
the same results 

It should be pointed out, however, that Niedenfuhr’s Eqs. 
(12) and (13) can be simplified appreciably if it is noted that the 
product of the two roots of Eq, (6) is unity. In fact, the two 
equations |(12)-(13)], reduce to the following: 
' 


5, = Goht;{1 — [(B* + BN*+!-*)/(1 + BY*1)]} 


T = 2Ge@hlt;; N — [28(8% — 1)]/[(8 — 1)(B%*! + 1)]j 
where 6 may be either 8; or Bo. 

The tube considered in reference 2 was slightly more general 
than that of Niedenfuhr in that the two outer webs were assumed 
to have a thickness different from that of the internal webs. 
With this in mind, however, the two equations above are essen- 
tially the same as Eqs. (10) and (11) of the writer’s paper. 


REFERENCES 
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Author’s Reply 


F. W. Niedenfuhr 
Associate Professor, The Ohio State University 
Columbus, Ohio 
November 19, 1960 
T IS A PLEASURE to receive Professor Wittrick’s very useful 
comment on the simplification which can be achieved in the 
equations, [(12) and (13)] as well as to acknowledge his priority 


* At present, Visiting Professor, College of Aeronautics, Cranfield, 
Bletchley, Buckinghamshire, England. 
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in publication. I regret that his earlier work had not come 
to my attention sooner 

Various other modifications of the problem can be made in 
addition to the one he suggests. It is simple, for instance, to 
consider a cross section, like that shown in Fig. 1, in which the 
webs are at an angle 8 to the face sheets. The governing differ- 


ence equation is then 


= 6, — (2 + Lyte/M ct )br+1 + 624 


—2G0A(t./t;) 


where A is the area enclosed by each cell. The rest of the 
analysis then proceeds as in the original example 


On the Maximum Static- Temperature 
Rise in Subsonic Flow 


H. E. Brandmaier 

Chief Project Engineer, Curtiss-Wright Corporation, 
Wright Aeronautical Div., Wood-Ridge, N.J 
November 2, 1960 


| td PROPERTIES resulting from subsonic, one-dimensional 
heat addition to a gas or mixture of gases can be determined 
from a simultaneous solution of the usual mass, momentum, and 
energy equations, the area rate of heat addition, the equation of 
state, and a relation between enthalpy, entropy and two of the 
state variables, say, pressure and temperature. Since enthalpy, 
entropy, and the state equation are generally not explicit func- 
tions of pressure P, and temperature 7, the solution of these 
equations without the aid of a digital computer is generally a 
long, tedious process. However, if only order-of-magnitude 
answers are required, limiting values may suffice. It is the 
purpose of this technical note to show that the maximum static- 
temperature conditions depend only on the initial conditions and 
the equation of state. 

Designating initial conditions by subscript 1 and maximum 
static-temperature conditions by subscript 2, and assuming no 


mass addition, conservation of mass can be written as 
mA,V; = poAeVe2 (1) 


where p is the mass density, A is the one-dimensional flow area, 


and V the corresponding velocity. Assuming an equation of 


state of the form, 
P = f(P,T)pRT (2 
where R is the universal gas constant, Eq. (1) can be written as 


6(T2/T1) = B(P2/P1)(V2/Vi1) (3) 


where 6 = f(P2,72)/f(P1,T;) and 8 = A2/A;. A maximum value 
of 672/T,; can be found by differentiating Eq. (3) with respect 


to V2/V;. Thus, 
(d 6T2/T,)/(d V2/Vi) = B(P2/Pi1) + 

(V2/Vi)(d P2/P1)/(d V2/Vi) (4) 
P./P,; can be found from momentum considerations as follows: 


Conservation of momentum between stations 1 and 2 can be 
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TABLE 1 


Limiting Values of V2/V;, 672/71, P2/P, A2/A1 





V/V P./P, 672/T; A2/A, 
: ] 2 ] 2 l 2 ] 2 
a 
0.0 x t 0.5000 x ] 2 
0.1 5.500 10.50 0.5500 3.025 5.513 1 1.906 
0.2. 3.000 5.500 0.6000 1.800 3.025 1 1.818 
0.3 2.167 3.833 0.6500 1.408 2.204 1 1.739 
0.4 1.750 3.000 0.7000 1.225 1.800 1. 1.667 
0.5 1.500 2.500 0.7500 0.3906 1.125 1.563 1 = 1.600 
0.6 1.333 2.167 O.8000 0.4225 1.067 1.409 1 = 1.538 
0.7 1.214 1.929 0.8500 0.4556 1.032 1.302 1 = 1.481 
0.8 1.125 1.750 0.9000 0.4900 1.0138 1.225 1 1.429 
0.9 1.056 1.611 0.9500 0.5256 1.003 1.168 1. 1.379 
1.0 1.000 1.500 1.0000 0.5625 1.000 1.125 1 = 1.3338 


written as 


en , a 8 
(Vi = Vs) Adp (5) 
g 1 


Assuming a pressure-area relation of the form suggested by 


Cr occ »,1 


Eq. (5) becomes 
(W/g) Vi — V: e(P2A2 — PA;) (7 


1 is equivalent to constant-area flow, and e 0 to 
WV,/gPiAi, Eq. (7) be 


Thus, ¢ 
constant-pressure flow. Letting a 
comes 
B( P2/P,) l — (a e)|( Ve Vi) — 1) (38) 
Substituting Eq. (8) and its derivative with respect to V2/V; 
jnto Eq. (4) yields 
(d 6T2/T;)/(d V2/V; [1 + (a@/e)| — 2(a/e)(V2/Vi) (9) 
Since V2/V; = 1 for e 0 the static-temperature function rise 
is directly proportional to the area ratio 8, and Eq. (9) does not 
apply In addition, « must be positive. 
Setting Eq. (9) equal to zero yields the maximum velocity 
ratio 
Vo/Vi = (e + a)/2a (10) 


Substituting Eq. (10) into (8) yields 


P2/P,; = (a Be) [(« + a)/2a} (11) 
or, after substituting Eq. (6) for 6, 
P2/P, = (a/e)*£[(e + a)/2aIf (12) 


Conversely, the area ratio is found to be 


B = (a/e)' *[(e + a)/2a]'* (13) 


Inserting Eq. (10) into Eq. (3), 
6 (T2/T;) = (a/e)[(e + a)/2a\? (14) 

If 6 can be assumed independent of pressure, these equations 
define conditions for the maximum static-temperature ratio 
For a real gas 6 is the ratio of compressibility factors, 22/21; 
for a mixture of dissociating ideal gases, it is the ratio of mixture 
molecular weights, 1/,/MJ.2, and for an ideal gas or constant- 
composition mixture of ideal gases, it equals 1. 

Knowing ¢ and a, the limiting values of V2/Vi, P2/P1, 672/T;, 
and A»/A, can be found from Eqs. (10), (12), (13), and (14) 

In conclusion, Table 1 gives these limiting values as func- 


tions of a for e = 1 and 2 


REFERENCE 
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Centerline Flow Gradients 
in Axisymmetric Ducts 


Solomon Ciolkowski 

Research Scientist, Grumman Aircraft Engineering Corp 
Bethpage, N.Y 

October 20, 1960 


iy STUDYING nonviscous, compressible-fluid flows through axi 
symmetric ducts of variable cross section, it is often desirable 
to describe the flow behavior along and away from the axis of 
symmetry The gradients of the Mach number and flow angu 
larity are often required to determine the usefulness of these 
ducts. This note is concerned with the analysis of the angu 
larity and Mach-number gradients along the centerline of ayi 
symmetric ducts of variable cross section. Axisymmetric wind 
tunnel nozzles, rocket exhaust nozzles, shock tubes, and shock 


tube wind tunnels are all examples of these ducts 


ANALYSIS OF THE PROBLEM 
The equations describing the isentropic flow of a compressible 


fluid are 


(Op/Ot) + V-(pq) = 0 (1 
p(Oq/dt) + p(q-V)q + Vp = 0 (2) 
p = P p) (3) 


where p is the density, p the pressure, and q the velocity vector 
This note treats steady axisymmetric flows. The coordinate 
system used is formed by the streamlines and the curves normal 
to them with the x-axis as the axis of symmetry Because of 
symmetry, the x-axis isa streamline. Let £ be the distance along 
the streamline, and n the distance normal to the direction of ¢ 


The equations in this coordinate system become (see reference 


l 
q(Op/0E) = —pl(O0g/E) + q(08/On + sin 6/r) (4 
q(og/0E) = —(1/p)(Op/ dé (5) 
g*(00/0E) = —(1/p)(Op/dn (6 
P= Fy (7 


where @ is the angle between the streamline and the x-axis, g is 
the absolute velocity, and ¢ is the distance from the x-axis 


On the centerline—i.e., the x-axis—the equations become 


g(dp/dx) = —p[(dg/dx) + 2q(d0/dr) (8 
q\ Og ox) = —-(1 p (Op Ox = —(a’ p)(Op/ Ox (9 
gd0/dx) = —(1/p)(dp/dr) (10) 


Combining Eqs. (8) and (9), 
[1 — (1 M?)|(Op/0x) = —2 (00/0) (11 
Thus a density gradient along the centerline will result in an 
angularity gradient normal to the centerline and vice versa 
Since on the centerline 6 = 0, Eq. (10) results in the relation 
Op or = 0 (12) 
Thus, the pressure gradient normal to the centerline is zero, and 
therefore the Mach-number gradient normal to the centerline is 
Zero 
Ducts OF VARYING CROSS SECTION, 
ONE-DIMENSIONAL APPROXIMATION 
Making use of Eq. (9), Eq. (11) may be written in the form 
2 


(.M? — 1)(1/q)(0g/Ox) = 2(06/dr) (13) 


If the flow is analyzed using one-dimensional theory to relate 
area change dA /dx to change in flow velocity g, it is found that 


(1/A)(dA/dx) = —(1 — M?)(1/q)(0q/dx) (14) 
Accordingly, after substitution of Eq. (14) into (13), 
00/dr = (1/2)(1/A )(0A /Ox) (15) 
Since A = 7rR?, with R equal to the radius of the duct, 
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00/dr = R,/R = tan Owan/R (16) 
where Owail is the local angle between the walls of the duct and 
the x-axis. Eq. (14) results in 

dM l q re dR 
dx (M2—1)(dq/0M) R dx 


CONICAL NOZZLE 
For the case of a conical nozzle, Eq. (16) becomes simply 
00/or = 1/L (18) 

where L is the length of the cone from its generating apex. Thus, 
the angularity gradient is a function only of the length of the 
conical nozzle. 

For a perfect gas, the Mach-number gradient can be obtained 
by using the relationship 


7’ => 
(4 


and Eq. (17). It follows that 


dq l l 
dM q (ly — 1)/2 + (1/M?)|? M3 


so that using Eq. (17) 


adM/dx = 2|M*/(M? — 1)|[(y — 1)/2 + 
(1/M?)|(1/R)(dR/dx) (20) 


so that for large MJ one obtains the approximate relationship 
dM/dx ~ (y — 1)M(1/R)dR/dx 


Accordingly, dM/dx ~ (y — 1)(M/L) 


In order to provide some feeling for the magnitude of the 
angularity and Mach-number gradients, consider a_ conical 
nozzle with exit diameter of 1.5 ft and 12.5° semiangle with exit 
M = 10, y = 1.4. Neglecting the boundary layer, this nozzle 
will have an angularity gradient at the centerline of 06/dor = 
1.41° per in. The Mach-number gradient is approximately 
dM/dx = 0.099 per in 


Buckling of Cylindrical Shell 


REFERENCE 
1 Chu, Boa-Teh, Wave Propagation and the Method of Characteristics in 
Reacting Gas Mixtures with Applications lo Hypersonic Flow, WADC TN- 
57-213, ASTIA Doc. AD118350, 1957 


Addendum to ‘‘The Characteristics of 
Two-Dimensional Sails in Hypersonic Flow’’ 


Walter Daskin and Lewis Feldman 

Scientific Supervisors, General Applied Science Laboratories, Inc., 
Westbury, L.I., N.Y. 

February 16, 1961 


— AUTHORS would like to point out that the steady shape 
assumed by a flexible member under the sine-squared law 
of forces which we displayed in reference 1 has been previously 
published. The hypersonic sail! is equivalent to a towing 
cable, a balloon cable, a kite string, or a paravane cable. Pub- 
lished work on the shapes assumed by these devices and the 
forces acting on them appeared before either of the authors were 
born. It is our pleasure to acknowledge that references 2-5 
contain the essentials of reference 1; some of these go further 
to include body forces, oscillations about the equilibrium posi- 
tion, and friction. These venerable tomes, lightly sprinkled with 
the dust of decades, wrestled successfully with the intricate prob- 


lems of hypersonic aerodynamics of the space age 


REFERENCES 


The Characteristics of Two-Dimensional 
25, 


! Daskin, W., and Feldman, L., 
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Comm. Aero., R&M 554, October 1918 

3 Glauert, H., The Stability of a Body Towed by a Light Wire, Brit. Ad- 
visory Comm. Aero., R&M 1312 (Ae 451), February 1930 

4 Glauert, H., The Form of a Heavy Flexible Cable Used for Towing a 
Heavy Body Below an Aeroplane, Brit. Advisory Comm. Aero., R&M 1592 
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Landweber, L., and Protter, M. H., The Shape and Tension of a Light 
Flexible Cable in a Uniform Current, J. Appl. Mech., Vol. 14, No. 2, pp. 


Al21-—126, June 1942 
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(Continued from page 578) 





in which the pressure is applied over the entire lateral 


surface. It may be seen that the end-rings in this 
example have little effect if the cylinder length is more 


than about five times the bandwidth. 
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